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Creating  a  virtual  shock  environment  by  use  of  a  computer  to  model  the  ship 
structure  and  the  surrounding  fluid  presents  a  valuable  design  tool  and  an  attractive 
alternative  to  these  tests.  The  research  summarized  in  this  report  investigated  the 
accuracy  of  shock  simulation  using  the  shock  trials  conducted  on  USS  WINSTON  S. 
CHURCHILL  (DDG  81)  in  2001.  All  three  explosions  DDG  81  was  subjected  to  are 
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of  fluid  volume  size,  mesh  density,  mesh  quality,  and  shot  location  are  investigated. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  BACKGROUND . 1 

B.  SCOPE  OF  RESEARCH . 3 

II.  UNDERWATER  SHOCK  MODELING  AND  SIMULATION . 5 

A.  STRUCTURAL  MODEL . 5 

1.  USS  WINSTON  S.  CHURCHILL  (DDG  81) . 5 

2.  Ship  System  Damping . 8 

B.  SIMULATION  CODE . 11 

1.  Model  Generation  and  Pre-Processing . 11 

2.  Simulation  Processing . 13 

a.  FLUMAS . 15 

b.  AUGMAT . 15 

c.  TIMINT . 16 

3.  Simulation  Post  Processing  and  Data  Extraction . 16 

4.  Data  Processing  and  Comparison . 16 

III.  FLUID  MESH  MODELING . 19 

A.  FLUID  MESH  CONSTRUCTION . 19 

1.  Inner  Liner . 19 

a.  USA  Code  Stability  Criteria . 20 

b.  Orthogonality . 21 

c.  Waterline . 21 

2.  Inner  Mesh . 22 

a.  Thomas-Middlecoff  Relaxation . 22 

b.  Equipotential  Relaxation . 23 

3.  Transitional  Mesh . 24 

4.  Outer  Mesh . 25 

B.  FLUID  MODEL  PROCESSING . 27 

IV.  DATA  COLLECTION  AND  ANALYSIS . 29 

A.  VELOCITY  RESPONSE  “DRIFT” . 29 

B.  HIGH  FREQUENCY  “NOISE” . 35 

V.  SIMULATION  RESULTS  AND  ANALYSIS . 37 

A.  ANALYSIS  AND  COMPARISON . 37 

1.  Sensor/Node  Location . 37 

2.  Error  Measurement . 38 

B.  SEVERE  SHOCK  (SHOT  3) . 41 

1.  Surrounding  Fluid  Depth  Study . 42 

a.  Two-Layer  Inner  Liner  Fluid  Mesh . 42 

b.  One-Half  Cavitation  Depth  Fluid  Mesh . 44 

c.  Cavitation  Depth  Fluid  Mesh . 45 


v 


d.  2X-Cavitation  Depth  Fluid  Mesh . 47 

2.  High  Pass  Filtered  Data  vs.  Drift  Compensated  Data . 49 

3.  Fluid  Mesh  Depth  Sensitivity . 50 

2.  Fluid  Mesh  Quality . 52 

3.  Fluid  Mesh  Density . 57 

a.  Error  Comparison . 59 

b.  Velocity  Plots . 63 

C.  SHOT  1 . 67 

1.  Fluid  Mesh  Depth . 67 

2.  Mesh  Quality . 68 

3.  Mesh  Density . 69 

a.  Error  Comparison . 70 

b.  Velocity  Plots . 75 

D.  SHOT  2 . 79 

1.  Error  Comparison . 80 

2.  Velocity  Plots . 84 

E.  QUALITATIVE  COMPARISON . 88 

F.  EXTENDED  DEPTH  FLUID  MODEL . 90 

G.  LATERAL  AND  LONGITUDINAL  VELOCITY  COMPARISON . 95 

1.  Lateral  Velocity  Response . 96 

2.  Longitudinal  Velocity  Response . 101 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS . 106 

APPENDIX  A.  DDG  81  VELOCITY  RESPONSE  PLOTS . 110 

A.  SHOT  1 . 110 

B.  SHOT  2 . 115 

C.  SHOT  3 . 121 

1.  Vertical  Velocity  Plots . 121 

2.  Lateral  (Athwartship)  Velocity  Plots . 126 

APPENDIX  B:  FLUID  MODELING  IN  TRUEGRID . 130 

A.  STRUCTURAL  MODEL . 130 

B.  INNER  LINER . 132 

C.  INNER  MESH . 139 

D.  TRANSITIONAL  MESH . 146 

E.  OUTER  MESH . 149 

F.  DAA  BOUNDARY . 151 

APPENDIX  C:  TRUEGRID  MODELING  COMMANDS . 154 

APPENDIX  D:  USA/LS-DYNA  INPUT  DECKS . 162 

A.  LS-DYNA  KEYWORD  FILE . 162 

B.  USA  INPUT  FILES . 164 

1.  FLUMAS . 164 

2.  AUGMAT . 164 

3.  TIMINT . 165 


LIST  OF  REFERENCES . 166 


INITIAL  DISTRIBUTION  LIST 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


LIST  OF  FIGURES 


Figure  1.  USS  WINSTON  S.  CHURCHILL  (DDG  81)  FINITE  ELEMENT 

MODEL  [from  Ref.  8] . 5 

Figure  2.  Equipment  Modeled  as  Rigid  Bodies  in  DDG  81  FEM  [from  Ref.  8] . 7 

Figure  3.  Vertical  Modal  Damping  Ratio  for  Single  Area  Group  of  DDG  53  [from 

Ref.  13] . 9 

Figure  4.  Damping  Ratio  Curve  Plot  for  Using  Recommended  Rayleigh  Damping 

Coefficients  for  DDG  53  [from  Ref.  13] . 11 

Figure  5.  Model  Generation  and  Simulation  Flow  Chart . 12 

Figure  6.  Fluid  Model  Two  Layer  Inner  Liner . 22 

Figure  7.  Fluid  Model  Inner  Mesh . 24 

Figure  8.  Fluid  Model  Transitional  Mesh . 25 

Figure  9.  Fluid  Model  Outer  Mesh . 26 

Figure  10.  Complete  Coupled  Fluid-Structure  Model  (Shot  1) . 27 

Figure  11.  Raw  Velocity  Data  Showing  Drift . 30 

Figure  12.  Velocity  Data  After  Drift  Compensation  Has  Been  Applied . 30 

Figure  13.  Comparison  of  High  Pass  Filtering  vs.  Drift  Compensation . 31 

Figure  14.  Two-Pole  High-Pass  Bessel  Filter  With  3Hz  Cut-off  Frequency . 32 

Figure  15.  Fourier  Spectra  Plot  Showing  Structural  Low  Frequency  Response 

(Vertical  Velocity  Input) . 33 

Figure  16.  Butterworth  vs.  Bessel  High  Pass  Filter . 33 

Figure  17.  Velocity  Response  After  Filtering  0.5Hz  High-Pass  Butterworth  Filter 

(Third-Order) . 34 

Figure  18.  Effect  of  Low-Pass  Filtering  Data . 35 

Figure  19.  Shot  3  Geometry.  Most  Severe  Shot,  Charge  Off  Port  Midships . 42 

Figure  20.  Fluid  Model  Two  Layer  Inner  Liner . 43 

Figure  2 1 .  One-Half  Cavitation  Depth  Fluid  Mesh  [from  Ref.  26] . 44 

Figure  22.  Cavitation  Depth  Fluid  Mesh  [from  Ref.  26] . 46 

Figure  23.  2X-Cavitation  Depth  Fluid  Mesh  [from  Ref.  26] . 48 

Figure  24.  2nd  Generation  Fluid  and  Structural  Model . 5 1 

Figure  25.  Sonar  Dome  Region  Showing  Improvement  in  Mesh  Unifonnity  of  3rd 

Generation  Fluid  Model . 53 

Figure  26.  Keel  Region  Showing  Improvement  in  Mesh  Uniformity  of  3rd  Generation 

Fluid  Model . 53 

Figure  27.  Comparison  of  2nd  and  3rd  Generation  Fluid  Model  Errors  in  Relation  to 

Longitudinal  Position  on  the  Structural  Model . 55 

Figure  28.  Sonar  Dome  Region  Velocity  Plot  Comparison . 56 

Figure  29.  Midship  Region  Velocity  Plot  Comparison . 56 

Figure  30.  Stem  Region  Velocity  Plot  Comparison . 57 

Figure  3 1 .  Natural  Frequency  Error  of  FEM  Pin-Pin  Beam  as  a  Function  of  Number 

of  Elements  Used . 58 

Figure  32.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison . 59 


IX 


Figure  33.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison . 59 

Figure  34.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison . 60 

Figure  35.  Russell’s  Comprehensive  Error  Trend  as  a  Function  of  Longitudinal 

Position . 62 

Figure  36.  Keel  Sensor  V2000V . 63 

Figure  37.  Keel  Sensor  V2016V . 64 

Figure  38.  Keel  Sensor  V2020V . 65 

Figure  39.  Bulkhead  Sensor  V2 124V . 66 

Figure  40.  Mast  Sensor  A2240V . 66 

Figure  41.  Shot  1  Geometry.  Most  Extreme  Aspect . 67 

Figure  42.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison . 70 

Figure  43.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison . 70 

Figure  44.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison . 7 1 

Figure  45.  Shot  1  Russell’s  Comprehensive  Error  Trends . 74 

Figure  46.  Keel  Sensor  V2000V . 76 

Figure  47.  Bulkhead  Sensor  V2 124V . 76 

Figure  48.  Keel  Sensor  V20 16V . 77 

Figure  49.  Keel  Sensor  V2020V . 78 

Figure  50.  Mast  Sensor  A2240V . 78 

Figure  5 1 .  Shot  2  Geometry.  Charge  Off  Starboard  Midships . 79 

Figure  52.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison . 80 

Figure  53.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison . 8 1 

Figure  54.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison . 8 1 

Figure  55.  Russell’s  Comprehensive  Error  Trend  as  a  Function  of  Longitudinal 

Position . 84 

Figure  56.  Keel  Sensor  V2000V . 85 

Figure  57.  Bulkhead  Sensor  V2 124V . 86 

Figure  58.  Keel  Sensor  V2016V . 86 

Figure  59.  Keel  Sensor  V2020V . 87 

Figure  60.  Mast  Sensor  A2237V . 88 

Figure  6 1 .  Average  Russell’s  Comprehensive  Error  Factor  Vertical  Trend . 89 

F igure  62 .  Bulkhead  S  ensor  A2 104V . 94 

Figure  63.  Keel  Sensor  V2034V . 95 

Figure  64.  Keel  Sensor  A2001A . 98 

Figure  65.  Keel  Sensor  A20 15 AI . 98 

Figure  66.  Bulkhead  Sensor  A23 1 1 A . 99 

Figure  67.  Keel  Sensor  A2021A . 100 

Figure  68.  Bulkhead  Sensor  A2105AI . 100 

Figure  69.  Mast  Sensor  A223 8 AI . 101 

Figure  70.  Bulkhead  Sensor  A2103F . 102 

Figure  7 1 .  Bulkhead  Sensor  A2 1 1  IF . 103 

Figure  72.  Bulkhead  Sensor  A2 106F . 104 

Figure  73.  Bulkhead  Sensor  A21 18F . 104 

Figure  74.  Bulkhead  Sensor  A2242F . 105 

Figure  75.  Bulkhead  Sensor  A2239F . 105 


x 


Figure  76.  Keel  Sensor  V2002V:  (RM=0.0772,  REM).  1755,  RC=0.1699) . 110 

Figure  77.  Keel  Sensor  V2007V:  (RM=-0.1 178,  RP=0. 1760,  RC=0. 1877) . Ill 

Figure  78.  Bulkhead  Sensor  A2109V:  (RM=-0.0679,  RP=0.0848,  R00.0962) . 11 1 

Figure  79.  Bulkhead  Sensor  V2013V:  (RM=-0.1366,  RP=0.1819,  RC=0.2016) . 112 

Figure  80.  Bulkhead  Sensor  A2104V:  (RM=-0.2235,  RP=0.1 184,  RC=0.2241) . 112 

Figure  81.  Bulkhead  Sensor  A21 16V:  (RM=-0.0679,  RP=0.0848,  RC=0.0962) . 113 

Figure  82.  Bulkhead  Sensor  A2237V  (RM=-0.0906,  REM).  1203,  RC=0.1334) . 113 

Figure  83.  Keel  Sensor  V2034V:  (RM=0.0296,  REM.2233,  RC=0. 1996) . 1 14 

Figure  84.  Keel  Sensor  V2019V:  (RMM.1463,  RP=0.2447,  RC=0.2526) . 114 

Figure  85.  Keel  Sensor  A2413V:  (RM=0.0171,  RP=0.3649,  RC=0.3238) . 115 

Figure  86.  Keel  Sensor  V2002V:  (RM=0.0579,  RP=0.2001,  RC=0.1846) . 116 

Figure  87.  Keel  Sensor  V2007V:  (RMM.0245,  RP=0.2092,  RC=0.1866) . 116 

Figure  88.  Bulkhead  Sensor  A2 109V:  (RAM.  1170,  RP=0.1647,  RC=0.1790) . 117 

Figure  89.  Bulkhead  Sensor  V2013V:  (RMM.1049,  RP=0.1420,  RC=0.1565) . 117 

Figure  90.  Bulkhead  Sensor  A2104V:  (RMM.0388,  RP=0.0953,  RC=0.0912) . 118 

Figure  91.  Bulkhead  Sensor  A21 16V:  (RM=0.0070,  RP=0.1003,  RC=0.0891) . 118 

Figure  92.  Mast  Sensor  A2240V:  (RMM.0626,  RP=0.1335,  RC=0.1307) . 119 

Figure  93.  Keel  Sensor  V2034V:  (RMM.1442,  RP=0.1874,  RC=0.2095) . 119 

Figure  94.  Keel  Sensor  V2019V:  (RMM.1327,  RP=0.2391,  RC=0.2423) . 120 

Figure  95.  Keel  Sensor  A2413:  (RMM.1245,  RP=0.2219,  RC=0.2255) . 120 

Figure  96.  Keel  Sensor  V2002V:  (RM=0.0762,  RP=0.2303,  RC=0.2150) . 121 

Figure  97.  Keel  Sensor  V2007V:  (RMM.1198,  RP=0.2883,  RC=0.2766) . 122 

Figure  98.  Bulkhead  Sensor  A2109V:  (RM=0.1353,  RP=0.1846,  RC=0.2028) . 122 

Figure  99.  Bulkhead  Sensor  V2013V:  (RMM.0342,  RP=0.1970,  RC=0.1772) . 123 

Figure  100.  Bulkhead  Sensor  A2104V:  (RM=0.0207,  RP=0.1760,  RC=0.1571) . 123 

Figure  101.  Bulkhead  Sensor:  (RM=0.0353,  RP=0.1789,  RC=0.1616) . 124 

Figure  102.  Mast  Sensor  A2237:  (RM=0.0613,  RP=0.1877,  RC=0.1750) . 124 

Figure  103.  Keel  Sensor  V2034V:  (RMM.0966,  RP=0.1537,  RC=0.1608) . 125 

Figure  104.  Keel  Sensor  V2019V:  (RMM.1212,  RP=0.1830,  RC=0.1945) . 125 

Figure  105.  Bulkhead  Sensor  A2413V:  (RMM.1315,  RP=0.2150,  RC=0.2233) . 126 

Figure  106.  Bulkhead  Sensor  A2102A . 126 

Figure  107.  Bulkhead  Sensor  A21 10A . 127 

Figure  108.  Bulkhead  Sensor  A21 17AI . 127 

Figure  109.  Mast  Sensor  A2241 A . 128 

Figure  1 10.  Keel  Sensor  A2033A . 128 

Figure  111.  Example  of  Curves  Used  for  Construction  of  Inner  Fluid  Mesh  Layer . 132 

Figure  1 12.  Face  Sets  Used  to  Generate  Inner  Fluid  Mesh  Layer  [from  Ref.  21] . 133 

Figure  113.  Inner  Mesh:  Lower  Port  and  Starboard  Hull,  Sonar  Dome,  and  Seams . 134 

Figure  1 14.  Transition  to  Test  Waterline:  Aft  Section  of  Upper  Inner  Liner . 136 

Figure  115.  Transition  to  Test  Waterline:  Forward  Section  of  Upper  Inner  Fluid  Mesh..  137 

Figure  116.  Completed  Inner  Mesh  Layer  of  Inner  Liner . 137 

Figure  117.  Completed  Outer  Mesh  Layer  of  Inner  Liner . 139 

Figure  118.  Fluid  Model  Projection  Illustration  Showing  “Shadow  Regions”  Created 

by  the  Contours  of  the  Ship  Model . 140 

Figure  119.  Curves  and  Surfaces  Used  to  Fonn  Outer  Boundary  of  Inner  Mesh . 142 


xi 


Figure  120. 
Figure  121. 
Figure  122. 
Figure  123. 
Figure  124. 


Completed  Inner  Mesh . 146 

Curves  and  Surfaces  Used  to  Form  Outer  Boundary  of  Transitional  Mesh  ..147 

Completed  Transitional  Fluid  Mesh . 148 

Completed  Outer  Mesh . 151 

Coupled  Fluid-Structural  Model . 152 


xii 


LIST  OF  TABLES 


Table  1.  DDG  81  Finite  Element  Model  Specifications  [from  Ref.  8] . 7 

Table  2.  Weighted  Mean  of  a  (DDG  53)  [from  Ref.  13] . 10 

Table  3.  Weighted  Mean  of  (5  (DDG  53)  [from  Ref.  13] . 10 

Table  4.  DDG  81  Wetted  Surface  Shell  Element  Properties . 20 

Table  5.  Statistical  Study  of  Low-Pass  Filtered  vs.  Unfiltered  Simulation  Data . 36 

Table  6.  Vertical  Velocity  Response  Sensor  Location  (N/A  =  no  data  available) . 38 

Table  7.  Russell’s  Comprehensive  Error  Factor  Acceptance  Criteria . 41 

Table  8.  Fluid  Mesh  Liner  120msec  Russell’s  Error  Comparison  [from  Ref.  26] . 43 

Table  9.  One-Half  Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error 

Comparison  [from  Ref.  26] . 45 

Table  10.  Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error  Comparison  [from 

Ref.  26] . 47 

Table  11.  2X-Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error  Comparison 

[from  Ref.  26] . 48 

Table  12.  Russell’s  Error  Factor  Using  High  Pass  Filtered  Data  (420  msec)  [from 

Ref.  26] . 49 

Table  13.  Russell’s  Error  Factor  Using  Drift  Compensated  Data  (500  msec) . 50 

Table  14.  Summary  of  Prediction  Error  Measures  Using  Dense  2nd  Generation 

Model . 52 

Table  15.  Summary  of  Simulation  Error  Measures  using  Dense  3rd  Generation 

Model . 54 

Table  16.  Shot  3  Fluid  Model  Properties . 58 

Table  17.  Shot  3:  Russell’s  Error  Factor . 61 

Table  18.  Shot  1  Russell  Error  Factors  With  Varying  Fluid  Model  Depth . 68 

Table  19.  Shot  1  Fluid  Model  Mesh  Properties . 69 

Table  20.  Shot  1 :  Russell’s  Error  Factor . 73 

Table  2 1 .  Shot  2  Fluid  Model  Properties . 80 

Table  22.  Shot  2:  Russell’s  Error  Factor . 83 

Table  23.  Accelerometer  vs.  Velocity  Meter  Correlation  Comparison . 90 

Table  24.  2X  Cavitation  Fluid  Mesh  Russell’s  Error  Comparison  (120  msec) . 91 

Table  25.  2X  Cavitation  Depth  Fluid  Mesh  Russell’s  Error  Comparison  (120  msec)  ....92 

Table  26.  Russell’s  Error  Comparison:  Cavitation  Depth  Model  vs.  2X  Cavitation 

Depth  Model  (500  msec) . 93 

Table  27.  Lateral  (Athwartship)  Accelerometer  Nodal  Locations . 96 

Table  28.  Longitudinal  (Fore-Aft)  Accelerometer  Nodal  Locations . 96 

Table  29.  Russell’s  Error  Factor  for  Lateral  Velocity  Response  Comparisons . 97 

Table  30.  Russell’s  Error  Factor  for  Longitudinal  Velocity  Response  Comparisons....  101 


xiii 


ACKNOWLEDGEMENTS 


The  authors  are  grateful  for  full  sponsorship  given  by  Naval  Sea  Systems 
Command,  AEGIS  Program  Office(PMS  400D).  The  authors  also  extend  their  gratitude 
to  CAPT  David  L.  Lewis,  CDR  Jeffery  S.  Riedel,  Mr.  Constintine  Constant  and  Mr. 
Michael  Winnette  for  their  continuous  discussions  and  encouragement  throughout  the 
course  of  DDG  8 1  ship  shock  modeling  and  simulation.  We  would  also  like  to  thank  Dr. 
R.  Rainsberger  for  his  support  in  fluid  volume  modeling  using  TrueGrid.  Thanks  are  also 
due  to  Mr.  T.  Christian  for  his  excellent  job  to  make  the  computer  always  on-line  and 
lastly  to  NPS  Shock  team  for  their  tireless  work. 


I.  INTRODUCTION 


A.  BACKGROUND 

In  World  War  II,  the  U.S.  Navy  experienced  the  highly  destructive  effects  of  near 
proximity  underwater  explosions  (UNDEX)  from  mines  and  torpedoes.  Many 
combatants  with  the  latest  in  combat  technology  for  the  time  were  rendered  helpless  due 
to  inadequate  shock  proofing  of  the  ship  systems.  Since  this  time,  extensive  work  has 
gone  into  the  research  and  study  of  the  effects  of  UNDEX.  A  major  goal  in  the  design  of 
modem  combatant  ships  has  been  to  eliminate  or  at  least  reduce  damage  caused  by 
UNDEX. 

Guidelines  and  specifications  have  been  developed  for  the  shock  testing  and 
hardening  of  shipboard  equipment  and  systems.  NAVSEA  0908-LP-000-3010A  [Ref.  1] 
and  MIL-S-901D  [Ref.  2]  set  forth  much  of  this  guidance.  As  directed  by  OPNAVINST 
9072.2  [Ref.  3],  the  complete  ship  system  is  tested  by  perfonning  what  are  called 
underwater  shock  trials.  These  shock  trials  attempt  to  test  the  ship  under  “near  combat 
conditions”  by  igniting  a  xxxx  pound  charge  of  HBX-1  underwater  at  varying  standoff 
distances  from  the  ship.  The  effect  of  the  shocks  to  ship  systems  is  observed  and  the 
response  of  the  ship  is  monitored  and  recorded  for  each  shot.  The  lead  ship  of  each  class, 
or  a  ship  substantially  deviating  from  other  ships  of  the  same  class,  is  required  to  undergo 
these  trials  in  order  to  correct  any  deficiencies  on  that  ship  as  well  as  the  follow  on  ships 
of  the  class. 

These  shock  trials,  while  beneficial  in  detennining  the  wartime  survivability  of 
surface  ships,  require  years  of  planning  and  preparation  and  are  extremely  expensive.  In 
the  Aegis  Destroyer  program  alone,  over  30  million  dollars  were  spent  for  the  shock  trials 
conducted  on  USS  JOHN  PAUL  JONES  (DDG  53)  in  1994  and  another  20  million 
dollars  for  the  shock  trials  conducted  on  USS  WINSTON  S.  CHURCHILL  (DDG  81)  in 
2001.  In  addition,  these  tests  present  an  obvious  danger  to  the  crew  onboard,  the  ship 
itself,  and  any  marine  life  in  the  vicinity  of  the  test.  Due  to  this  inherent  safety  risk, 
shock  trials  do  not  test  up  to  the  ship’s  design  limits  or  even  the  true  combat  shock 
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environment.  This  has  raised  the  question  as  to  whether  or  not  the  information  gleaned 
from  doing  the  tests  is  worth  the  high  cost  of  conducting  them  [Ref.  4], 

Advances  in  the  capabilities  of  computers  in  the  last  few  decades  have  allowed 
many  events  to  be  tested  in  a  virtual  world.  This  has  allowed  engineers  not  only  to  save 
millions  of  dollars,  but  it  has  also  allowed  for  faster  advancements  in  design.  Creating  a 
virtual  UNDEX  environment  for  ship  system  underwater  shock  testing  presents  an 
extremely  useful  design  tool  and  an  attractive  alternative  to  future  shock  trials.  This 
would  save  the  Navy  millions  of  dollars  while  at  the  same  time  not  harm  our  sailors  or 
the  environment.  To  date  however,  virtual  UNDEX  testing  of  ship  systems  has  not  been 
considered  sufficiently  reliable  and  the  U.S.  Navy,  though  very  interested  in  the 
simulations,  has  not  considered  it  an  acceptable  alternative  to  actual  shock  trials. 

To  obtain  accurate  results  from  a  simulation,  the  structural  finite  element  model 
must  be  extremely  detailed.  Additionally,  the  surrounding  acoustic  medium  must  be 
coupled  with  the  wetted  surface  of  the  structural  model  almost  perfectly.  The  resulting 
coupled  fluid  and  structural  model  necessary  for  an  accurate  simulation  is  enormous 
computationally.  The  UNDEX  environment  is  very  complex,  made  up  of  an  initial 
“kick”  from  the  incident  shock  wave  followed  by  the  effects  of  cavitation,  bubble  pulse, 
and  structural  whipping  [Ref.  4],  While  the  computational  time  step  must  be  very  small, 
on  the  order  of  microseconds,  the  entire  UNDEX  ship  response  lasts  on  the  order  of 
seconds.  Because  of  this,  current  computer  simulations  cannot  realistically  predict  the 
complete  response  from  start  to  finish. 

Current  technology  is  capable  of  creating  a  simulation  to  predict  the  initial  “early 
time”  peak  response  of  a  surface  ship  subjected  to  an  underwater  shock  [Ref.  4].  By 
capturing  the  peak  response,  simulations  can  provide  very  useful  insight  into  ship  system 
design.  Because  they  are  virtual  in  nature,  shock  simulations  can  be  conducted  at  or 
beyond  design  limits,  providing  more  useful  design  information  than  that  which  is  gained 
by  conducting  real  world  shock  trials.  By  validating  the  predictions  made  using  UNDEX 
simulation,  virtual  testing  can  be  used  to  improve  and  accelerate  combatant  ship  system 
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design.  Through  this  validation  and  further  advancements  in  computer  processing 
technology,  virtual  shock  trials  may  one  day  reduce  or  eliminate  the  requirement  for  live 
fire  testing  of  entire  ship  systems. 


B.  SCOPE  OF  RESEARCH 

This  paper  investigates  the  feasibility  of  creating  a  virtual  UNDEX  environment 
to  model  the  response  of  a  surface  ship  exposed  to  an  underwater  shock.  This  is 
accomplished  using  the  shock  trial  data  from  the  shock  trial  conducted  on  USS 
WINSTON  S.  CHURCHILL  (DDG  81)  in  June  of  2001.  DDG  81  was  subjected  to  three 
different  underwater  explosions  during  these  trials,  each  explosion  having  a  different 
geometry.  Comparisons  between  the  simulated  response  predictions  and  the  actual 
measured  response  are  made  for  all  three  shots.  The  effects  of  charge  location  in  relation 
to  the  ship,  size  of  the  fluid  mesh,  fluid  mesh  density,  and  the  quality  of  the  fluid  mesh 
are  investigated.  By  creating  an  easily  modified  fluid  mesh  model,  a  wide  range  of  shot 
and  ship  loading  conditions  can  be  simulated.  This  creates  unlimited  potential  for  future 
UNDEX  simulation  research.  The  modifiability  of  the  fluid  models  investigated  in  this 
paper  and  the  modeling  techniques  that  pennit  this  easy  modification  are  discussed. 

These  underwater  shock  simulations  were  conducted  using  the  most  accurate  and 
state  of  the  art  finite  element  ship  model  to  date.  The  fluid  volume  models  were 
generated  using  the  TrueGrid  finite  element  mesh  generator  program  [Ref.  5].  The 
simulation  processing  was  conducted  using  the  LS-DYNA/USA  (Underwater  Shock 
Analysis)  coupled  computer  code  [Ref.  6,7]. 
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II.  UNDERWATER  SHOCK  MODELING  AND  SIMULATION 


A.  STRUCTURAL  MODEL 

1.  USS  WINSTON  S.  CHURCHILL  (DDG  81) 

A  finite  element  model  (FEM)  of  USS  WINSTON  S.  CHURCHILL  (DDG  81) 
built  by  Gibbs  &  Cox,  Inc  [Ref.  8]  was  used  for  the  research  summarized  in  this  report. 
This  model,  shown  in  Figure  1,  is  a  modification  of  the  USS  JOHN  PAUL  JONES  (DDG 
53)  model  used  in  the  UNDEX  simulations  conducted  at  the  Naval  Postgraduate  School 
(NPS)  from  1999  to  2000  [Ref.  9,  10.]. 


Figure  1.  USS  WINSTON  S.  CHURCHILL  (DDG  81)  FINITE  ELEMENT 

MODEL  [from  Ref.  8]. 
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DDG  81  is  a  Flight  IIA  Arleigh  Burke  Class  destroyer  having  significant 
modifications  from  prior  destroyers  in  the  class.  The  major  structural  modifications 
include: 

•  the  addition  of  dual  helicopter  hangars 

•  six  additional  vertical  launch  system  (VLS)  cells,  three  forward  and  three 
aft 

•  installation  of  a  5 ’762  caliber  gun 

•  a  five  foot  extension  of  the  transom 

•  addition  of  five  blast  hardened  bulkheads 

Delivered  in  MSC/NASTRAN  input  deck  form,  this  model  represents  the  largest 
and  most  complex  finite  element  ship  model  to  date  [Ref.  11].  Though  based  on  the 
DDG  53  model,  this  model  was  greatly  improved.  Structural  beam  elements  were  refined 
to  reflect  true  areas,  volumes,  and  densities  and  end  constraints  were  modified  to  better 
reflect  actual  structural  boundary  conditions.  The  mass  distribution  was  improved 
significantly.  Following  completion  of  the  model,  nodal  masses  were  added  throughout 
the  model  to  reflect  such  masses  such  as  weapons  loadout,  liquid  loading,  equipment,  and 
even  personnel  loading.  Rather  than  the  nonnal  smearing  of  this  mass  arbitrarily  onto  the 
entire  model,  masses  were  distributed  onto  nodes  in  the  model  that  corresponded  to  the 
associated  mass  on  the  actual  ship.  Since  the  three  underwater  shocks  were  accomplished 
on  separate  days,  three  separate  models  were  delivered,  each  with  a  unique  mass 
distribution  applicable  to  the  actual  ship  condition  when  the  shot  was  conducted.  Table  1 
is  a  summary  of  this  finite  element  model. 

Complex  equipment  such  as  the  gas  turbine  engines,  shafts,  generators,  reduction 
gear,  VLS  launchers,  and  the  gun  mount  were  modeled  as  rigid  body  elements.  Since  the 
internal  response  of  this  equipment  is  not  measured,  modeling  these  individual 
components  in  detail  is  unnecessary.  Each  of  these  element  groups  was  attached  to  the 
nodal  points  on  the  ship  structure  at  their  foundations  with  their  spring-mass 
representation  applied  to  the  supporting  nodes.  To  separate  the  propulsion  shafts  from 
the  reduction  gear,  a  beam  element  was  added  as  a  link  between  the  two  groupings.  The 
shaft  mass  was  then  represented  by  each  of  the  shaft  bearing  assemblies.  Figure  8  shows 
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the  significant  rigid  body  element  groups.  A  total  of  107  total  rigid  body  element  groups 
were  in  the  final  models  used  in  the  simulations.  In  addition  to  those  shown,  many  of 
these  groups  were  small  groups  of  one,  two,  or  three  elements  that  make  up  antennas  on 
the  mast  and  the  ship's  Recovery,  Assist,  Securing,  and  Traversing  (RAST)  system 
located  on  the  flight  deck. 


Table  1.  DDG  81  Finite  Element  Model  Specifications  [from  Ref.  8] 


Nodes 

40,514 

Degrees  of  Freedom 

243,084 

Beam  Elements 

49,397 

Shell  Elements 

48,662 

Spring  Elements 

416 

Rigid  Body  Element  Groups 

107 

Lumped  Masses 

92,541 

Propulsion  Shafting 

Systems 


Fwd  M  S 
Module 


5”  Gun 


Aft  VLS  Module  PGT  #2  A  &  #2B  MRG  #2 


Radar 

Array 


SSGTG  #3 


Ships  Service  Gas 
Turbine  Generator 
(SSGTG)  #1 
Propulsion  Gas  Turbine 
(PGT)  #1A  &  #1B  & 

Main  Reduction  Gear 
(MRG)  #1 


Figure  2.  Equipment  Modeled  as  Rigid  Bodies  in  DDG  81  FEM  [from  Ref.  8] 
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2.  Ship  System  Damping 

Ninety  percent  of  damping  in  a  structure  occurs  due  to  frictional  energy 
dissipation  in  bolted  or  riveted  mechanical  joints.  The  vast  majority  of  joints  on  naval 
combatant  ship  structure  systems  are  welded  while  all  stiffeners  are  welded  directly  to  the 
hull  plates,  decks,  and  bulkheads.  These  welded  joints  provide  minimal  energy 
dissipation  [Ref.  12].  Regardless,  many  energy  dissipation  sources  such  as  long  cable 
trays,  hangers,  snubbers,  and  the  surrounding  fluid  coupled  with  the  hull  are  inherent  to 
all  ship  systems.  While  ship  system  damping  is  measurable,  it  is  difficult  to  quantify  for 
modeling  purposes.  Extensive  study  was  conducted  at  NPS  in  conjunction  with  the  ship 
shock  simulation  of  USS  JOHN  PAUL  JONES  (DDG  53)  using  2000  msec  of  measured 
shock  trial  data  [Ref.  13].  The  results  of  this  study  were  applied  to  the  Rayleigh  damping 
coefficients  used  for  the  DDG  8 1  model 

Rayleigh  damping  is  a  general  fonn  of  proportional  damping  representing  the 
damping  matrix,  [C],  in  the  structural  equation  of  motion. 

[M]{x}+[C]{x}+[K]x={F}  (i) 

The  [C]  matrix  is  represented  by 

[C]  =<x[M]+p[K]  (2) 

where  a  and  P  are  constants.  This  equation  can  be  normalized  using  the  mass 
nonnalization  matrix  [</;] 

MT[C]M  =  [2^T]diag  =  a[I]  +  p[cj]diag  (3) 

Determining  these  damping  coefficients  for  a  simple  system  having  only  two 
modes  with  two  modal  frequencies  of  interest  is  academic.  Determining  them  for 
complex  systems  such  as  a  ship  structure  with  more  than  two  modes  of  interest  presents  a 
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much  greater  challenge.  In  this  case,  Equation  (3)  becomes  over  determined,  having 
more  equations  than  unknowns.  These  coefficients  can  be  determined  using  measured 
data  and  a  least  squares  curve  fitting  technique. 


For  comparison  with  measured  ship  response  the  modal  damping  ratio  is 
calculated.  This  is  determined  for  each  mode  by 


1  CO; 

Ci  =  a - +  p  —  (4) 

2c0j  2 


DDG  53  was  divided  into  67  area  groups  for  the  comparison  study.  Measured 
modal  response  was  recorded  over  the  entire  frequency  spectrum  of  interest,  0  to  250  Hz, 
for  both  athwartship  and  vertical  response.  A  least  squares  curve  fit  was  then  applied  and 
the  corresponding  a  and  P  were  recorded  for  each  area.  Figure  3  shows  an  example  of 
this  curve  fitting  for  one  area  group.  This  particular  plot  includes  the  measured  data 
outliers  and  “noise”  that  were  subsequently  discarded  in  performing  the  final  curve  fit. 


Figure  3.  Vertical  Modal  Damping  Ratio  for  Single  Area  Group  of  DDG  53  [from 

Ref.  13] 
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Each  area  group  was  given  a  weight  based  upon  the  number  of  modes  used  in 
determining  the  least  squares  curve  fit  for  that  area  group.  These  values  were  summed 
and  divided  by  the  sum  of  the  modes  used  to  find  both  the  weighted  amean  and  the 
weighted  pmean.  These  values  are  shown  in  Tables  2  and  3.  The  recommended  values  for 
a  and  p  for  DDG  53  were  determined  to  be  19.2  and  2.09e-6  respectively.  These  values 
were  chosen  due  to  the  fact  that  the  vertical  response  is  much  larger  in  magnitude  than 
the  athwartship  response.  The  resulting  damping  ratio  curve  for  both  response  directions 
is  shown  in  Figure  4  [Ref.  13]. 


Table  2.  Weighted  Mean  of  a  (DDG  53)  [from  Ref.  131 


Athwartship  Direction 

Vertical  Direction 

18.4 

19.2 

Table  3.  Weighted  Mean  of 


B  (DDG  53)  [from  Ref.  131 


Athwartship  Direction 

Vertical  Direction 

2.82E-6 

2.09E-6 

From  a  structural  damping  perspective,  DDG  53  and  DDG  81  are  virtually 
identical.  Due  to  this  similarity  these  values  were  used  for  the  DDG  81  model  [Ref.  13]. 
Application  of  Rayleigh  damping  to  the  model  is  performed  in  the 
“GFOBAF  DAMPING”  card  of  FS-DYNA  [Ref.  6].  All  structural  shell  and  beam 
elements  in  the  model  were  assigned  the  same  proportional  damping  value. 
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Figure  4.  Damping  Ratio  Curve  Plot  for  Using  Recommended  Rayleigh  Damping 

Coefficients  for  DDG  53  [from  Ref.  13] 

B.  SIMULATION  CODE 

Several  analysis  programs  were  used  in  the  modeling  and  simulation  of  DDG  8 1 . 
The  model  building,  simulation,  and  data  analysis  procedure  is  outlined  on  the  following 
page  in  Figure  5. 

1.  Model  Generation  and  Pre-Processing 

The  USS  WINSTON  S.  CFIURCHILL  structural  finite  element  model  was 
originally  constructed  by  Gibbs  &  Cox,  Inc  using  FEMAP.  FEMAP  is  a  Windows-based 
pre-  and  post-processor  used  for  engineering  finite  element  analysis.  It  is  CAD  neutral 
and  solver  neutral,  making  it  a  very  useful  modeling  tool  for  a  wide  range  of  applications 
[Ref.  14].  The  model  was  output  from  FEMAP  in  MSC/NASTRAN  format  and 
delivered  in  this  form.  MSC/NASTRAN/PATRAN  was  then  initially  used  to  identify 
elements  and  element  properties  and  as  a  quality  control  check  to  verify  all  element  and 
material  properties  were  properly  defined. 
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Figure  5.  Model  Generation  and  Simulation  Flow  Chart 


The  actual  simulations  were  performed  using  LS-DYNA  coupled  with  USA.  This 

required  that  the  model  be  translated  from  NASTRAN  format  to  LS-DYNA  keyword 

format.  This  was  a  complicated,  but  critical  step  that  was  done  manually  and  through  a 
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conversion  code  written  in  FORTRAN.  For  construction  of  the  surrounding  fluid  model, 
the  structural  model  was  imported  directly  into  the  TrueGrid  mesh  generator  in 
NASTRAN  format.  More  will  be  discussed  on  the  construction  of  the  fluid  mesh  using 
TrueGrid  in  Chapter  III. 

2.  Simulation  Processing 

LS-DYNA  (Version  960)  and  USA  were  used  to  simulate  the  underwater  shock 
and  solve  for  the  response  of  the  structural  model.  LS-DYNA  is  a  computational  tool 
very  popular  in  the  automotive  industry  where  it  is  commonly  used  to  simulate  such 
events  as  automobile  crashes  and  airbag  deployment.  It  is  an  explicit  and  implicit  finite 
element  program  capable  of  analyzing  the  nonlinear  dynamic  response  of  three- 
dimensional  inelastic  structures  [Ref.  6], 

Working  in  concert  with  LS-DYNA,  USA  was  used  to  calculate  the  transient 
response  of  the  ship’s  wetted  surface  structure  to  the  incident  shock  wave.  USA  uses  the 
Doubly  Asymptotic  Approximation  (DAA)  to  solve  the  fluid  structure  interaction 
equations.  This  approximation  method  begins  first  with  the  structural  equation  of  motion 

[M]{x|  +  [C]{x}  +  [K]{x|  =  {f}  (5) 

where  [M],  [C],  and  [K]  make  up  the  symmetric  linear  structural  mass,  damping 
and  stiffness  matrices  respectively  and  {f}  is  the  external  forcing  vector.  For  a 
submerged  structure  excited  by  an  acoustic  wave  this  external  forcing  vector  is 

{f}  = -[G][Af  ]{(pi+ps)}  +  {fD}  (6) 

where  {pi}  and  {ps}  are  the  nodal  incident  and  scattered  pressure  vectors 
respectively  for  the  wetted  surface  structure,  {fd}  is  the  dry  structure  applied  force  vector, 
[Af]  is  the  diagonal  area  matrix  pertaining  to  the  fluid  elements,  and  [G]  is  the 
transfonnation  matrix  that  relates  the  structural  and  fluid  nodal  surface  forces.  While  the 
incident  pressure  vector  is  known,  the  scattered  pressure  vector  is  unknown  [Ref.  15]. 
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The  DAA  is  used  to  determine  the  fluid  structure  interaction.  The  DAA  is  so 
named  as  it  approaches  exactness  in  both  the  high-frequency  (early  time)  and  low 
frequency  (late  time)  limits  [Ref.  16].  This  approximation  uses  two  equations,  the 
structural  equation  of  motion  and  a  fluid  particle  equation, 

[Mf  ]{ps}+pc[Af  ]{ps}  =pc[Mf]{us}  (7) 

where  {us}  is  the  scattered-wave  fluid  particle  velocity  vector  nonnal  to  the 
structural  wetted  surface,  [Mf]  is  the  symmetric  fluid  mass  matrix  for  the  fluid  mesh  on 
the  wet  surface,  and  c  is  the  acoustic  velocity.  The  DAA  is  able  to  represent  the  entire 
fluid  surrounding  a  submerged  structure  through  variables  associated  only  with  the 
wetted  surface. 

The  fluid  particle  velocities  are  related  to  the  structural  response  through  the 
compatibility  relation 


[G]T  {x}  =  ]u,}+{us}  (8) 

This  equation  contains  the  constraint  that  the  nonnal  fluid  particle  velocity  must  match 
the  normal  structural  velocity  on  the  wet  surface  of  the  structure  [Ref.  16]. 

Equations  (6)  through  (8)  can  now  be  combined  to  form  two  equations 

[Ms]{x}+[Cs]{x}+[Ks]{x]  =-[G][Af]{(Pl+Ps)}  (9) 

[Mf  ]{Ps}+Pc[Af  ]{Ps}=  Pc[Mf  ]([G]T  {x}-{u,})  (10) 

Equations  (9)  and  (10)  are  solved  using  a  staggered  solution  scheme  in  the 
Underwater  Shock  Analysis  (USA)  code  in  order  to  approximate  the  forces  at  the  fluid 
structure  interface  [Ref.  7].  These  forces  can  then  be  input  into  a  finite  element  solver  to 
determine  the  structural  response  throughout  the  ship  structure. 
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This  method  allows  the  problem  to  be  solved  without  requiring  a  large  fluid 
volume  mesh.  This  works  very  well  for  incident  shock  wave  predictions,  but  does  not 
consider  cavitation,  making  it  ideal  for  structures  not  affected  by  cavitation  such  as  a 
deeply  submerged  submarine. 

Cavitation  has  a  major  impact  on  the  response  of  a  surface  ship  subjected  to  an 
underwater  shock,  particularly  in  the  late  time  response.  A  fluid  volume  finite  element 
model  (FEM)  must  surround  the  ship  an  adequate  distance  in  order  to  capture  the  effects 
of  bulk  and  local  cavitation.  Previous  studies  have  shown  that  the  minimum  depth  of  this 
fluid  mesh  should  be  at  least  one  half  the  cavitation  depth  for  a  reasonable  approximation 
[Ref.  9,  10].  This  fluid  volume  model  should  be  extruded  from  the  wetted  surface 
structure,  matching  the  structural  element  faces  and  node  points  as  perfectly  as  possible 
[Ref.  17].  The  DAA  boundary  is  then  placed  on  the  outer  surface  of  the  fluid  volume. 
This  requires  considerably  more  computational  resources,  but  the  improvement  in  the 
predicted  response  is  substantial,  making  it  a  necessary  requirement  [Ref.  18].  The  USA 
code  consists  of  three  components,  FLUMAS;  AUGMAT;  and  T1MINT. 
ci.  FLUMAS 

This  component  of  USA  constructs  the  fluid  mass  matrix  for  the  wetted 
surface  structure  using  the  boundary  element  technique,  ft  sets  the  fluid  mesh  properties, 
such  as  density,  acoustic  velocity,  atmospheric  pressure,  and  gravity.  The  set  of 
transfonnation  coefficients  that  relate  to  the  structural  and  fluid  degrees  of  freedom  on 
the  wetted  surface  structure  are  determined  in  this  component.  It  also  has  the  capability 
to  solve  the  fluid  boundary  eigenvalue  problem  and  automatically  calculates  the  rigid 
body  added  mass  coefficients  [Ref.  7,  19]. 
b.  AUGMAT 

The  purpose  of  this  component  is  to  use  the  data  detennined  by  the 
FLUMAS  component  to  determine  the  compatibility  relation  matrix  used  in  the  staggered 
solution  procedure  described  in  Chapter  II.  This  solution  procedure  is  used  to  detennine 
the  transient  response  of  a  submerged  structure  [Ref.  19]. 
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c.  TIMINT 

The  TIMINT  component  of  USA  is  where  the  majority  of  the  simulation 
takes  place.  It  conducts  a  step-by-step  direct  numerical  time  integration  of  the  wetted 
surface  structure  exposed  to  the  spherical  shock  wave.  This  component  utilizes  the  DAA 
described  in  Chapter  II  where  a  staggered  solution  technique  is  used  solving  the  structural 
and  fluid  response  equations  separately  [Ref.  19].  Time  histories  of  the  structural  and 
fluid  response  are  output  from  this  component.  Also  output  from  this  component  are  the 
plot  files  necessary  to  perform  the  simulation  animation.  The  vast  majority  of  the 
computational  resources  used  for  the  simulations  were  used  by  this  component. 

3.  Simulation  Post  Processing  and  Data  Extraction 

After  LS-DYNA  and  USA  have  completed  processing  the  model,  the  data  stored 
in  the  output  files  were  extracted  using  GLview.  GLview  is  a  very  capable  3D 
visualization  and  animation  post-processor  [Ref.  20].  The  simulations  required 
tremendous  quantities  of  data  to  be  processed  due  to  the  extremely  small  calculation  time 
step  necessary.  GLview  is  capable  of  handling  these  large  quantities,  producing  a 
spectacular  interactive  animation  of  the  virtual  shock  environment.  In  addition,  it  was 
used  to  download  and  display  the  structural  response  time  history  graphically.  It  served 
as  an  interface  between  the  raw  simulation  data  output  from  USA  and  the  data 
comparison  tools  used  to  compare  the  simulation  predictions  with  measured  test  results. 
These  raw  files  were  loaded  into  GLview  and  then  structural  response  parameters  such  as 
nodal  velocities  and  accelerations  were  exported  as  ASCII  files  for  further  analysis. 

4.  Data  Processing  and  Comparison 

The  nodal  velocities  and  accelerations  exported  from  GLview  required  filtering  to 
remove  undesirable  frequencies.  This  was  accomplished  by  importing  the  data  files  into 
a  MATLAB  code  which  processed  the  data  through  a  low  pass  Butterworth  filter.  This 
process  and  the  reasons  for  doing  it  will  be  discussed  in  greater  detail  in  Chapter  V.  The 
data  was  then  imported  into  UERDTools  for  comparison  with  measured  shock  trial  data. 
UERDTools  was  developed  by  the  Underwater  Explosions  Research  and  Development 
(UERD)  project  office  in  Bethesda,  Maryland  for  the  purpose  of  analyzing  shock  trial 
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data.  It  has  numerous  features  specific  to  ship  shock  analysis  such  as  unit  scaling  of  data, 
filtering,  interpolation,  time  shift,  and  error  criteria  calculations.  More  will  be  discussed 
about  many  of  these  features  in  Chapter  IV. 
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III.  FLUID  MESH  MODELING 


Fluid  mesh  modeling  is  the  primary  focus  of  this  report.  Several  fluid  models 
were  constructed  to  model  the  three  different  shock  environments.  These  models  were  all 
constructed  using  TrueGrid.  TrueGrid  is  a  powerful  mesh  generator  that  allows  meshes 
of  all  types  to  be  created  interactively.  While  its  full  capabilities  cannot  be  overstated, 
one  key  feature  that  makes  it  ideal  for  modeling  the  underwater  shock  environment  is  its 
ability  to  extrude  the  fluid  mesh  through  the  structural  mesh.  This  results  in  a  fluid  mesh 
that  aligns  perfectly  with  the  structural  mesh  from  the  fluid- structure  interface  out  to  the 
Doubly  Asymptotic  Approximation  (DAA)  boundary.  While  several  models  were  used 
in  the  simulations,  all  of  these  models  were  variations  of  just  three  models.  The 
construction  of  the  third  model,  an  improvement  over  the  first  two  models,  is  highlighted 
in  this  chapter  and  described  in  greater  detail  in  Appendix  B.  A  list  of  many  of  the  useful 
commands  used  to  generate  this  model  in  TrueGrid  is  contained  in  Appendix  C. 


A.  FLUID  MESH  CONSTRUCTION 

The  fluid  mesh  was  constructed  using  essentially  four  subdivisions  merged 
together;  the  inner  liner,  the  inner  mesh,  the  transitional  mesh,  and  the  outer  mesh.  By 
building  the  fluid  mesh  in  this  manner,  several  goals  were  accomplished: 

•  USA  code  stability  criteria  was  met 

•  The  fluid  volume  was  meshed  perfectly  orthogonal  to  the  structural  mesh 

•  A  quality  mesh  was  maintained  throughout  the  fluid  volume  from  the  fluid 
structure  interface  out  to  the  boundary 

•  The  fluid  model  was  easily  modified 

•  The  proper  waterline  was  achieved 

1.  Inner  Liner 

The  inner  liner  of  the  fluid  mesh  serves  as  the  fluid  structure  interface.  It  is  the 
most  critical  part  of  the  fluid  mesh.  The  liner  must  mesh  perfectly  with  the  structural 
mesh,  it  must  be  as  orthogonal  to  the  hull  as  possible,  and  it  must  meet  USA  code 
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stability  criteria.  Regardless  of  the  quality  of  the  remainder  of  the  fluid  mesh,  a  poorly 
constructed  inner  liner  will  result  in  a  poor  approximation.  This  segment  of  the  fluid 
mesh  consists  of  37  subparts  and  accounts  for  approximately  50%  of  the  fluid  mesh 
model  code  written. 

a.  USA  Code  Stability  Criteria 

The  nodal  spacing  of  this  inner  liner  is  critical  to  the  stability  of  the  USA 
code.  The  code  requires  that  the  nodal  spacing  of  this  inner  layer  cannot  exceed  a 
defined  thickness,  D,  given  by  the  following  equation: 


<5  (li) 

Psls 

where  p  is  the  density  of  water  and  ps  and  ts  are  the  density  and  thickness 
respectively  of  the  wetted  surface  shell  elements. 

The  DDG  81  model  wetted  surface  was  constructed  using  a  variety  of 
quadrilateral  and  triangle  shell  elements.  A  listing  of  these  elements  and  their  associated 
thickness  is  shown  in  Table  4. 


Table  4.  DDG  81  Wetted  Surface  Shell  Element  Properties 


Shell  Part  ID 

Material  ID 

Thickness  (in) 

1062 

2 

0.375 

1072 

2 

0.4375 

1073 

3 

0.4380 

1082 

2 

0.5 

1102 

2 

0.625 

1112 

2 

0.6875 

1122 

2 

0.75 

1123 

3 

0.75 

1124 

3 

0.75 

1163 

3 

1.0 

1202 

2 

1.25 
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As  shown  in  Table  4,  the  thinnest  shell  element  used  is  Shell  Part  ID  1062 
which  is  0.375  inches  thick.  The  material  density  of  this  shell  element  is  ps=7.324E-4 
(lbf-s2/in4)  while  the  density  used  for  the  water  is  p=9.3601E-5  (lbf-s2/in4).  Applying  the 
values  in  Equation  (11)  results  in  the  maximum  allowable  thickness  of  the  inner  liner 
being  7.33  inches.  To  ensure  that  this  was  not  exceeded,  a  thickness  of  seven  inches  was 
used  for  the  entire  inner  liner. 

b.  Orthogonality 

The  majority  of  the  inner  liner  was  made  orthogonal  to  the  hull  by  using 
TrueGrid’s  “BLUDE”  command  along  with  the  “NORMAL”  command.  The  BLUDE 
command  extrudes  the  fluid  mesh  through  the  structural  mesh  resulting  in  a  fluid  mesh 
which  is  merged  perfectly  with  the  structural  mesh.  While  the  BLUDE  command  does 
generally  achieve  nonnality,  the  NORMAL  command  was  used  to  ensure  that  the  fluid 
mesh  was  nonnal  to  the  hull  over  its  entirety  and  only  seven  inches  in  thickness. 

The  abrupt  structural  transitions  such  as  the  keel,  bow,  sonar  dome,  and 
stern  areas  did  not  allow  for  a  single  continuous  fluid  liner  to  be  extruded  from  the  hull. 
In  these  areas  the  liner  was  first  extruded  on  all  faces  around  the  transition  area.  The 
gaps  that  were  left  between  these  meshes  were  filled  in  using  wedge  shaped  elements. 
This  ensured  that  orthogonality  was  maintained  at  these  transitions  while  filling  in  the 
gap  to  allow  for  further  construction  of  the  fluid  mesh. 

c.  Waterline 

The  proper  waterline  is  essential  for  two  reasons,  to  ensure  proper  global 
orientation  of  the  fluid  model  and  to  accurately  model  the  fluid  structural  interface.  The 
DDG  81  structural  model  was  constructed  without  mesh  lines  at  the  test  waterline.  To 
ensure  that  the  fluid  mesh  was  merged  properly  with  and  orthogonal  to  the  structural 
mesh  at  the  waterline,  wedge  shaped  elements  were  constructed  to  transition  from  the 
structural  mesh  line  to  the  actual  waterline.  Fluid- structure  orthogonality  was 
compromised  in  the  aft  section  of  the  ship  due  to  the  extreme  transition  from  the 
structural  mesh  line  to  the  test  waterline.  To  minimize  this  compromise,  an  additional 
layer  was  constructed  to  encapsulate  the  inner  liner  and  continue  the  transition  to  the  test 
waterline  in  this  area.  This  second  layer  was  eight  inches  in  thickness,  bringing  the  total 
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thickness  of  the  inner  liner  to  fifteen  inches  nominally.  Figure  6  shows  the  completed 
fluid  mesh  two  layer  inner  liner. 


Figure  6.  Fluid  Model  Two  Layer  Inner  Liner 
2.  Inner  Mesh 

Many  of  the  contours  of  the  ship,  particularly  the  sonar  dome  region  and  areas 
around  the  keel,  form  concave  surfaces  and  have  sharp  transition  areas.  The  purpose  of 
this  inner  mesh  is  to  provide  a  transition  from  these  complex  contours  of  the  structural 
model  to  a  surface  convex  in  nature.  The  quality  of  the  entire  fluid  mesh  depends  greatly 
upon  the  quality  of  this  region.  Complex  control  of  the  fluid  mesh  through  the  use  of 
proper  vertices  placement,  curve  attachments,  and  surface  projections  is  key  to 
constructing  this  region  of  the  mesh.  Proper  use  of  mesh  interpolation  methods  is 
essential  in  this  region  to  obtain  a  quality  mesh.  TrueGrid  uses  several  interpolation 
techniques  to  ensure  a  quality,  uniform  mesh. 

a.  Thomas-Middlecoff  Relaxation 

Using  linear  iteration  as  described  in  Reference  [21],  Thomas-Middlecoff 
relaxation  (TME)  improves  the  mesh  by  solving  a  Poisson  equation.  This  relaxation 
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method  excels  in  propagating  the  edge’s  nodal  distribution  into  the  interior  of  the  mesh. 
In  doing  this,  it  attempts  to  maintain  orthogonality  throughout  the  mesh  [Ref.  5]. 
b.  Equipotential  Relaxation 

This  method  of  mesh  interpolation  treats  the  mesh  lines  as  solutions  to  the 
Laplace  equation  with  the  boundary  conditions  formed  by  the  boundary  nodes.  This 
allows  the  interior  volume  of  a  mesh  to  be  relaxed  while  maintaining  the  boundary  nodal 
distributions.  The  primary  purpose  of  using  this  interpolation  method  is  to  maintain  a 
defined  boundary  while  producing  a  uniform  interior  mesh  [Ref.  5], 

The  inner  mesh  was  constructed  to  allow  for  variation  of  mesh  density,  allowing 
the  number  of  elements  between  the  inner  and  outer  boundaries  to  be  changed  using 
parameters.  In  its  entirety,  five  separate  subparts  were  constructed  separately  and  merged 
together  to  form  this  region  of  the  fluid  model.  The  finished  inner  mesh  extended 
outward  a  nominal  distance  of  five  feet  from  the  inner  liner.  It  is  shown  in  Figure  7. 
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Figure  7.  Fluid  Model  Inner  Mesh 

3.  Transitional  Mesh 

The  transitional  mesh  encapsulates  the  inner  mesh  and  extends  out  to 
approximately  ten  feet  from  the  hull.  Conical  shaped  mesh  regions  were  used  to  extend 
the  bow  and  stem  areas  outward  to  merge  perfectly  with  the  contoured  outer  boundary  of 
this  mesh. 

The  fluid  model  has  the  ability  to  model  the  underwater  shock  environment  for 
DDG  81  under  varying  loading  conditions  by  allowing  the  waterline  to  be  adjusted.  The 
ability  of  the  model  to  adjust  to  these  different  conditions  occurs  in  this  region  of  the 
mesh.  This  is  accomplished  by  allowing  the  inner  mesh  to  “float”  along  with  the  ship 
model  inside  this  transitional  mesh.  The  transitional  mesh  then  transitions  from  the  fixed 
258  inch  waterline  of  the  inner  mesh  to  the  actual  test  waterline.  The  simulations 
investigated  in  this  thesis  all  had  a  test  waterline  of  258  inches  fore  and  aft  which  the 
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inner  liner  and  mesh  was  initially  designed  for.  This  feature  was  not  necessary  in  these 
simulations,  but  was  added  to  allow  for  future  underwater  shock  simulation  modeling  of 
DDG81. 

In  addition  to  allowing  for  variation  in  waterline,  this  region  also  allows  for 
variation  in  mesh  density.  It  consists  of  three  parts  constructed  separately  and  merged 
together  and  with  the  previously  constructed  fluid  model.  The  finished  transitional  mesh 
is  shown  in  Figure  8. 


Figure  8.  Fluid  Model  Transitional  Mesh 
4.  Outer  Mesh 

The  outer  mesh  extends  the  fluid  volume  out  to  the  DAA  boundary.  The  DAA 

boundary  is  formed  on  the  outer  surface  of  this  mesh  by  selecting  the  outer  faces  of  this 

mesh.  Constructed  to  allow  for  the  simulation  of  virtually  any  underwater  shock 

environment,  the  depth,  fluid  mesh  density,  waterline,  and  contour  of  this  outer  mesh  can 
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be  modified  with  just  a  few  simple  parameter  inputs.  Though  comprising  the  majority  of 
the  fluid  model,  this  region  was  constructed  using  only  one  part.  The  completed  outer 
mesh  is  shown  in  Figure  9.  The  complete  coupled  fluid  structure  model  is  shown  in 
Figure  10. 


Figure  9.  Fluid  Model  Outer  Mesh 
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Figure  10.  Complete  Coupled  Fluid-Structure  Model  (Shot  1) 
B.  FLUID  MODEL  PROCESSING 


The  completed  fluid  models  were  output  from  TrueGrid  in  LS-DYNA  “Keyword” 
format.  The  fluid  model  output  file  consists  of  the  nodal  coordinate  definitions,  the  fluid 
volume  element  nodal  assignments,  and  the  face  set  constituting  the  outer  DAA 
boundary.  Of  particular  importance  here  is  that  the  fluid  nodes  corresponding  to  the 
same  nodes  on  the  ship  cannot  be  redefined.  Additionally,  the  DAA  boundary  face  set 
must  be  defined  with  the  normal  pointing  outward  using  a  counter-clockwise  nodal 
assignment  [Ref.  7]. 

The  fluid  models  were  all  assigned  as  material  property  Material  Type  90  in  LS- 
DYNA,  an  acoustic  material.  The  LS-DYNA  acoustic  pressure  element  has  a  cavitation 
flag  which  can  be  turned  on  and  off.  When  on,  this  flag  prevents  negative  pressures  from 
occurring  by  cutting  off  any  negative  pressure  at  zero.  This  indicates  the  presence  of 
cavitation  [Ref.  6].  This  flag  was  turned  on  for  all  simulations  presented  in  this  thesis. 
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The  simulations  began  with  an  initial  pressure  field  established  and  defined  using 
a  card  in  LS-DYNA  specifically  for  explosive  detonations.  This  pressure  field  is 
determined  with  the  incident  pressure  wave  already  propagated  out  to  one  element  away 
from  the  closest  structural  node  to  the  charge.  This  fluid  node,  tenned  the  reference 
node,  was  identified  using  a  Microsoft  Excel  spreadsheet  and  defined  in  the  LS-DYNA 
keyword  file.  In  addition,  the  closest  point  on  the  DAA  boundary,  termed  the  standoff 
point,  was  identified  using  Excel  and  the  coordinates  defined  in  LS-DYNA.  The 
coordinates  of  the  charge,  standoff  point,  and  the  reference  node  were  defined  in  the 
TIMINT  component  of  USA  along  with  the  slant  distance,  the  radial  distance  from  the 
charge  to  the  standoff  point  on  the  boundary.  Examples  of  the  LS-DYNA  keyword  file 
and  USA  component  files  are  included  in  Appendix  D. 


28 


IV.  DATA  COLLECTION  AND  ANALYSIS 


Many  variables  exist  in  the  underwater  shock  environment,  one  of  those  being  the 
method  in  which  the  measured  data  is  collected  and  analyzed.  Prior  to  making  any 
comparisons  with  simulation  predictions,  the  data  must  be  properly  analyzed.  The 
existence  of  low  frequency  “drift”  and  high  frequency  “noise”  in  shocktrial  data  present  a 
formidable  challenge  to  overcome. 

A.  VELOCITY  RESPONSE  “DRIFT” 

The  velocity  meter  and  accelerometer  were  originally  intended  for  measuring 
short  duration  transient  motion  signals  over  a  limited  frequency  bandwidth.  The 
accelerometer  is  ideal  for  frequency  response  measurement  at  frequencies  far  below  its 
own  resonant  frequency  while  the  velocity  meter  is  more  suitable  for  measuring  the 
response  to  frequencies  far  above  its  own  resonant  frequency  [Ref.  22].  For  purposes  of 
underwater  shock  simulation  measurement  it  is  desirable  to  extend  the  measurement  time 
out  to  capture  more  of  the  response.  This  results  in  requiring  the  sensors  to  measure  a 
considerably  larger  range  of  frequencies  than  they  were  intended  for.  Velocity  and 
displacement  time  histories  acquired  by  integrating  acceleration  time  histories  have  a 
tendency  to  drift  due  to  errors  in  the  zero  reference  portion  of  the  time  history.  This 
problem  exists  as  well  in  velocity  meter  time  history  as  velocity  meters  require  seismic 
correction,  which  involves  integration  as  well.  As  a  result,  the  measured  data  displays  a 
substantial  amount  of  drift  as  shown  in  Figure  1 1 . 

A  significant  amount  of  study  was  conducted  to  determine  how  to  eliminate  this 
drift.  UERDTools  has  a  built-in  function  called  “Drift  Compensation”  which  uses  an 
algorithm  to  eliminate  this  drift.  Figure  12  shows  the  same  velocity  data  after  this  drift 
compensation  has  been  used.  The  advantages  of  using  this  technique  are  that  the 
magnitude  of  the  response  is  maintained  and  there  is  no  shift  in  phase.  The  disadvantages 
are  that  it  does  not  always  eliminate  the  drift  as  shown  in  Figure  12.  Further  curve  fitting 
where  the  analyzer  uses  his/her  own  judgment  to  detennine  where  the  data  begins  to  drift 
is  often  required.  This  requires  the  user  to  detennine  a  best-fit  polynomial  to  the  spurious 
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trend.  This  polynomial  curve  is  subsequently  subtracted  from  the  data.  This  technique 
has  proven  to  be  effective  in  “cleaning”  the  data,  but  it  invites  a  fair  amount  of  subjective 
interpretation  and  potential  for  biased  modification  of  the  data.  Ideally,  a  method 
involving  less  “hands-on”  manipulation  of  the  data  would  be  preferable. 
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Figure  11.  Raw  Velocity  Data  Showing  Drift 
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Figure  12.  Velocity  Data  After  Drift  Compensation  Has  Been  Applied 
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An  alternate  method  of  dealing  with  this  drift  is  to  process  the  data  through  a 
high-pass  filter.  While  velocity  meters  record  the  global  response,  it  is  most  desirable  to 
determine  only  the  local  response  at  a  particular  location.  To  determine  the  local 
response  it  is  therefore  desirable  to  eliminate  any  rigid  body  modes  from  the  sensor 
response.  This  can  be  measured  by  analyzing  the  late  time  displacement  response  of  the 
structure  where  rigid  body  modes  reveal  themselves  in  the  form  of  a  non-zero  steady  state 
displacement.  Doing  this  rigid  body  motion  study  for  a  number  of  sensors,  it  was 
initially  determined  that  the  measured  data  should  be  filtered  using  a  Bessel  two-pole 
high-pass  filter  set  at  a  cut-off  frequency  of  3Hz.  The  same  was  done  with  the  simulation 
data  where  a  cut-off  frequency  of  1.8Hz  was  necessary.  A  two-pole  Bessel  filter  was 
chosen  as  it  is  the  only  built-in  high-pass  filter  in  UERDTools.  Figure  13  shows  a 
comparison  of  filtered  versus  drift  compensated  velocity  data.  As  shown  in  Figure  14, 
filtering  results  in  a  velocity  plot  with  no  significant  drift  apparent,  while  the  drift 
compensated  data  still  shows  strong  evidence  of  low  frequency  drift. 


Vertical  Velocity 


Time  (msec) 

Drift  Compensated  - 3  Hz  High  Pass  Bessel  Filtered 


Figure  13.  Comparison  of  High  Pass  Filtering  vs.  Drift  Compensation 
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The  disadvantage  of  using  this  high-pass  filter  is  that  desirable  frequencies  above  3  Hz  are 
filtered  out  due  to  the  characteristic  of  the  two-pole  Bessel  filter.  Figure  14  shows  a 
frequency  response  of  a  two-pole  high-pass  Bessel  filter  set  at  a  cut-off  frequency  of  3Hz. 
As  is  shown  in  the  Fourier  spectra  shown  in  Figure  15,  a  significant  portion  of  the 
structural  response  occurs  in  the  5-10Hz  frequency  spectrum.  Figure  20  shows  that  with 
the  cut-off  set  at  3Hz  the  two-pole  Bessel  filter  actually  filters  frequencies  well  above 
lOFIz.  While  not  easily  apparent  in  this  comparison,  this  results  in  a  significant  loss  of 
magnitude  in  the  structural  velocity  response,  particularly  after  50  msec.  Lowering  the 
cut-off  frequency  to  0.5Hz  results  in  less  unwanted  filtering,  but  as  with  drift 
compensation,  evidence  of  rigid  body  motion  still  exist  afterwards. 


Bessel  High  Pass  Filter 


Frequency  (Hz) 


Figure  14.  Two-Pole  High-Pass  Bessel  Filter  With  3Hz  Cut-off  Frequency. 
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"Vertical  "Velocity 


Figure  15.  Fourier  Spectra  Plot  Showing  Structural  Low  Frequency  Response 

(Vertical  Velocity  Input) 


High  Pass  Filter  Comparison  (3Hz  Cut-off) 
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A  third-order  high-pass  Butterworth  filter  has  a  preferable  high-pass 
characteristic.  As  shown  by  Figure  16,  a  third-order  Butterworth  filter  reduces  the 
amount  of  desirable  frequency  response  that  is  filtered,  but  still  has  a  small  effect  on  the 
low  frequency  response.  Figure  17  shows  the  result  after  applying  a  0.5Hz  cut-off  to  the 
data  using  a  Butterworth  filter.  As  with  the  0.5Hz  cut-off  Bessel  filter,  a  small  amount  of 
low  frequency  drift  is  still  apparent. 


Vertical  Velocity 


Time  (msec) 

Figure  17.  Velocity  Response  After  Filtering  0.5Hz  High-Pass  Butterworth  Filter 

(Third-Order) 

Previous  UNDEX  research  has  utilized  only  drift  compensated  data  and  it  is 
currently  the  recognized  method  of  shock  trial  data  analysis.  Simulation  comparisons 
made  using  filtered  data  versus  using  drift  compensated/polynomial  fitted  data  were 
comparable.  While  the  majority  of  the  comparisons  made  in  this  report  use  drift 
compensated  data,  high-pass  filtering  the  data  offers  an  easier  alternative  worthy  of 
consideration  in  the  future. 
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B.  HIGH  FREQUENCY  “NOISE” 


In  addition  to  drift,  raw  shock  trial  data  has  a  fair  amount  of  high  frequency 
“noise”  which  clutters  the  data.  The  high  frequency  response  is  not  of  concern  in 
UNDEX.  These  high  frequencies  are  filtered  out  using  a  250Hz  cut-off  low-pass  filter. 
Figure  18  shows  an  example  of  the  effects  this  filtering  has  on  the  data. 
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Figure  18.  Effect  of  Low-Pass  Filtering  Data 


While  applying  this  low-pass  filter  to  trial  data  has  been  widely  accepted,  some 
debate  has  risen  over  the  validity  of  applying  the  same  low-pass  filter  to  simulation  data. 
To  answer  this,  a  statistical  study  was  conducted  using  223  accelerometer  measurements 
low-pass  filtered  at  250Hz.  These  were  compared  with  33  nodal  acceleration  predictions 
from  a  simulation  of  the  same  shot.  The  simulation  data  was  compared  before  and  after 
being  low-pass  filtered.  The  results  of  this  study  showed  an  unfiltered  simulation 
prediction  mean  much  higher  than  the  measured  values  and  an  excessively  large 
variation.  The  filtered  simulation  predictions  not  only  displayed  a  more  accurate  mean, 
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but  had  a  more  reasonable  variation  as  well.  As  recommended  by  this  study,  all 
simulation  response  predictions  compared  in  this  thesis  were  low-pass  fdtered  using  the 
same  cut-off  frequency  as  was  used  with  the  measured  data.  A  summary  of  the  results  of 
this  study  is  shown  in  Table  5. 

Table  5.  Statistical  Study  of  Low-Pass  Filtered  vs.  Unfiltered  Simulation  Data 


Shock  Trial  Data 
(Filtered) 

Simulation  Data 
(Unfdtered) 

Simulation  Data 
(Filtered) 

Mean 

26.225 

82.985 

34.297 

Variance 

520.229 

5775.711 

606.426 

Standard  Deviation 

22.809 

75.998 

24.626 
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V.  SIMULATION  RESULTS  AND  ANALYSIS 


Three  underwater  explosions,  shots,  were  conducted  during  the  DDG  81  shock 
trials.  The  most  severe  of  these  shots  was  the  third  shot,  “Shot  3”,  while  the  least  severe 
shot  was  the  first  shot,  “Shot  1”.  All  three  of  these  shots  have  been  simulated  using  the 
structural  model  described  in  Chapter  III  coupled  with  several  fluid  models.  The  primary 
focus  of  this  report  is  to  investigate  the  effects  of  fluid  mesh  density,  the  number  of 
elements  between  the  outer  DAA  boundary  and  the  fluid- structure  interface,  and  the 
ability  to  accurately  simulate  different  shot  geometries.  In  addition  to  this,  the  effects  of 
varying  mesh  quality  and  fluid  volume  depth  are  investigated. 


A.  ANALYSIS  AND  COMPARISON 

1.  Sensor/Node  Location 

Over  600  sensors  were  mounted  onboard  DDG  8 1  to  measure  the  response  during 
each  shot  conducted.  These  sensors  consisted  of  accelerometers,  velocity  meters  and 
strain  gauges.  For  comparison  purposes,  corresponding  nodes  were  built  into  the 
structural  model  representing  several  of  these  sensors.  These  nodes  were  defined  in  LS- 
DYNA  allowing  a  transient  response  time  history  to  be  kept  for  each  one.  This  time 
history  was  then  downloaded  into  GLview  and  subsequently  imported  into  UERDTools 
for  comparison  with  the  measured  sensor  response. 

During  an  underwater  shock  event,  the  primary  structural  response  is  in  the 
vertical  direction.  For  this  reason,  the  vertical  velocities  were  used  primarily  for  the 
comparisons  discussed  in  this  paper.  For  vertical  velocity  comparisons,  32  sensors  were 
compared  for  Shot  3,  30  for  Shot  2,  and  28  for  Shot  1.  While  the  majority  of  these 
sensors  were  vertical  velocity  meters,  several  accelerometers  were  used  as  well.  The 
accelerometer  data  was  subsequently  integrated  for  comparison  with  predicted  velocity 
response.  A  list  of  all  of  the  sensors  used  in  the  comparisons  is  shown  in  Table  6.  Those 
sensors  beginning  with  an  “A”  are  accelerometers  while  those  beginning  with  a  “V”  are 
velocity  meters. 
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Table  6.  Vertical  Velocity  Response  Sensor  Location  (N/A  =  no  data  available) 


Sensor 

NODE 

X(Node) 

Y(Node) 

Z(Node) 

Compartment/Location 

Compt#/Area 

Shot  1 

Shot  2 

Shot  3 

V2000V 

120217 

5328.0 

0.0 

82.0 

PASSAGE 

4-022-0-L 

X 

X 

X 

V2002V 

142489 

4656.0 

27.0 

82.0 

SONAR  COOLING  EQ  RM 

4-042-0-L 

X 

X 

X 

A4701V 

211371 

4080.0 

0.0 

280.0 

CREW  LVG  SPACE  #2 

3-097-2 -L 

N/A 

N/A 

X 

A8516V 

210993 

4080.0 

238.5 

202.0 

CREW  LVG  SPACE  #2 

3 -097-2 -L 

X 

X 

X 

V2107V 

211904 

4080.0 

-312.0 

390.0 

PASSAGE 

1-078-1 -L 

N/A 

N/A 

X 

V2108V 

212196 

4080.0 

312.0 

390.0 

PASSAGE 

1-078-1-L 

N/A 

X 

X 

V2007V 

210430 

4080.0 

0.0 

82.0 

AMR  #1 

4-126-0-E 

X 

X 

X 

V2009V] 

210808 

4080.0 

-174.0 

177.0 

AMR  #1 

4-126-0-E 

X 

X 

X 

V2008V1 

210894 

4080.0 

174.0 

177.0 

AMR  #1 

4-126-0-E 

X 

X 

X 

A2109V 

414953 

4059.1 

0.0 

722.8 

RR  #1 

03-128-0-C 

X 

X 

X 

V2011VJ 

221188 

3504.0 

216.0 

177.0 

AMR  #1 

4-126-0-E 

X 

X 

X 

V2010V 

220589 

3504.0 

0.0 

85.5 

AMR  #1 

4-126-0-E 

X 

X 

X 

V2012VJ 

221102 

3504.0 

-216.0 

177.0 

AMR  #1 

4-126-0-E 

X 

X 

X 

V2124V 

222060 

3504.0 

-375.4 

390.0 

PASSAGE 

1-158-1-L 

X 

X 

X 

V2125V 

222436 

3504.0 

375.4 

390.0 

PASSAGE 

1-158-8-L 

N/A 

X 

X 

V2013V 

221601 

3504.0 

0.0 

280.0 

ENGINE  RM  #1 

4-174-0-E 

X 

X 

X 

A2116V 

414367 

3504.0 

0.0 

702.0 

RR  #2 

03-142 -0-C 

X 

X 

X 

A2237V 

416419 

3504.0 

135.0 

848.0 

PORT  MAST  LEG 

MAST 

X 

X 

X 

A2240V 

416269 

3504.0 

-135.0 

848.0 

STBD  MAST  LEG 

MAST 

X 

X 

X 

A2104V 

222240 

3504.0 

0.0 

390.0 

CIC  ANNEX 

1- 126-0 -C 

X 

X 

X 

V2014V 

230461 

2952.0 

0.0 

85.5 

AMR  #2 

4-220-0-E 

X 

X 

X 

A3565V 

231696 

2952.0 

-81.0 

316.7 

ACCESS  TRUNK 

3-220-0-T 

X 

X 

X 

V2016V 

242399 

2544.0 

0.0 

116.0 

ENGINE  RM  #2 

4-254-0-E 

X 

X 

X 

V2026V 

312302 

1992.0 

0.0 

55.0 

ESCAPE  TRUNK 

4-296- 1-T 

X 

X 

X 

A2310V 

320746 

1536.0 

0.0 

177.0 

A/C  MCHY  &  PUMP  RM 

5-300-01-E 

X 

X 

X 

V2034V 

330759 

1152.0 

-135.0 

193.3 

FAN  ROOM 

3-362-0-E 

X 

X 

X 

V2032V 

330764 

1152.0 

0.0 

196.0 

FAN  ROOM 

3-362-0-E 

X 

X 

X 

V2035V 

330769 

1152.0 

135.0 

193.3 

FAN  ROOM 

3-362-0-E 

X 

X 

X 

V2019V 

340167 

672.0 

0.0 

196.5 

GENERATOR  RM 

3-370-0-E 

X 

X 

X 

V2018V1 

340992 

672.0 

0.0 

364.8 

PASSAGE 

2-410-0-L 

X 

X 

X 

A2413V 

350220 

288.0 

0.0 

273.4 

RAST  MACHINERY  RM 

2-442 -0-E 

X 

X 

X 

V2020V 

350052 

288.0 

0.0 

211.0 

STEERING  GEAR  RM 

4-442 -0-E 

X 

X 

X 

2.  Error  Measurement 

Quantifying  how  well  a  calculated  transient  response  compares  to  a  measured 
response  is  very  subjective.  One  way  to  eliminate  any  bias  from  the  comparison  is  by 
using  an  impartial  error  measure  such  as  Russell’s  error  factor.  Russell’s  error  factor 
evaluates  the  magnitude  and  phase  errors  separately,  combining  the  two  to  form  a  single 
comprehensive  error  factor  [Ref.  23]. 

The  relative  magnitude  error  is  found  by  first  forming  two  variables, 

N 

A=Ifi(i)2 

i=l 
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(12) 


and 


B 


Zf2(i)2 


i=l 


(13) 


where  fi(i)  and  f2(i)  are  the  measured  and  predicted  response  magnitudes 
respectively  at  a  given  time  step,  i.  The  two  variables  can  then  be  used  to  find  the 
relative  magnitude  error, 


m  = 


(A-B) 

Vab 


(14) 


The  phase  error  contributes  significantly  to  the  error  between  two  transient 
responses.  If  (f>  is  allowed  to  be  the  nonnalized  unit  vector  of  a  transient  response,  the 
phase  correlation  between  the  measured  and  predicted  responses  can  be  defined  by, 

/v  /v 

p  =  <h'<k  <15> 

Since  the  unit  vectors  are  nonnalized,  the  values  of  p  can  range  from  -1.0  to  1.0 
where  -1.0  indicates  that  the  two  responses  are  completely  out  of  phase,  while  1.0 
indicates  that  they  are  completely  in  phase. 

A  measure  of  the  phasing  between  two  transient  response  vectors  in  terms  of 
conelation  can  be  found  by  defining  a  new  tenn, 

N 

c  =  If 

i=l 

The  phase  correlation  between  the  two  responses  can  then  be  computed  by 

c 

p  -Jab 


(i)f2(i) 


(16) 
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It  is  important  to  note  that  p  represents  the  phasing  correlation  between  the  two 
responses,  it  is  not  a  measure  of  phase  error.  To  calculate  phase  error,  the  following 
equation  is  used: 


RP  = 


(18) 


The  phase  error  factor  has  an  error  range  of  0.0  to  1.0  where  0.0  indicates  both 
responses  are  completely  in  phase  while  1 .0  indicates  they  are  completely  out  of  phase. 


Although  the  phase  error  factor  has  a  maximum  value  of  1.0,  the  relative 
magnitude  error  factor  is  unbounded.  Since  the  two  are  combined  to  form  the 
comprehensive  error,  it  is  easy  to  see  that  the  magnitude  error  could  easily  dominate  the 
comprehensive  error,  presenting  an  undesirable  bias.  To  apply  a  similar  bound  to  the 
magnitude  error  factor  the  following  magnitude  error  factor  is  defined: 


RM 


sign(m)log10(l+  m ) 


(19) 


This  maintains  the  sign  unbiased  nature  of  m  while  effectively  artificially 
bounding  the  magnitude  error  factor  since  a  RM  value  of  1.0  represents  an  order  of 
magnitude  error  between  the  two  responses.  The  comprehensive  error  factor  can  now  be 
determined  using  the  following  equation: 

RC  = 


-(RM2  +RP2 
41 


(20) 


where  the  —  tenn  is  a  scale  factor  found  by  calculating  the  area  of  a  square  with  a 
width  of  length  RM  and  height  of  length  RP.  A  circle  with  a  corresponding  area  has  a 
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radius  equal  to  —  times  the  diagonal  of  the  square  [Ref.  23].  The  comprehensive  error 

factor  is  not  bounded,  but  errors  in  excess  of  1.0  indicate  substantial  error  between  data 
sets  and  virtually  no  correlation. 

Russell’s  error  factor  allows  an  unbiased  error  value  to  be  assigned  to  the 

correlation  between  the  measured  and  predicted  transient  responses.  The  determination 
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as  to  what  justifies  an  acceptable  or  unacceptable  prediction  is  more  subjective  however. 
It  has  been  suggested  that  the  values  outlined  in  Table  7  should  be  used  for  acceptance 
criteria  of  a  500  msec  processed  velocity  data  comparison  [Ref.  24].  Processed  data  in 
this  sense  means  both  measured  and  predicted  response  data  low  pass  filtered  with  a  cut¬ 
off  frequency  of  250  Hz  and  measured  data  drift  compensated  to  remove  gauge  drift. 
This  acceptance  criteria  is  not  valid  for  shorter  duration  comparisons  or  data  which  has 
been  processed  using  other  techniques  [Ref.  25], 


Table  7.  Russell’s  Comprehensive  Error  Factor  Acceptance  Criteria 


RC  <  0.15 

EXCELLENT 

0.15  <RC<  0.28 

ACCEPTABLE 

RC>0.28 

POOR 

It  should  be  noted  that  this  acceptance  criteria  is  only  an  agreed  upon  means  of 
determining  an  acceptable  or  unacceptable  prediction.  While  these  error  factors  present  a 
useful  tool  for  statistical  comparison  between  models,  there  are  no  magical  lines  of 
distinction  drawn  at  the  levels  of  acceptance.  This  is  to  say  that  a  comparison  resulting  in 
a  comprehensive  error  factor  of  0.148  placing  it  in  the  “Excellent”  category  is  not 
significantly  better  than  an  “Acceptable”  comparison  resulting  in  a  comprehensive  error 
factor  of  0.152.  Similarly,  though  unacceptable,  a  comparison  resulting  in  a  0.285  error 
factor  is  not  significantly  worse  than  one  whose  comprehensive  error  factor  is  0.275,  an 
“Acceptable”  value.  With  this  in  mind,  comparisons  should  be  summarized  using  the 
actual  error  factor  values,  not  the  subjective  ratings  applied  to  those  numbers. 

B.  SEVERE  SHOCK  (SHOT  3) 

The  most  severe  shock  conducted  on  DDG  81  was  Shot  3.  A  diagram  showing 
the  relative  aspect  of  this  shot  to  the  ship  is  shown  in  Figure  19. 
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Figure  19.  Shot  3  Geometry.  Most  Severe  Shot,  Charge  Off  Port  Midships 
1.  Surrounding  Fluid  Depth  Study 

The  initial  focus  of  the  investigations  discussed  in  this  report  was  to  determine  the 
effects  of  surrounding  fluid  depth  on  the  accuracy  of  the  simulation  predictions  [Ref.  26]. 
Four  models  were  developed,  a  15  inch  thick  two  layer  inner  liner  fluid  mesh,  a  half¬ 
cavitation  depth  fluid  mesh,  a  cavitation  depth  fluid  mesh,  and  a  fluid  mesh  having  a 
depth  twice  the  cavitation  depth.  This  comparison  used  Shot  3  exclusively  as  it  was  the 
most  severe  shock,  producing  the  greatest  response  magnitudes. 

a.  Two-Layer  Inner  Liner  Fluid  Mesh 

The  two-layer  inner  liner  fluid  mesh  shown  in  Figure  20  is  the  same  fluid 
mesh  shown  previously  in  Figure  6.  It  consists  of  a  seven  inch  thick  inner  layer  which  is 
encased  in  an  eight  inch  thick  outer  layer.  This  fluid  mesh  constitutes  the  smallest  fluid 
model  used  in  the  simulations.  Table  8  shows  the  Russell’s  error  factors  when  comparing 
the  simulation  predictions  against  measured  data  over  the  first  120  msec  of  response  time. 
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Figure  20.  Fluid  Model  Two  Layer  Inner  Liner 
Table  8.  Fluid  Mesh  Liner  120msec  Russell’s  Error  Comparison  [from  Ref.  26] 


Sensor 

RM 

RP 

RC 

V2000V 

-0.3284 

0.2920 

0.3894 

V2002V 

-0.0953 

0.3096 

0.2871 

V2007V 

-0.2074 

0.2989 

0.3224 

V2008VI 

-0.1413 

0.3070 

0.2995 

V2009VI 

-0.2038 

0.2996 

0.3211 

V2010V 

-0.1796 

0.2776 

0.2930 

V2011VI 

-0.2351 

0.2723 

0.3188 

V2012VI 

-0.2193 

0.3052 

0.3331 

V2013V 

-0.2517 

0.2790 

0.3330 

V2014V 

0.0409 

0.2803 

0.2510 

V2016V 

-0.0838 

0.2462 

0.2306 

V2018VI 

-0.2055 

0.2405 

0.2803 

V2019V 

-0.2193 

0.2787 

0.3143 

V2020V 

-0.2213 

0.2682 

0.3081 

V2026V 

-0.0328 

0.2810 

0.2508 

V2032V 

-0.1792 

0.2799 

0.2945 

V2034V 

-0.1850 

0.2599 

0.2828 

V2035V 

-0.1208 

0.3076 

0.2929 

V2108V 

-0.2230 

0.3371 

0.3582 

V2124V 

-0.1629 

0.3625 

0.3523 

V2125V 

-0.2443 

0.3068 

0.3476 
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It  should  be  noted  that  in  this  comparison  the  measured  data  was  high-pass 
fdtered  using  a  2-pole  Bessel  fdter  set  at  a  cut-off  frequency  of  3  Hz  while  the  predicted 
data  was  high-pass  filtered  using  the  same  filter  but  with  a  cut-off  frequency  of  1.8  Hz. 
Additionally,  the  focus  of  this  investigation  was  to  detennine  how  well  each  fluid  model 
captured  the  initial  peak  response.  This  early  time  response  is  captured  completely 
within  the  first  120  msec  after  arrival  of  the  incident  pressure  wave. 

As  shown  in  Table  8,  the  inner  liner  fluid  model  does  not  accurately  predict  the 
structural  response  to  the  shock  wave.  Over  half  of  the  21  sensors  compared  had  a 
comprehensive  error  factor  over  0.30  and  none  had  comprehensive  error  factors  below 
0.20.  The  errors  were  significant  in  both  phase  and  magnitude,  though  the  phase  error 
did  dominate  the  majority  of  the  comparisons.  In  all  but  one  case  the  simulation 
predictions  underestimated  the  response  of  the  structure. 

b.  One-Half  Cavitation  Depth  Fluid  Mesh 

The  second  fluid  model  compared  extended  the  depth  of  the  fluid  out  to 
approximately  one-half  of  the  cavitation  depth.  This  model  is  shown  in  Figure  21. 


IUII. 
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Table  9  lists  the  Russell’s  error  factors  using  this  model  in  the  simulation. 
As  with  the  inner  liner  fluid  mesh  comparison,  both  data  sets  were  high-pass  filtered  and 
the  comparison  is  for  the  initial  120  msec  of  the  response. 


Table  9.  One-Half  Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error  Comparison 

[from  Ref.  26] 


Sensor 

RM 

RP 

RC 

V2000V 

-0.3978 

0.2966 

0.4397 

V2002V 

-0.8131 

0.3959 

0.8015 

V2007V 

-0.5015 

0.1983 

0.4779 

V2008VI 

-0.5190 

0.4017 

0.5816 

V2009VI 

-0.4634 

0.1357 

0.4279 

V2010V 

-0.2777 

0.1873 

0.2969 

V2011VI 

-0.2966 

0.1386 

0.2901 

V2012VI 

-0.3243 

0.1026 

0.3015 

V2013V 

-0.2881 

0.1473 

0.2868 

V2014V 

0.0663 

0.1722 

0.1635 

V2016V 

-0.0858 

0.1228 

0.1328 

V2018VI 

-0.0807 

0.2748 

0.2538 

V2019V 

0.0459 

0.1796 

0.1643 

V2020V 

0.1331 

0.1877 

0.2039 

V2026V 

-0.0380 

0.1802 

0.1632 

V2032V 

0.0228 

0.1198 

0.1081 

V2034V 

0.0922 

0.0881 

0.1130 

V2035V 

0.0382 

0.1219 

0.1132 

V2108V 

-0.5072 

0.3714 

0.5571 

V2124V 

-0.3295 

0.0879 

0.3023 

V2125V 

-0.3143 

0.1457 

0.3070 

The  comparisons  of  this  fluid  model  simulation  show  significant 
improvement  over  the  inner  liner  fluid  model  simulation  comparison.  Significant  errors 
still  exist  in  the  comparison  with  9  of  the  2 1  comparisons  having  a  comprehensive  error 
factor  over  0.30.  A  third  of  the  comparisons  however  had  a  comprehensive  error  factor 
less  than  0.20,  over  half  of  these  below  0.15.  Though  better  than  the  inner  liner  fluid 
model  simulation  predictions  for  the  most  part,  using  this  fluid  model  in  the  simulation 
does  not  result  in  predictions  which  are  consistently  accurate  as  Table  9  suggests. 
c.  Cavitation  Depth  Fluid  Mesh 

The  third  fluid  model  compared  had  a  depth  approximately  equal  to  the 
cavitation  depth  of  the  actual  shot.  Although  the  outer  boundary  of  the  fluid  model  is 
perfectly  cylindrical,  the  cavitation  depth  changes  along  the  length  of  the  ship.  Due  to 
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this,  a  perfectly  cylindrical  fluid  model  can  only  approximate  the  cavitation  depth  and,  in 
this  case,  was  modeled  having  a  sufficient  depth  to  include  the  cavitation  region 
throughout  the  length  of  the  ship.  This  fluid  model  is  shown  in  Figure  22. 


i 

140  ft. 


Figure  22.  Cavitation  Depth  Fluid  Mesh  [from  Ref.  26] 

The  120  msec  comparison  using  this  fluid  model  in  the  simulation  is  shown  in 
Table  10.  As  with  the  other  fluid  model  comparisons,  all  data  has  been  high-pass  filtered 
in  these  comparisons.  These  comparisons  show  substantial  improvement  over  the 
previous  two  simulation  prediction  comparisons.  Only  one  comparison  has  a 
comprehensive  error  factor  greater  than  0.30  and  19  of  the  21  comparisons  have 
comprehensive  error  factors  less  than  0.20,  14  below  0.15,  and  4  below  0.10.  The 
majority  of  the  error  is  due  to  the  phase  error  as  the  majority  of  the  magnitude  errors  are 
very  small.  From  these  comparisons,  it  is  clear  that  the  fluid  model  should  be  extended 
out  to  the  cavitation  depth  to  accurately  predict  the  initial  120  msec  response  of  the 
structure. 
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Table  10.  Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error  Comparison  [from 

Ref.  26] 


Sensor 

RM 

RP 

RC 

V2000V 

-0.2445 

0.2909 

0.3368 

V2002V 

0.0565 

0.1370 

0.1314 

V2007V 

-0.0988 

0.1478 

0.1576 

V2008VI 

-0.0134 

0.1580 

0.1405 

V2009VI 

0.0060 

0.1084 

0.0962 

V2010V 

-0.0447 

0.1199 

0.1134 

V2011VI 

-0.0310 

0.1204 

0.1102 

V2012VI 

0.0100 

0.0848 

0.0756 

V2013V 

0.0303 

0.0854 

0.0803 

V2014V 

0.0656 

0.1526 

0.1472 

V2016V 

-0.0077 

0.1152 

0.1023 

V2018VI 

-0.1336 

0.1525 

0.1796 

V2019V 

-0.1872 

0.1532 

0.2144 

V2020V 

-0.1024 

0.1437 

0.1564 

V2026V 

-0.0789 

0.1255 

0.1314 

V2032V 

-0.1261 

0.1562 

0.1779 

V2034V 

-0.1006 

0.1092 

0.1316 

V2035V 

-0.1153 

0.1272 

0.1521 

V2108V 

0.0109 

0.1262 

0.1123 

V2124V 

-0.0233 

0.0992 

0.0903 

V2125V 

-0.0596 

0.1222 

0.1205 

d.  2X-Cavitation  Depth  Fluid  Mesh 

The  exceptional  correlation  using  the  cavitation  depth  fluid  mesh  and  the 
poor  correlation  using  the  smaller  fluid  mesh  models  indicates  that  there  is  at  least  a 
minimum  depth  constraint  for  the  fluid  model.  To  determine  if  there  is  a  maximum  depth 
constraint  as  well,  a  fluid  model  having  a  depth  equal  to  approximately  twice  the 
cavitation  depth  was  used  to  conduct  a  simulation.  This  model  is  shown  in  Figure  23. 

Table  1 1  shows  the  resulting  error  factors  from  the  simulations  conducted 
with  this  model.  This  comparison  shows  the  worst  overall  correlation  of  all  of  the  models 
used  in  the  simulations.  Nineteen  of  the  21  sensors  compared  had  comprehensive  errors 
over  0.30  for  this  model.  Of  interest  in  this  comparison  however  is  that  the 
comprehensive  error  is  not  dominated  by  the  phase  error  as  in  the  previous  three 
comparisons,  but  in  the  magnitude  error.  This  is  indicative  of  incident  pressure  wave 
losses  occurring  in  the  model.  This  is  investigated  further  later  on  in  this  chapter. 
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Table  11. 


ft. 


i _ i 

280  ft. 

Figure  23.  2X-Cavitation  Depth  Fluid  Mesh  [from  Ref.  26] 


2X-Cavitation  Depth  Fluid  Mesh  120msec  Russell’s  Error  Comparison  [from 

Ref.  26] 


Sensor 

RM 

RP 

RC 

V2000V 

-0.6022 

0.3148 

0.6022 

V2002V 

-0.3407 

0.1664 

0.3360 

V2007V 

-0.4032 

0.2202 

0.4071 

V2008VI 

-0.3683 

0.2736 

0.4066 

V2009VI 

-0.3466 

0.1847 

0.3481 

V2010V 

-0.3472 

0.1959 

0.3533 

V2011VI 

-0.3287 

0.1901 

0.3365 

V2012VI 

-0.3097 

0.1485 

0.3044 

V2013V 

-0.3048 

0.1543 

0.3027 

V2014V 

-0.2080 

0.1621 

0.2337 

V2016V 

-0.2350 

0.1687 

0.2564 

V2018VI 

-0.3806 

0.2692 

0.4132 

V2019V 

-0.3924 

0.2887 

0.4317 

V2020V 

-0.4265 

0.2606 

0.4430 

V2026V 

-0.2738 

0.2147 

0.3084 

V2032V 

-0.4255 

0.1587 

0.4024 

V2034V 

-0.3917 

0.1392 

0.3684 

V2035V 

-0.4147 

0.1442 

0.3891 

V2108V 

-0.3856 

0.1992 

0.3846 

V2124V 

-0.3149 

0.1734 

0.3186 

V2125V 

-0.3355 

0.2027 

0.3474 
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2.  High  Pass  Filtered  Data  vs.  Drift  Compensated  Data 

In  the  previous  surrounding  fluid  model  comparisons  the  Russell’s  error  factors 
were  detennined  after  high-pass  filtering  the  measured  data  using  a  cut-off  frequency  of  3 
Hz  and  filtering  the  predicted  response  data  using  a  high-pass  filter  cut-off  frequency  of 
1.8  Hz  This  filtering  was  performed  using  a  two-pole  Bessel  filter.  The  resulting 
Russell’s  error  factors  found  using  this  data  processing  technique  for  the  cavitation  depth 
model  comparison  of  420  msec  of  response  time  are  shown  in  Table  12. 


Table  12.  Russell’s  Error  Factor  Using  High  Pass  Filtered  Data  (420  msec)  [from  Ref. 

26] 


Sensor 

RM 

RP 

RC 

V2013V 

0.0263 

0.1329 

0.1200 

V2012VI 

0.0583 

0.1424 

0.1363 

V2108V 

0.0038 

0.1543 

0.1368 

V2124V 

0.0386 

0.1562 

0.1425 

V2016V 

0.0291 

0.1592 

0.1434 

V2009VI 

0.1113 

0.1388 

0.1577 

V2034V 

0.0173 

0.1787 

0.1591 

V2026V 

0.0542 

0.1849 

0.1707 

V2014V 

0.0372 

0.1949 

0.1759 

V2011VI 

0.0906 

0.1972 

0.1923  i 

V2002V 

0.1240 

0.1873 

0.1991  ! 

V2035V 

0.0513 

0.2243 

0.2039 

V2018VI 

0.1039 

0.2055 

0.2041 

V2008VI 

0.1290 

0.1937 

0.2063 

V2125V 

0.0488 

0.2334 

0.2113 

V2010V 

0.0449 

0.2456 

0.2213 

V2032V 

0.1399 

0.2083 

0.2224 

V2020V 

0.1334 

0.2165 

0.2254 

V2019V 

0.1059 

0.2331 

0.2269 

V2007V 

0.1464 

0.2122 

0.2285  i 

V2000V 

0.2419 

0.3164 

0.3529 

The  same  predicted  response  data  were  used  to  find  the  Russell’s  error  factors 
summarized  in  Table  13  where  no  high-pass  filtering  has  been  applied.  The  measured 
response  data  has  instead  been  corrected  using  the  drift  compensation  method  described 
previously  in  this  chapter.  As  shown  in  the  two  comparisons,  the  filtered  data  yields 
artificially  low  Russell’s  comprehensive  error  factors,  primarily  due  to  the  magnitude 
error  factor.  This  is  consistent  with  the  discussion  in  the  previous  chapter  where  filtering 
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results  in  a  significant  reduction  in  the  response  magnitude.  This  reduction  tends  to 
reduce  the  error  between  the  two  responses  as  indicated  by  the  reduction  in  the  magnitude 
error  factor.  As  discussed,  the  Russell’s  error  factor  acceptance  criteria  is  only  valid  for 
velocity  response  comparison  using  drift  compensated  measured  data.  For  that  reason  all 
further  comparisons  highlighted  in  this  report  were  conducted  using  drift  compensated 
data. 


Table  13.  Russell’s  Error  Factor  Using  Drift  Compensated  Data  (500  msec) 


Sensor 

RM 

RP 

RC 

V2013V 

-0.1123 

0.1859 

0.1925 

V2124V 

-0.1543 

0.1716 

0.2045 

V2012VI 

-0.1435 

0.1864 

0.2085 

V2034V 

0.1857 

0.1450 

0.2088 

V2016V 

-0.1481 

0.1948 

0.2169 

V2002V 

0.0310 

0.2466 

0.2203 

V2014V 

-0.1086 

0.2414 

0.2346 

V2018VI 

0.1306 

0.2304 

0.2347 

V2108V 

-0.1372 

0.2274 

0.2354 

V2026V 

-0.1723 

0.2219 

0.2490 

V2008VI 

-0.1394 

0.2476 

0.2519 

V2009VI 

-0.1511 

0.2543 

0.2621 

V2019V 

0.1780 

0.2420 

0.2662 

V2125V 

-0.1982 

0.2277 

0.2675 

V2035V 

0.2364 

0.2082 

0.2791 

V2032V 

-0.2280 

0.2200 

0.2807 

V2011VI 

-0.2121 

0.2436 

0.2863 

V2010V 

-0.2358 

0.2425 

0.2997 

V2007V 

-0.2296 

0.2676 

0.3125 

V2020V 

0.2165 

0.2843 

0.3167 

V2000V 

0.2891 

0.3109 

0.3763 

3.  Fluid  Mesh  Depth  Sensitivity 

The  cavitation  depth  models  used  in  the  previous  simulations  were  constructed 
having  a  depth  of  840  inches  as  a  rough  approximation  to  the  actual  cavitation  depth. 
This  was  to  account  for  the  greatest  cavitation  depth,  that  depth  at  a  point  on  the  ship’s 
keel  furthest  from  the  charge.  The  cavitation  depth  at  the  point  on  the  keel  closest  to  the 
charge  was  found  to  be  only  738  inches  for  this  shot.  Simulations  conducted  using  a 
similar  model,  but  with  a  reduced  depth,  resulted  in  a  decrease  of  over  0.01  in  the  average 
Russell’s  comprehensive  error  factor.  This  suggests  that  the  simulation  is  particularly 
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sensitive  to  small  changes  in  the  depth  of  the  fluid  modeled.  This  model  was 
subsequently  modified,  retaining  the  initial  inner  liner  and  majority  of  the  inner  mesh 
while  replacing  the  flat  end  caps  with  spherical  end  caps.  This  modification  eliminated 
the  need  for  multiple  parts  in  the  outer  fluid  mesh,  facilitating  more  efficient  modification 
of  the  outer  fluid  mesh,  and  improved  the  fluid  model’s  applicability  to  a  wider  variety  of 
underwater  explosion  geometries.  No  appreciable  change  in  the  average  error  factors 
occurred  from  this  modification.  This  2nd  generation  model  is  shown  in  Figure  24. 


Table  14  summarizes  the  best  correlation  achieved  using  the  model  shown  in 
Figure  24.  This  was  done  using  a  densely  packed  fluid  mesh  where  the  average  hull 
nonnal  spacing  between  elements  was  less  than  12  inches.  Compared  with  the  original 
1st  generation  model,  the  improvement  in  the  average  comprehensive  error  was  over  0.02, 
a  significant  reduction. 
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Table  14.  Summary  of  Prediction  Error  Measures  Using  Dense  2nd  Generation  Model 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Fluid  Mesh  (2nd 
Generation  Model) 

X 

n 

z 

RM 

RP 

RC 

V2000V 

120217 

Keel 

5328 

0 

0.2980 

0.3165 

V2002V 

142489 

Keel 

4656 

27 

0.2403 

0.2197 

V2007V 

210430 

Keel 

4080 

0 

82 

-0.1274 

0.2875 

0.2787 

V2009VI 

210808 

Bulkhead 

4080 

-174 

177 

-0.0793 

0.2145 

0.2026 

V2008VI 

210894 

Bulkhead 

4080 

174 

177 

-0.0124 

0.2623 

0.2327 

A8516V 

210993 

Bulkhead 

4080 

239 

0.1128 

0.2551 

0.2472 

A4701V 

211371 

Bulkhead 

4080 

0 

0.0080 

0.2042 

0.1811 

V2107V 

211904 

Bulkhead 

4080 

-312 

-0.0427 

0.2205 

0.199 

V2108V 

212196 

Bulkhead 

4080 

312 

-0.0344 

0.2394 

0.2143 

V2010V 

220589 

Keel 

3504 

0 

-0.1732 

0.2630 

0.2791 

V2012VI 

221102 

Bulkhead 

3504 

-216 

-0.0803 

0.2054 

0.1954 

V2011VI 

221188 

Bulkhead 

3504 

216 

-0.1594 

0.2454 

0.2593 

V2013V 

221601 

Bulkhead 

3504 

0 

-0.0366 

0.2088 

0.1878 

V2124V 

222060 

Bulkhead 

3504 

-375 

-0.1013 

0.1833 

0.1856 

A2104V 

222240 

Bulkhead 

3504 

0 

0.0109 

0.1797 

0.1596 

V2125V 

222436 

Bulkhead 

3504 

375 

-0.1453 

0.2280 

0.2396 

V2014V 

230461 

Keel 

2952 

0 

86 

-0.0992 

0.2565 

0.2438 

A3565V 

231696 

Bulkhead 

2952 

-81 

317 

0.0492 

0.1847 

0.1694 

V2016V 

242399 

Keel 

2544 

0 

116 

-0.1050 

0.2235 

0.2189 

V2026V 

312302 

Keel 

1992 

0 

55 

-0.1334 

0.2228 

0.2302 

A2310V 

320746 

Bulkhead 

1536 

0 

177 

-0.0970 

0.1496 

0.1580 

V2034V 

330759 

Keel 

1152 

-135 

193 

-0.0704 

0.1439 

0.1419 

V2032V 

330764 

Keel 

1152 

0 

196 

-0.1305 

0.2247 

0.2303 

V2035V 

330769 

Keel 

1152 

135 

193 

-0.1447 

0.2189 

0.2325 

V2019V 

340167 

Keel 

672 

0 

197 

0.0053 

0.2771 

0.2456 

V2018VI 

340992 

Bulkhead 

672 

0 

365 

0.0539 

0.2710 

0.2449 

V2020V* 

350052 

Keel 

288 

0 

211 

-0.0801 

0.2422 

0.2261 

A2413V 

350220 

Bulkhead 

288 

0 

273 

-0.0602 

0.2351 

0.2151 

A2116V 

414367 

Bulkhead 

3504 

0 

0.0175 

0.1749 

0.1557 

A2109V 

414953 

Bulkhead 

4059 

0 

0.1189 

0.1801 

0.1913 

A2240V 

416269 

Bulkhead 

3504 

-135 

848 

-0.0034 

0.1697 

0.1504 

A2237V 

416419 

Bulkhead 

3504 

135 

848 

0.0382 

0.1855 

0.1678 

Sum(E(X)) 

-1.63710 

7.09560 

6.82010 

Sum(E(XA2)) 

0.29793 

1.62116 

1.50719 

Mean 

-0.05116 

0.22174 

0.21313 

Variation 

0.00691 

0.00154 

0.00173 

*  - ; - 1  .  1  1 

350  msec  comparison  due  to  long  time  sensor  data  inaccuracy 


2.  Fluid  Mesh  Quality 

After  several  simulations  were  conducted  and  compared  using  the  2nd  generation 

model,  a  3ld  generation  model  was  developed  in  an  attempt  to  improve  the  simulation 

predictions.  Retaining  only  the  inner  liner  from  the  first  two  models,  this  model  focused 

on  eliminating  the  “trouble  spots”  on  the  1st  and  2nd  generation  models  where  the  fluid 
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mesh  was  extremely  non-uniform,  particularly  the  areas  surrounding  the  sonar  dome  and 
skeg  as  shown  in  Figures  25  and  26.  This  fluid  mesh  also  utilized  the  interpolation 
methods  discussed  in  Chapter  IV  more  heavily  than  the  first  two  models.  The  end  result 
is  that  this  model  had  an  improved  mesh  quality  over  the  previous  two  models. 


Figure  25.  Sonar  Dome  Region  Showing  Improvement  in  Mesh  Uniformity  of  3rd 

Generation  Fluid  Model. 


Figure  26.  Keel  Region  Showing  Improvement  in  Mesh  Uniformity  of  3  Generation 

Fluid  Model 
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The  comparison  of  the  2nd  generation  model  simulation  was  summarized  in  Table 
14.  Table  15  is  a  summary  of  the  comparable  3ld  generation  model  comparison.  A 
graphical  comparison  of  the  comprehensive  error  trend  lines  is  shown  with  respect  to 
location  on  the  ship  in  Figure  27. 


Table  15.  Summary  of  Simulation  Error  Measures  using  Dense  3rd  Generation  Model 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Fluid  Mesh  (3rd 
Generation  Model) 

x  y  z 

It  VI 

RP 

RC 

V2000V 

120217 

Keel 

5328  0  82 

-0.1795 

0.2842 

0.2979 

V2002V 

142489 

Keel 

4656  27  82 

0.0762 

0.2303 

0.2150 

V2007V 

210430 

Keel 

4080  0  82 

-0.1198 

0.2883 

0.2766 

V2009VI 

210808 

Bulkhead 

4080  -174  177 

-0.0690 

0.2144 

0.1996 

V2008VI 

210894 

Bulkhead 

4080  174  177 

0.0014 

0.2561 

0.2270 

A8516V 

210993 

Bulkhead 

4080  239  202 

0.1217 

0.2547 

0.2502 

A4701V 

211371 

Bulkhead 

4080  0  280 

0.0282 

0.2149 

0.1921 

V2107V 

211904 

Bulkhead 

4080  -312  390 

-0.0295 

0.2420 

0.2161 

V2108V 

212196 

Bulkhead 

4080  312  390 

-0.0216 

0.2348 

0.2089 

V2010V 

220589 

Keel 

3504  0  86 

0.1353 

0.1846 

0.2028 

V2012VI 

221102 

Bulkhead 

3504  -216  177 

-0.1729 

0.2549 

0.2730 

V2011VI 

221188 

Bulkhead 

3504  216  177 

-0.0777 

0.1935 

0.1848 

V2013V 

221601 

Bulkhead 

3504  0  280 

-0.1468 

0.2398 

0.2468 

V2124V 

222060 

Bulkhead 

3504  -375  390 

-0.0342 

0.1970 

0.1772 

A2104V 

222240 

Bulkhead 

3504  0  390 

-0.0872 

0.1882 

0.1838 

V2125V 

222436 

Bulkhead 

3504  375  390 

0.0207 

0.1760 

0.1571 

V2014V 

230461 

Keel 

2952  0  86 

-0.1264 

0.2264 

0.2298 

A3565V 

231696 

Bulkhead 

2952  -81  317 

0.0353 

0.1789 

0.1616 

V2016V 

242399 

Keel 

2544  0  116 

0.0122 

0.1747 

0.1552 

V2026V 

312302 

Keel 

1992  0  55 

0.0613 

0.1877 

0.1750 

A2310V 

320746 

Bulkhead 

1536  0  177 

-0.0861 

0.2501 

0.2344 

V2034V 

330759 

Keel 

1152  -135  193 

0.0646 

0.2003 

0.1865 

V2032V 

330764 

Keel 

1152  0  196 

-0.0879 

0.2085 

0.2005 

V2035V 

330769 

Keel 

1152  135  193 

-0.0910 

0.2312 

0.2202 

V2019V 

340167 

Keel 

672  0  197 

-0.1037 

0.1635 

0.1716 

V2018VI 

340992 

Bulkhead 

672  0  365 

-0.0966 

0.1537 

0.1608 

V2020V* 

350052 

Keel 

288  0  211 

-0.1534 

0.2103 

0.2307 

A2413V 

350220 

Bulkhead 

288  0  273 

-0.1598 

0.2233 

0.2434 

A2116V 

414367 

Bulkhead 

3504  0  702 

-0.1344 

0.2266 

0.2335 

A2109V 

414953 

Bulkhead 

4059  0  723 

-0.1212 

0.1830 

0.1945 

A2240V 

416269 

Bulkhead 

3504  -135  848 

-0.1499 

0.2178 

0.2343 

A2237V 

416419 

Bulkhead 

3504  135  848 

-0.1315 

0.2150 

0.2233 

Sum(E(X)) 

-1.82320 

6.90470 

6.76420 

Sum(E(XA2)) 

0.34947 

1.52398 

1.47020 

Mean 

-0.05698 

0.21577 

0.21138 

Variation 

0.00792 

0.00110 

0.00130 

350  msec  comparison  due  to  long  time  sensor  data  inaccuracy 


54 


Shot  3  Russell's  Error  vs.  Longitudinal  Position 


Figure  27.  Comparison  of  2nd  and  3rd  Generation  Fluid  Model  Errors  in  Relation  to 

Longitudinal  Position  on  the  Structural  Model. 

From  Tables  14  and  15  and  the  comparison  in  Figure  27,  the  results  from  this 
quality  comparison  are  inconclusive.  The  average  comprehensive  error  difference  is  less 
than  0.002,  a  relatively  insignificant  amount.  Notable  improvement  in  the  error  factor  is 
shown  in  the  regions  around  the  sonar  dome  and  the  stern  area,  all  areas  where  significant 
improvements  were  made  to  the  mesh.  However,  the  2nd  generation  model  had  smaller 
errors  in  other  regions  of  the  structure.  These  inconsistencies  indicate  that  other 
variables,  such  as  data  processing,  sensor  error,  and  environmental  factors,  have  more  of 
an  effect  on  simulation  error  factors  than  small  variations  in  the  quality  of  the  fluid  model 
mesh.  Figures  28  through  30  show  three  of  the  areas  with  the  largest  variations  in  error 
factors.  With  the  exception  of  the  sonar  dome  region,  where  the  3ld  generation  model 
does  a  better  job  of  capturing  the  initial  peak  magnitude,  the  primary  difference  lies  in  the 
late  time  response.  This  is  to  be  expected,  since  improving  the  mesh  quality  should  have 
an  affect  primarily  on  the  ability  to  simulate  cavitation  which  occurs  after  this  initial  peak 
response. 
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Velocity  (  ft/ see )  ^  Vertical  Velocity  (ft/sec) 


SHOT  3 

Grid  120217-vz  (V2000V) 


Keel  (x=5328  y=0  z=82) 


Time  (msec) 


Shock;  Trial  Data  - 3rd  Generation  Fluid  IVIodcl 

2nd  Generation  Fluid  Model 


igure  28. 


Sonar  Dome  Region  Velocity  Plot  Comparison 


SHOT  3 

Grid  23  1696  -vz  (A3565V) 


Bulkhead  (x=2952  y=-81  z=3  16.7) 


Time  (msec) 

Ship  Shock  Trial  Data  - 3rd  Generation  Fluid  Ivlodel 

2nd  Generation  Fluid  Pvlodel 


Figure  29.  Midship  Region  Velocity  Plot  Comparison 
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SHOT  3 

Grid  340992-vz  (V2018VI) 


Bulkhead  (x=672  y=0  z=364.8) 


Ship  Shock:  Trial  - 3rd  Generation  Fluid  iVTesh 

2nd  Generation  Fluid  Ivtesh 


Figure  30.  Stern  Region  Velocity  Plot  Comparison 

3.  Fluid  Mesh  Density 

The  accuracy  of  a  finite  element  model  (FEM)  in  predicting  the  dynamic  response 
of  a  structure  is  improved  by  increasing  the  number  elements  used  to  model  that 
structure.  As  illustrated  in  Figure  31,  this  improvement  in  accuracy  due  to  increasing  the 
number  of  elements  is  usually  quite  large  with  a  relatively  small  number  of  elements  and 
decreases  to  becoming  insignificant  with  a  relatively  large  number  of  elements.  It  would 
stand  to  reason  then  that  the  accuracy  of  the  fluid  mesh  FEM  would  be  improved  by 
increasing  the  number  of  elements  in  the  model.  For  purposes  of  this  discussion,  the 
number  of  elements  is  not  the  total  number  of  elements  in  the  model.  If  looking  at  a  two- 
dimensional  cross  sectional  cut  of  the  fluid-structure  model,  the  number  of  elements  in 
the  fluid  would  be  the  number  of  elements  in  each  mesh  line  extending  outward  in  a 
normal  direction  from  the  hull. 

Three  fluid  mesh  FEM’s  of  varying  densities  were  modeled  for  Shot  3.  All  of  the 
models  maintained  the  same  thickness  of  the  inner  two  layers.  The  inner  mesh  and 
transitional  mesh  were  varied  from  each  having  five  to  seven  elements  in  the  normal 
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Mode  5  FEM  Natural  Frequency  Error 


Figure  31.  Natural  Frequency  Error  of  FEM  Pin-Pin  Beam  as  a  Function  of  Number 

of  Elements  Used. 

direction.  The  outer  mesh  was  varied  from  having  as  few  as  8  elements  to  as  many  as  23 
elements  in  the  normal  direction.  Table  16  lists  a  summary  of  the  properties  for  each 
fluid  model  constructed  for  Shot  3. 

Table  16.  Shot  3  Fluid  Model  Properties 


Property 

Coarse  Mesh 

Medium  Mesh 

Dense  Mesh 

Inner  Mesh  Elements 

5 

6 

7 

Inner  Mesh  Element 
Spacing  (Ave.) 

11.4  in 

10.0  in 

9.3  in 

Transition  Mesh 
Elements 

5 

6 

7 

Transition  Mesh 
Element  Spacing  (Ave.) 

13.4  in 

1 1.6  in 

10.1  in 

Outer  Mesh  Elements 

8 

14 

23 

Outer  Mesh  Element 
Spacing  (Ave.) 

36  in 

20  in 

12  in 

Fluid  Nodes 

150,607 

197,429 

250,181 

DAA  Elements 

5048 

5090 

5132 

Simulation  Run  Time 

32.5  hrs 

34.1  hrs 

37.3  hrs 

To  avoid  a  sharp  transition  in  the  spacing  of  the  elements  in  the  coarse  and 
medium  fluid  meshes,  the  outer  mesh  element  spacing  was  graduated  from  smallest  on 
the  inner  boundary  to  largest  on  the  outer  boundary.  This  was  not  required  for  the  dense 
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mesh.  Table  16  lists  the  simulation  run  time.  As  shown,  it  took  nearly  five  additional 
hours  to  run  the  dense  model  than  the  coarse  model. 
a.  Error  Comparison 

The  Russell’s  error  factors  are  summarized  in  Figures  32  through  34. 


Russell's  Error  Factor  (Shot  3  Medium) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  Excellent 


Figure  32.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison 


Russell's  Error  Factor  (Shot  3  Medium) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  Excellent 


Figure  33.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison 
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Russell's  Error  Factor  (Shot  3  Dense) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  -  -  Excellent 

Figure  34.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison 

The  most  significant  change  occurs  between  the  dense  and  the  medium 
mesh  fluid  model  comparisons.  The  large  number  of  negative  magnitude  error  factors  is 
indicative  of  the  predicted  response  magnitude  being  lower  than  the  measured  response 
magnitude.  The  trend  when  increasing  the  density  of  the  fluid  mesh  is  that  the  errors  tend 
to  move  primarily  to  the  right  and  slightly  downward.  This  indicates  that  the  predicted 
response  magnitude  is  increasing  as  the  density  of  the  mesh  is  increased.  At  the  same 
time,  though  not  as  pronounced,  the  phase  error  is  decreasing  as  the  mesh  density  is 
increased. 

Table  17  lists  all  of  the  sensors  compared  and  their  corresponding 
Russell’s  error  factors.  The  averages  for  each  error  factor  are  shown  at  the  bottom  of 
their  respective  columns  along  with  the  variation.  These  averages  show  that  the  dense 
mesh  does  in  fact  have  lower  overall  magnitude,  phase,  and  comprehensive  error  factor. 
In  addition,  the  largest  reduction  in  these  error  factors  occurs  between  the  medium  and 
dense  fluid  mesh  comparisons  where  the  most  significant  improvement  is  in  the 
magnitude  error. 
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Table  17.  Shot  3:  Russell’s  Error  Factor. 


Sensor  NODE 


V2000V 

V2002V 

V2007V 

V2009VI 

V2008VI 

A8516V 

A4701V  ' 

V21Q7V 

V2108V 

V2010V 

V2012VI 

V2011VI 

V2013V 

V2124V 

A2104V 

V2125V 

V2014V 

A3565V 

V2016V 

V2026V 

A2310V 

V2034V 

V2032V 

V2035V 

V2019V 

V2018VI 

V2020V* 

A2413V 

A2116V 

A2109V 

A2240V 

A2237V 


120217 

142489 

210430 

210808 

210894 

210993 

211371 

211904 

212196 

220589 

221102 

221188 

221601 

222060 

222240 

222436 

230461 

231696 

242399 

312302 

320746 

330759 

330764 

330769 

340167 

340992 

350052 

350220 

414367 

414953 

416269 

416419 


Bulkhead 

Keel 


Ke 
Keel 
Keel 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Keel 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Keel 

Bulkhead 

Keel 

Keel 

Bulkhead 

Keel 

Keel 

Keel 

Keel 

Bulkhead 

Keel 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 

Bulkhead 


Location  (in) 

x  I  y  I  z 


5328 
4656  : 

4080 
4080  -1 
4080  1 
4080  2 
4080 
4080  -2 
4080  3 
3504 
3504  -2 
3504  2 
3504 
3504  -• 
3504 
3504  3 
2952 
2952  - 
2544 
1992 
1536 
1152  -1 
1152 
1152  1 
672 
672 
288 
288 
3504 
4059 
3504  -1 
3504  1 


Coarse  Mesh 
1  RP 


Medium  Mesh 


-0.2872 

0.2836 

0.3577 

0.0362 

0.2411 

0.2160 

-0.1806 

0.2952 

0.3067 

-0.1273 

0.2280 

0.2315 

-0.0789 

0.2530 

0.2348 

0.0520 

0.2609 

0.2358 

-0.0464 

0.2229 

0.2017 

-0.0813 

0.2794 

0.2579 

-0.0886 

0.2329 

0.2208 

-0.2085 

0.2582 

0.2941 

-0.1123 

0.2027 

0.2054 

-0.1675 

0.2473 

0.2647 

-0.0563 

0.1975 

0.1820 

-0.1146 

0.1961 

0.2013 

-0.0020 

0.1886 

0.1672 

-0.1480 

0.2335 

0.2450 

-0.0989 

0.2522 

0.2401 

0.0498 

0.2130 

0.1939 

-0.0935 

0.2085 

0.2025 

-0.1092 

0.2318 

0.2271 

-0.1177 

0.1789 

0.1898 

-0.1162 

0.1713 

0.1835 

-0.1666 

0.2039 

0.2333 

-0.1724 

0.2242 

0.2507 

-0.1588 

0.2319 

0.2491 

-0.1352 

0.1858 

0.2036 

-0.2000 

0.2134 

0.2592 

-0.1804 

0.2132 

0.2475 

0.0195 

0.1918 

0.1709 

0.1052 

0.2052 

0.2052 

-0.0059 

0.1888 

0.1674 

0.0464 

0.1962 

0.1787 

-0.2507 

0.2925 

0.3414 

0.0468 

0.2378 

0.2148 

-0.1647 

0.2923 

0.2973 

-0.1145 

0.2167 

0.2172 

-0.0593 

0.2520 

0.2294 

0.0667 

0.2571 

0.2354 

-0.0284 

0.2181 

0.1949 

-0.0710 

0.2596 

0.2385 

-0.0766 

0.2320 

0.2165 

-0.1973 

0.2589 

0.2885 

-0.0983 

0.2016 

0.1987 

-0.1666 

0.2435 

0.2615 

-0.0563 

0.1975 

0.1820 

-0.1059 

0.1938 

0.1957 

0.0022 

0.1834 

0.1625 

-0.1473 

0.2302 

0.2422 

-0.0966 

0.2484 

0.2362 

0.0525 

0.2103 

0.1921 

-0.0982 

0.2035 

0.2003 

-0.1018 

0.2306 

0.2234 

-0.1159 

0.1653 

0.1789 

-0.1073 

0.1648 

0.1743 

-0.1634 

0.2041 

0.2317 

-0.1701 

0.2226 

0.2483 

-0.1647 

0.2166 

0.2411 

-0.1435 

0.1711 

0.1979 

-0.1697 

0.2165 

0.2438 

-0.1507 

0.2144 

0.2323 

0.0211 

0.1865 

0.1664 

0.1093 

0.1949 

0.1981 

-0.0024 

0.1831 

0.1623 

0.0460 

0.1920 

0.1750 

Dense  Mesh 
I  RP 


.2842 

.2303 

.2883 

.2144 

.2561 

.2547 

.2149 

.2420 

.2420 

.2549 

.1935 

.2398 

.1970 

.1882 

.1760 

.2264 

.2501 

.2003 

.2085 

.2312 

.1635 

.1537 

.2103 

.2233 

.2266 

.1830 

.2178 

.2150 

.1789 

.1846 

.1747 

.1877 

Sum(E(X)) 

-2.94520 

7.13100 

7.22510 

-2.67660 

6.99170 

7.01860 

. 1 

-1.82320  6.91190 

6.76420 

Sum(E(XA2)) 

0.52795 

1.61991 

1.68734 

0.46380 

1.56144 

1.59066 

0.34947  1.52742 

1.47020 

Mean 

-0.09204 

0.22284 

0.22578 

-0.08364 

0.21849 

0.21933 

-0.05698  0.21600 

0.21138 

Variation 

0.00829 

0.00099 

0.00181 

0.00774 

0.00109 

0.00165 

I 

0.00792  0.00111 

0.00130 

350msec  comparison  due  to  long  time  sensor  data  inaccuracy. 
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Of  particular  interest  in  Table  17  are  the  errors  associated  with  the 
accelerometer  data  (indicated  by  a  sensor  number  preceded  by  an  “A”).  As  a  whole,  the 
errors  associated  with  these  sensors  are  significantly  lower  than  the  velocity  meters.  This 
gives  evidence  that  perhaps  a  large  portion  of  the  error  measurement  is  actually  error  in 
the  measured  data. 

Figure  35  is  a  plot  of  the  cumulative  Russell’s  error  factor  trend  for  each 
fluid  model  as  a  function  of  the  longitudinal  position.  . 


Shot  3  Russell's  Error  vs.  Longitudinal  Position 


Figure  35.  Russell’s  Comprehensive  Error  Trend  as  a  Function  of  Longitudinal 

Position 


For  Shot  3,  the  shock  wave  initially  impacts  the  structure  at  midships. 
Since  the  shock  wave  initiates  at  a  node  one  element  away  from  this  point  of  impact,  it 
propagates  through  very  little  of  the  fluid  model  before  reaching  the  structure.  The 
remainder  of  the  structure,  particularly  the  stern  and  bow  sections,  are  not  impacted  until 
the  shock  wave  has  propagated  through  more  of  the  fluid  and  structural  model.  The 
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resulting  error  is  reduced  by  increasing  the  normal  direction  mesh  density,  but  not 
eliminated  as  shown  in  Figure  38.  In  terms  of  fluid  mesh  density  comparison,  it  is 
apparent  that  the  mesh  density  has  the  largest  effect  in  these  outer  regions.  Very  little 
noticeable  improvement  is  shown  in  the  regions  closest  to  the  charge  location. 

b.  Velocity  Plots 

It  is  useful  to  observe  the  velocity  plots  to  help  visualize  how  the  predicted 
response  varies  as  the  fluid  mesh  density  is  increased.  Figure  36  shows  the  predicted  and 
measured  500  msec  velocity  response  plots  for  sensor  V2000V.  This  sensor  is  the 
forward  most  sensor,  located  above  the  sonar  dome.  This  sensor  resulted  in  the  greatest 
comprehensive  error  in  all  three  cases.  From  this  plot  it  is  evident  that  at  least  some  of 
this  error  is  due  to  the  late  time  drift  in  the  sensor.  This  drift  was  significant  and  could 
not  be  removed. 
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Grid  120217-vz  (V2000V) 


KLeel  (x=5328  y=0  z=82 


Time  (msec) 


Ship  Shock  Trial  -  Simulation  (COARSE) 

Simulation  (MEDIUM^  - Simulation  (DENSE) 


Figure  36.  Keel  Sensor  V2000V. 


Figure  37  shows  the  velocity  response  plots  for  keel  sensor  V2016V.  This 
sensor  is  significant  as  it  is  the  closest  sensor  to  the  charge.  This  plot  shows  much  better 
agreement  between  the  predicted  and  measured  response,  having  a  comprehensive  error 
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factor  of  around  0.20  for  all  three  fluid  models.  This  plot  also  shows  how  the  late  time 
predicted  response  settles  out  similarly  to  the  measured  response.  This  is  a  measure  of 
the  quality  of  the  structural  damping  coefficients  which  appear  to  be  very  accurate  in  this 
case. 


SHOT  3 

Grid  242 3 99- vz  (V2016V) 


KLeel  (x=2544  y=0  z=l  16) 


Ship  Shock:  Trial  - Simulation  (COARSE) 

Simulation  (MEDIUM)  - Simulation  (DENSE) 


Figure  37.  Keel  Sensor  V2016V 

The  third  sensor  shown  is  keel  sensor  V2020V  shown  in  Figure  38.  This 
sensor  is  one  of  the  furthest  aft  sensors,  located  at  frame  442.  Due  to  excessive  drift  in 
the  last  portion  of  the  sensor  data  only  350  msec  of  data  was  compared.  Though  better 
than  sensor  V2000V,  this  sensor  shows  how  there  are  apparent  transmission  losses  in  the 
shock  wave  as  it  propagates  outward  through  the  fluid  and  structural  model.  Here  the 
predicted  late  time  response  is  in  better  agreement  than  the  early  time  response.  As  in 
Figure  37,  this  is  evidence  of  accurate  structural  damping. 
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SHOT  3 

Grid  3  5 0052- vz  (V2020V) 


Keel  (x=288  y=0  z=21  1) 


Ship  Shock:  Trial  -  Simulation  (COARSE) 

Simulation  (MEDIUM)  -  Simulation  (DENSE) 


Figure  38.  Keel  Sensor  V2020V 

Figure  39  shows  bulkhead  sensor  V2124V.  This  sensor  is  bulkhead 
mounted  on  the  starboard  side,  just  forward  of  midships.  Located  higher  in  the  structure 
than  the  four  sensors  shown  above,  this  plot  shows  how  the  structure  dampens  the  shock 
wave  as  it  passes  through.  The  response  is  more  linear  and  settles  down  more  rapidly 
than  the  keel  response.  This  plot  shows  good  agreement,  both  in  phase  and  magnitude, 
between  the  measured  and  predicted  velocity  responses.  The  cumulative  Russell’s  error 
factor  for  this  comparison  was  0.1838  for  the  dense  fluid  mesh,  0.1957  for  the  medium 
density  fluid  mesh,  and  0.2013  for  the  coarse  fluid  mesh. 

The  velocity  plot  shown  in  Figure  40  uses  integrated  acceleration  data 
from  sensor  A2240V  located  on  the  starboard  mast  leg.  This  sensor  had  the  best 
correlation  between  measured  and  predicted  response  of  all  sensors  compared  for  this 
shot.  The  comprehensive  error  factor  of  the  comparison  using  the  dense  mesh  was 
0.1552.  Most  of  the  comparisons  utilizing  data  gathered  using  accelerometers  showed 
results  similar  to  this.  Several  additional  Shot  3  velocity  plots  are  included  in  Appendix 
A. 
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Vertical  Velocity  (ft/sec)  Vertical  Velocity  (ft/sec) 
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Figure  39.  Bulkhead  Sensor  V2124V 
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Figure  40. 


Mast  Sensor  A2240V 
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c. 


SHOT  1 


The  first  shot,  though  the  least  severe,  presented  an  interesting  comparison  due  to 
its  geometry.  This  shot  was  detonated  outside  of  its  intended  box  off  the  ships  port  bow, 
its  longitudinal  position  actually  forward  of  the  ship.  Figure  41  shows  the  relative  aspect 
of  this  shot  to  the  ship. 


Figure  41.  Shot  1  Geometry.  Most  Extreme  Aspect 

1.  Fluid  Mesh  Depth 

The  first  Shot  1  simulation  was  conducted  and  compared  using  the  1st  generation 
fluid  model  shown  in  Figure  22.  This  model  was  constructed  having  a  maximum  depth 
of  70  feet.  The  actual  cavitation  depth  for  Shot  1  at  the  point  closest  to  the  charge 
location  was  found  to  be  75  feet. 

The  average  RusselFs  comprehensive  error  factor  was  found  to  be  0.27179  using 
the  1st  generation  70-foot  depth  fluid  model.  The  2nd  generation  Shot  1  fluid  model  was 
constructed  with  a  maximum  depth  of  75  feet.  A  comparison  of  this  model  showed  that 
the  average  Russell’s  comprehensive  error  factor  was  reduced  to  0.21627.  The  majority 
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of  this  comprehensive  improvement  was  due  to  the  magnitude  error  which  was  reduced 
from  0.14389  down  to  -0.01656.  The  average  phase  error  was  reduced  from  0.25849 
down  to  0.21886,  a  significant  improvement  as  well.  Further  study  showed  that  the 
majority  of  this  improvement  was  due  primarily  due  to  the  change  in  depth.  The 
modification  of  the  shape  of  the  outer  boundary  from  having  flat  ends  to  having  spherical 
ends  also  contributed  to  this  improvement.  Table  18  lists  the  error  factors  for  all  of  the 
sensors  compared  using  both  models. 


Table  18.  Shot  1  Russell  Error  Factors  With  Varying  Fluid  Model  Depth 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Shot  1  -  840  inch  | 

j  Shot  1-  900  inch 

X 

y 

Z 

RM 

RP 

RC 

■  RM 

RP 

RC 

V2000V 

120217 

Keel 

0 

ESI 

mmm 

0.1803 

0.2022 

M 

0.1744 

0.1597 

V2002V 

142489 

Keel 

27 

m 

mm 

0.4053 

0.3648 

Bstesgii 

0.2899 

0.2711 

A8516V 

210993 

Bulkhead 

EH 

0.2414 

0.2964 

0.0062 

0.2055 

0.1822 

V2008VI 

210894 

Bulkhead 

WBsm 

0.3099 

0.3517 

0.0578 

0.2692 

0.2440 

V2007V 

210430 

Keel 

« 

0.3368 

0.3735 

0.0715 

0.3498 

0.3164 

V2009VI 

210808 

Bulkhead 

0.3533 

0.3966 

0.1959 

0.3595 

0.3628 

A2109V 

414953 

Bulkhead 

■a 

EB 

0.2083 

0.1891 

-0.0555 

0.0958 

0.0981 

V2010V 

220589 

Keel 

EBI 

0.1993 

0.1910 

-0.1086 

0.1513 

0.1650 

A2104V 

222240 

Bulkhead 

RBI 

0.2173 

0.2938 

0.1422 

0.1718 

0.1976 

V2011VI 

221188 

Bulkhead 

EH 

0.2528 

0.2611 

-0.1129 

0.2062 

0.2083 

V2013V 

221601 

Bulkhead 

En 

0.3119 

0.3865 

0.1532 

0.3042 

0.3018 

V2012VI 

221102 

Bulkhead 

EH 

0.3712 

0.3696 

-0.0186 

0.2786 

0.2475 

V2124V 

222060 

Bulkhead 

mm 

0.3750 

0.3723 

0.0986 

0.3887 

0.3554 

A2237V 

416419 

Bulkhead 

0.1794 

0.1892 

0.2310 

0.0676 

0.1256 

0.1265 

A2116V 

414367 

Bulkhead 

0 

0.0598 

0.1807 

0.1687 

-0.0462 

0.0922 

0.0914 

A2240V 

416269 

Bulkhead 

3504 

0.0176 

0.2288 

0.2034 

-0.0801 

0.1231 

0.1302 

V2014V 

230461 

Keel 

2952 

0.1698 

0.1772 

0.2175 

0.0253 

0.1730 

0.1550 

A3565V 

231696 

Bulkhead 

2952 

-0.1198 

0.2494 

0.2452 

-0.2076 

0.1380 

0.2209 

V2016V 

242399 

Keel 

2544 

0 

0.2568 

0.3535 

0.3872 

0.1142 

0.2836 

0.2710 

V2026V 

312302 

Keel 

1992 

0 

0.0195 

0.2213 

0.1969 

-0.1100 

0.1846 

0.1904 

A2310V 

320746 

Bulkhead 

1536 

0 

0.1251 

0.2677 

0.2619 

-0.0808 

0.2309 

0.2168 

V2035V 

330769 

Keel 

1152 

0.1703 

0.1771 

0.2177 

0.0141 

0.1634 

0.1454 

V2032V 

330764 

Keel 

1152 

0 

196 

0.1001 

0.2548 

0.2426 

-0.0847 

0.2287 

0.2162 

V2034V 

330759 

Keel 

1152 

-135 

0.0704 

0.2084 

0.1950 

-0.0418 

0.2625 

0.2356 

V2019V 

340167 

Keel 

672 

0 

0.0671 

0.2264 

0.2093 

-0.1403 

0.1901 

0.2094 

V2018VI 

340992 

Bulkhead 

672 

0 

0.1931 

0.2554 

0.2837 

0.0838 

0.2426 

0.2275 

V2020V 

350052 

Keel 

288 

0.1866 

0.2658 

0.2878 

-0.2590 

0.2646 

0.3282 

A2413V 

350220 

Bulkhead 

288 

0.1000 

0.2192 

0.2136 

-0.0959 

0.1804 

0.1811 

Sum(E(X)) 

4.02880 

7.23770 

7.61010 

-0.46380 

6.12820 

6.05550 

Sum(E(XA2)) 

0.83054 

1.99418 

2.21855 

0.34044 

1.50779 

1.45173 

Mean 

0.14389 

0.25849 

0.27179 

-0.01656 

0.21886 

0.21627 

Variation 

0.00929 

0.00457 

0.00556 

0.01232 

0.00617 

0.00526 

2.  Mesh  Quality 

Comparing  Shot  1  predictions,  the  3rd  generation  model  showed  minor 
improvement  over  the  1st  and  2nd  generation  model  overall,  reducing  the  average 
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comprehensive  Russell’s  error  factor  by  approximately  0.005.  Similar  to  what  was  found 
in  the  Shot  3  mesh  quality  comparison  discussed  previously,  individual  sensor 
comparisons  showed  mixed  results.  Despite  having  a  slightly  more  significant 
improvement  overall  than  Shot  3,  this  mesh  quality  comparison  proved  to  be  inconclusive 
due  to  these  inconsistencies. 

3.  Mesh  Density 

Shot  1  was  the  furthest  shot  from  the  ship.  As  a  result  the  fluid  models  were  the 
largest  cavitation  depth  models  constructed,  extending  to  a  depth  of  75  feet.  For  these 
models  the  number  of  normal  fluid  elements  in  the  outer  mesh  ranged  from  13  in  the 
coarse  mesh  to  38  in  the  dense  mesh.  Table  19  summarizes  the  mesh  properties  of  three 
of  the  fluid  models  constructed  for  the  Shot  1 . 


Table  19.  Shot  1  Fluid  Model  Mesh  Properties 


Property 

Coarse  Mesh 

Medium  Mesh 

Dense  Mesh 

Inner  Mesh  Elements 

5 

6 

7 

Inner  Mesh  Element 
Spacing  (Ave.) 

11.4  in 

10.0  in 

9.3  in 

Transition  Mesh 
Elements 

5 

6 

7 

Transition  Mesh 
Element  Spacing  (Ave.) 

13.4  in 

1 1.6  in 

10.1  in 

Outer  Mesh  Elements 

14 

24 

38 

Outer  Mesh  Element 
Spacing  (Ave.) 

36  in 

20  in 

12  in 

Fluid  Nodes 

176,287 

244,049 

328,541 

DAA  Elements 

5048 

5090 

5132 

Simulation  Run  Time 

34.75  hrs 

40.48  hrs 

47.33  hrs 

The  simulations  of  Shot  1  took  considerably  longer  to  run  than  Shot  3  due  to  the 
larger  volume  of  the  fluid  mesh  requiring  more  nodes.  While  all  of  the  Shot  3  runs  could 
be  completed  in  a  day  and  a  half  timeframe,  the  medium  and  dense  mesh  Shot  1  runs  took 
nearly  two  days  to  process.  While  less  than  5  hours  separated  the  run  times  for  all  of  the 
Shot  3  fluid  meshes,  there  is  over  a  12  hour  time  difference  between  the  coarse  Shot  1 
mesh  and  the  dense  Shot  1  mesh  run  times.  This  time  difference  is  indicative  of  how 
quickly  the  model  grows  in  size  when  extending  the  depth  of  the  fluid  model. 
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a.  Error  Comparison 

A  summary  of  the  Russell’s  error  factors  found  for  each  model 
comparison  are  shown  in  Figures  42  through  44. 


Russell's  Error  Factor  (Shot  1  Coarse) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  •  Excellent 


Figure  42.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison 


Russell's  Error  Factor  (Shot  1  Medium) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  •  Excellent 


Figure  43.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison 
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Russell's  Error  Factor  (Shot  1  Dense) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  •  Excellent 

Figure  44.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison 

Unlike  the  Shot  3  comparison  error  plots,  there  are  no  apparent  trends  in 
the  Shot  1  comparisons.  The  plots  suggest  that  there  is  a  large  variation  in  both  phase 
and  magnitude  error  as  the  error  factors  are  very  scattered  in  all  three  plots.  Interestingly, 
the  coarse  model  had  smaller  errors  than  the  other  two  models,  opposite  of  what  was 
observed  in  the  Shot  3  case. 

Table  20  lists  all  of  the  sensors  compared  and  their  corresponding 
Russell’s  error  factors.  The  averages  for  each  error  factor  are  shown  at  the  bottom  of 
their  respective  columns  along  with  the  variation.  These  averages  show  that,  the  coarse 
mesh  does  in  fact  have  lower  average  comprehensive  and  phase  error  factor  at  0.21019 
and  0.20767  respectively,  while  the  dense  mesh  has  the  lowest  average  magnitude  error 
factor  at  0.00094.  It  should  be  noted  that  the  average  magnitude  error  factor  can  portray 
a  better  correlation  than  what  really  exists  since  negative  magnitude  error  factors  exist. 
The  average  magnitude  error  then  indicates  how  much  less  or  more  the  predicted 
response  magnitudes  were  overall.  The  variation  is  a  good  indicator  of  how  consistent 
the  predictions  were.  In  this  case,  the  variation  of  the  dense  mesh  comparison  is  largest 
for  three  error  factors  indicating  that  not  only  was  the  dense  mesh  the  least  accurate 
overall,  but  it  was  also  the  least  consistent. 
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Table  20.  Shot  1:  Russell’s  Error  Factor 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Mesh  j 

|  Medium  Mesh  j 

j  Coarse  Mesh 

X 

y 

Z 

RM 

RP 

RC 

RM 

RP 

RC 

RM 

RP 

RC 

V2026V 

312302 

Keel 

lEEH 

HI 

m 

0.1275 

0.2417 

0.2422 

0.0811 

0.2575 

0.2392 

0.0534 

0.2684 

0.2425 

V2000V 

120217 

Keel 

5328 

m 

-0.2684 

0.2621 

0.3325 

-0.2713 

0.2670 

0.3374 

-0.2670 

0.2654 

0.3336 

V2002V 

142489 

Keel 

grass 

m 

0.0797 

0.1723 

0.1682 

0.0729 

0.1711 

0.1648 

0.0772 

0.1755 

0.1699 

V2007V 

210430 

Keel 

m 

-0.0706 

0.1914 

0.1808 

-0.1041 

0.1833 

0.1868 

-0.1178 

0.1760 

0.1877 

V2010V 

220589 

Keel 

m 

0.0073 

0.2670 

0.2367 

0.0420 

0.2710 

0.2431 

-0.0620 

0.2697 

0.2452 

V2014V 

230461 

Keel 

EI 

0.1370 

0.2433 

0.2474 

0.0911 

0.2395 

0.2271 

0.0653 

0.2440 

0.2238 

V2016V 

242399 

Keel 

nisi 

-0.0723 

0.2032 

0.1912 

-0.1185 

0.2009 

0.2067 

-0.1396 

0.2013 

0.2171 

A2310V 

320746 

Bulkhead 

1536 

0.1492 

0.2711 

0.2742 

0.0905 

0.2525 

0.2377 

0.0584 

0.2458 

0.2239 

V2008VI 

210894 

Bulkhead 

gran 

-0.0580 

0.2289 

0.2092 

-0.0772 

0.2264 

0.2120 

-0.0887 

0.2284 

0.2171 

V2009VI 

210808 

Bulkhead 

0.0482 

0.1607 

0.1486 

0.0048 

0.1541 

0.1367 

-0.0121 

0.1480 

0.1316 

V2011VI 

221188 

Bulkhead 

-0.1420 

0.1996 

0.2171 

-0.1822 

0.1834 

0.2291 

-0.1863 

0.1772 

0.2279 

V2012VI 

221102 

Bulkhead 

0.0584 

0.1745 

0.1631 

0.0095 

0.1784 

0.1583 

-0.0116 

0.1778 

0.1579 

V2034V 

330759 

Keel 

0.1365 

0.2649 

0.2641 

0.0699 

0.2393 

0.2209 

0.0296 

0.2233 

0.1996 

V2035V 

330769 

Keel 

0.1192 

0.3301 

0.3110 

0.1720 

0.2726 

0.2856 

0.1321 

0.2592 

0.2578 

V2032V 

330764 

Keel 

1152 

0.1891 

0.2754 

0.2960 

0.1407 

0.2611 

0.2628 

1022 

0.2501 

0.2394 

V2019V 

340167 

Keel 

672 

-0.1004 

0.2598 

0.2469 

-0.1539 

0.2463 

0.2574 

-0.1463 

0.2447 

0.2526 

A8516V 

210993 

Bulkhead 

-0.0497 

0.2324 

0.2107 

-0.0676 

0.2363 

0.2173 

-0.0751 

0.2389 

0.2220 

V2020V 

350052 

Keel 

-0.0370 

0.3792 

0.3376 

-0.0909 

0.3782 

0.3447 

-0.0827 

0.4013 

0.3631 

A2413V 

350220 

Bulkhead 

0.0781 

0.3461 

0.3144 

0.0200 

0.3449 

0.3062 

0.0171 

0.3649 

0.3238 

V2013V 

221601 

Bulkhead 

-0.0933 

0.1866 

0.1849 

-0.1288 

0.1827 

0.1981 

-0.1366 

0.1819 

0.2016 

A3565V 

231696 

Bulkhead 

0.1805 

0.1773 

0.2242 

0.1324 

0.1623 

0.1856 

1093 

0.1557 

0.1686 

V2018VI 

340992 

Bulkhead 

672 

-0.0775 

0.2370 

0.2210 

-0.1242 

0.2230 

0.2262 

-0.1149 

0.2197 

0.2198 

A2104V 

222240 

Bulkhead 

-0.1963 

0.1496 

0.2187 

-0.2204 

0.1297 

0.2266 

-0.2235 

0.1184 

0.2241 

V2124V 

222060 

Bulkhead 

-0.0741 

0.1653 

0.1605 

-0.1200 

0.1509 

0.1708 

-0.1330 

0.1473 

0.1759 

A2116V 

414367 

Bulkhead 

-0.0402 

0.1149 

0.1079 

-0.0679 

0.0915 

0.1010 

-0.0679 

0.0848 

0.0962 

A2109V 

414953 

Bulkhead 

0.0904 

0.1239 

0.1359 

0.0709 

0.1281 

0.1298 

0.0682 

0.1342 

0.1334 

A2237V 

416419 

Bulkhead 

-0.0634 

0.1386 

0.1351 

-0.0915 

0.1205 

0.1341 

-0.0906 

0.1203 

0.1334 

A2240V 

416269 

Bulkhead 

-0.0316 

0.1122 

0.1033 

-0.0567 

0.0943 

0.0975 

-0.0556 

0.0926 

0.0957 

Sum(EfX)) 

0.02630 

6.10910 

6.08340 

-0.87740 

5.84680 

5.94350 

-1.29850 

5.81480 

5.88520 

Sum(E(XA2)) 

0.36994 

1.45976 

1.43693 

0.39283 

1.35133 

1.36958 

1  0.36280 

1.35485 

1.34893 

Mean 

0.00094 

0.21818 

0.21726 

-0.03134 

0.20881 

0.21227 

I  -0.04638 

0.20767 

0.21019 

Variation 

0.01370 

0.00470 

0.00427 

0.01353 

0.00483 

0.00400 

1  0.01121 

0.00545 

0.00415 
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Figure  45  shows  the  Russell’s  comprehensive  error  factor  trend  in  relation 
to  the  longitudinal  position  on  the  ship  for  each  fluid  mesh.  As  is  shown,  all  three 
meshes  have  poor  correlation  in  the  sonar  dome  region  as  was  the  case  for  Shot  3.  For 
Shot  1  however,  the  sonar  dome  region  was  the  closest  region  on  the  ship  to  the  charge. 
This  should  have  resulted  in  much  better  accuracy  in  the  prediction.  This  anomaly  may 
be  due  to  other  factors  such  as  gauge  error  or  the  effects  of  the  sonar  dome  which  has 
been  omitted  from  the  structural  model. 


Shot  1  Russell's  Error  vs.  Longitudinal  Position 


Figure  45.  Shot  1  Russell’s  Comprehensive  Error  Trends 


Very  good  correlation  exists  in  the  forward  areas  of  the  ship,  particularly 
for  those  sensors  located  in  the  superstructure  of  the  ship.  The  trends  begin  to  diverge  aft 
of  midships  where  the  coarse  mesh  correlation  is  significantly  better  than  the  medium  and 
dense  mesh.  This  is  contradicted  by  the  stem  region  correlation  which  is  poor  for  all 
meshes,  but  poorest  for  the  coarse  mesh. 
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The  interaction  at  the  hull  is  a  result  of  the  pressure  wave  as  well  as  the 
fluid  particle  velocities.  While  pressure  always  acts  nonnal  to  the  structure,  the  fluid 
particle  velocities  do  not.  In  Shot  3,  the  fluid  particle  velocities  acted  primarily  nonnal  to 
the  majority  of  the  structure  since  the  charge  was  located  off  the  beam.  In  Shot  1,  the 
charge  was  located  very  far  forward  off  the  bow.  A  large  portion  of  the  structure, 
primarily  the  aft  portion  of  the  ship,  is  subjected  to  fluid  particle  velocities  which  are 
nearly  tangent  to  the  hull.  This,  coupled  with  the  fact  that  the  incident  pressure  wave 
decreases  as  the  wave  propagates  aft,  tends  to  magnify  any  prediction  error. 

The  increased  error  in  the  dense  fluid  mesh  in  the  aft  sections  of  the  ship 
may  also  be  due  to  this  nearly  tangential  flow.  In  a  shot  off  the  beam,  the  element  aspect 
ratios  of  concern  are  those  which  make  up  the  DAA  boundary  and  the  wetted  surface 
structure  interface.  When  the  shot  is  a  bow  or  stern  shot,  the  cross-sectional  aspect  ratios 
of  the  fluid  element  become  important  in  calculating  the  wave  propagation.  While  Shot  1 
is  not  a  direct  end-on  shot,  as  the  wave  propagates  aft,  the  propagation  becomes 
increasingly  tangential  to  the  hull.  As  the  density  of  the  mesh  is  increased,  the  aspect 
ratios  become  much  smaller  and  the  resulting  approximation  less  accurate  [Ref.  27]. 

b.  Velocity  Plots 

To  help  visualize  the  Russell’s  error  factor  correlation  discussed,  five 
velocity  plots  are  shown  in  Figures  46  through  50.  Additional  velocity  plots  can  be  found 
in  Appendix  A.  The  first  plot  is  sensor  V2000V,  located  in  the  sonar  dome  region.  This 
plot  shows  little  correlation  between  the  predicted  and  measured  response.  The  measured 
response  is  oscillating  at  a  higher  frequency  and  magnitude  than  the  predicted  response. 
This  poor  correlation  is  similar  to  that  of  Shot  3  indicative  of  the  possible  contributing 
factors  mentioned  earlier. 

Figure  47  shows  bulkhead  sensor  V2124V,  located  forward  of  midships 
and  above  the  waterline.  This  shows  exceptional  correlation  in  all  three  cases.  The  peak 
response  is  captured  very  well  as  is  the  recovery  region  immediately  following  this  peak. 
The  predicted  response  settles  out  faster  than  the  measured  response,  suggesting  that  the 
model  may  be  over-damped  in  this  case. 
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Vertical  Velocity  (ft/s)  Vertical  Velocity  (ft/s) 


SHOT  1 


Grid  120217-vz  (V2000V) 


ICeel  (x=5328  y=0  z=82) 


Time  (msec) 

Ship  Shock.  Trial  -  Simulation  (COARSE) 

Simulation  (MEDIUM)  -  Simulation  (DENSE) 


Figure  46.  Keel  Sensor  V2000V. 

SHOT  1 

Grid  222060-vz  (V2124V) 


Bulkhead  (x=3504  y=-375  z=390) 


Time  (msec) 

Ship  Shock  Trial  -  Simulation  (COARSE) 

Simulation  (MEDIUM)  -  Simulation  (DENSE) 


Figure  47.  Bulkhead  Sensor  V2124V. 
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Figure  48  is  keel  sensor  V2016V  located  just  aft  of  midships.  This  shows 
good  correlation  in  the  early  and  late  time  response  with  a  bit  of  a  phase  miscorrelation 
after  the  initial  peak.  This  is  most  apparent  in  the  dense  mesh  prediction.  This  region  of 
the  ship  is  where  the  predictions  from  each  mesh  began  to  diverge. 


SHOT  1 

Grid  242 3 99- vz  (V2016V) 


Keel  (x=2544  y=0  z=l  16) 


Time  (  msec  ) 

Ship  Shock:  Trial  - Simulation  (COARSE) 

Simulation  ( Pvl EDIU Wl )  - Simulation  (DENSE) 


Figure  48.  Keel  Sensor  V2016V 


Moving  to  the  furthest  aft  location  compared,  Figure  49  shows  keel  sensor 
V2020V.  This  shows  poor  correlation  both  in  phase  and  in  magnitude.  The  comparisons 
for  all  of  the  meshes  were  poor  for  this  sensor.  As  discussed  earlier,  this  can  be  attributed 
to  the  tangential  propagation  of  the  shock  wave  in  this  region. 

Shown  in  Figure  50  is  sensor  A2240V,  an  accelerometer,  was  located  high 
in  the  structure  on  the  starboard  mast  leg.  Excellent  phase  and  magnitude  correlation  is 
shown  in  this  comparison.  This  excellent  correlation  is  evidence  of  proper  structural 
damping  in  the  model  and  a  minimal  amount  of  measured  data  error.  As  with  Shot  3, 
most  of  the  comparisons  made  with  accelerometer  data  showed  significantly  lower 
predictions  errors  than  comparisons  made  with  velocity  meter  data. 
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Vertical  Velocity  (ft/s)  Vertical  Velocity  (ft/s) 


SHOT  1 

Grid  350052  -vz  (V2020V) 


Reel  (x=288  y=0  z=21  1) 


Time  (msec) 

Ship  Shock:  Trial  -  Simulation  (COARSE) 

Simulation  ( IV1  ED  I LJ IVI)  -  Simulation  (DENSE) 


Figure  49.  Keel  Sensor  V2020V. 


SHOT  1 


Grid  416269-vz  (A2240V) 

Ivlast  (x=3504  y=-135  z=848) 


lOO 


200  300 

Time  (msec) 


Ship  Shock  Trial 
Simulation  (MEDIUM^ 


Simulation  (COARSE) 
Simulation  (DENSE) 


500 


Figure  50.  Mast  Sensor  A2240V. 
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D. 


SHOT  2 


The  final  shot  compared  is  Shot  2.  This  shot  was  the  only  shot  detonated  on  the 
starboard  side  of  DDG  81.  In  this  shot,  the  charge  was  located  off  the  beam,  slightly 
forward  of  where  the  Shot  3  charge  was  detonated  on  the  port  side  and  with  a  larger 
stand-off  distance.  Figure  5 1  shows  the  relative  aspect  of  the  Shot  2  charge  to  the  ship. 


Figure  51.  Shot  2  Geometry.  Charge  Off  Starboard  Midships 


Three  3ld  generation  models  of  varying  density  were  constructed  for  the  Shot  2 
simulation.  A  summary  of  the  fluid  model  properties  is  shown  in  Table  21. 
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Phase  Erroi 


Table  21.  Shot  2  Fluid  Model  Properties 


Property 

Coarse  Mesh 

Medium  Mesh 

Dense  Mesh 

Inner  Mesh  Elements 

5 

6 

7 

Inner  Mesh  Element 
Spacing  (Ave.) 

11.4  in 

10.0  in 

9.3  in 

Transition  Mesh 
Elements 

5 

6 

7 

Transition  Mesh 
Element  Spacing  (Ave.) 

13.4  in 

1 1.6  in 

10.1  in 

Outer  Mesh  Elements 

10 

20 

32 

Outer  Mesh  Element 
Spacing  (Ave.) 

36  in 

20  in 

12  in 

Fluid  Nodes 

150,879 

213,329 

287,197 

DAA  Elements 

5048 

5090 

5132 

Simulation  Run  Time 

31.5  hrs 

37.2  hrs 

42.8  hrs 

1.  Error  Comparison 

A  summary  of  the  Russell’s  error  factors  found  in  the  comparison  of  all  three 
meshes  is  shown  in  Figures  52  through  54. 


Russell's  Error  Factor  (Shot  2  Coarse) 


Figure  52.  Russell’s  Error  Factors  for  Coarse  Mesh  Fluid  Model  Comparison 
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Russell's  Error  Factor  (Shot  2  Medium) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  •  Excellent 


Figure  53.  Russell’s  Error  Factors  for  Medium  Mesh  Fluid  Model  Comparison 


Russell's  Error  Factor  (Shot  2  Dense) 


-0.4000  -0.3000  -0.2000  -0.1000  0.0000  0.1000  0.2000  0.3000  0.4000 

Magnitude  Error 

•  Poor  •  Acceptable  •  Excellent 


Figure  54.  Russell’s  Error  Factors  for  Dense  Mesh  Fluid  Model  Comparison 
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The  trends  in  these  comparisons  are  very  similar  to  those  observed  in  Shot  3.  The 
dense  mesh  shows  significant  improvement  in  the  overall  error  factors.  This 
improvement  is  most  significant  in  terms  of  the  magnitude  error.  As  in  Shot  3,  the 
general  trend  of  the  error  factor  plots  is  that  the  points  move  primarily  to  the  right, 
indicative  of  a  larger  velocity  magnitude  response  being  predicted  as  the  density  of  the 
fluid  mesh  is  increased.  The  phase  error  remains  fairly  constant  throughout. 

Table  22  lists  all  of  the  sensors  compared  and  their  corresponding  Russell’s  error 
factors.  The  averages  for  each  error  factor  are  shown  at  the  bottom  of  their  respective 
columns  along  with  the  variation.  As  shown  in  the  table,  the  dense  mesh  has  a 
significantly  lower  magnitude  error  than  the  other  two  meshes.  As  was  the  case  with 
Shot  3,  the  most  significant  improvement  is  shown  between  the  medium  and  dense 
meshes.  This  shot  comparison  represents  the  best  overall  comparison  of  the  three  shots 
compared.  Every  mesh  compared  for  this  shot  had  a  mean  comprehensive  error 
significantly  less  than  the  other  two  shot  comparisons. 

Figure  55  shows  the  comprehensive  error  trends  of  each  mesh  comparison  over 
the  length  of  the  ship.  This  trend  plot  shows  that  the  most  accurate  predictions  occur 
around  the  midship  regions,  while  the  least  accurate  predictions  are  at  either  end.  The 
dense  mesh  shows  a  significantly  better  prediction  at  these  ends  and  in  the  aft  section  of 
the  ship.  By  comparison,  the  Shot  3  trend  was  very  similar. 
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Table  22.  Shot  2:  Russell’s  Error  Factor 


Sensor 

NODE 

Bulkhead 

Location  (in) 

Shot  2:  Coarse  Mesh  | 

Shot  2 

:  Medium  Mesh  | 

[  Shot  2:  Dense  Mesh 

Keel 

X 

y 

Z 

RM 

RP 

RC 

RM 

RP 

RC 

RM 

RP 

RC 

V2000V 

120217 

Keel 

5328 

0 

82 

-0.3953 

0.2976 

0.4385 

-0.3671 

0.2936 

0.4165 

-0.3385 

0.2971 

0.3991 

V2002V 

142489 

Keel 

4656 

27 

82 

0.0236 

0.2101 

0.1874 

0.0309 

0.2032 

0.1822 

0.0579 

0.2001 

0.1846 

V2009VI 

210808 

Bulkhead 

4080 

-174 

177 

-0.1204 

0.2212 

0.2232 

-0.1120 

0.2307 

0.2272 

-0.0882 

0.2272 

0.2160 

V2008VI 

210894 

Bulkhead 

4080 

174 

177 

-0.1190 

0.2037 

0.2091 

-0.1136 

0.1968 

0.2014 

-0.0714 

0.1994 

0.1877 

V2007V 

210430 

Keel 

4080 

0 

82 

-0.0701 

0.2098 

0.1961 

-0.0623 

0.2095 

0.1937 

-0.0245 

0.2092 

0.1866 

V2108V 

212196 

Bulkhead 

4080 

312 

390 

-0.1227 

0.1883 

0.1992 

-0.1186 

0.1843 

0.1943 

-0.0788 

0.1872 

0.1800 

A8516V 

210993 

Bulkhead 

4080 

238.5 

202 

-0.0303 

0.1625 

0.1465 

-0.0217 

0.1571 

0.1405 

0.0311 

0.1623 

0.1464 

A2109V 

414953 

Bulkhead 

4059 

0 

722.8 

0.0851 

0.1507 

0.1534 

0.0830 

0.1549 

0.1557 

0.1170 

0.1647 

0.1790 

V2124V 

222060 

Bulkhead 

3504 

-375.4 

390 

-0.1929 

0.1981 

0.2451 

-0.1956 

0.2055 

0.2514 

-0.1793 

0.2113 

0.2456 

V2012VI 

221102 

Bulkhead 

3504 

-216 

177 

-0.1512 

0.2135 

0.2319 

-0.1507 

0.2172 

0.2343 

-0.1299 

0.2211 

0.2273 

V2011VI 

221188 

Bulkhead 

3504 

216 

177 

-0.0104 

0.2293 

0.2034 

-0.0015 

0.2259 

0.2002 

0.0412 

0.2240 

0.2018 

V2010V 

220589 

Keel 

3504 

0 

85.5 

-0.1194 

0.2077 

0.2123 

-0.1147 

0.2082 

0.2106 

-0.0827 

0.2070 

0.1975 

V2125V 

222436 

Bulkhead 

3504 

375.4 

390 

-0.0627 

0.194 

0.1807 

-0.0549 

0.1937 

0.1785 

-0.0214 

0.1914 

0.1707 

V2013V 

221601 

Bulkhead 

3504 

0 

280 

-0.1344 

0.1413 

0.1728 

-0.1326 

0.1403 

0.1710 

-0.1049 

0.1420 

0.1565 

A2240V 

416269 

Bulkhead 

3504 

-135 

848 

-0.0700 

0.1292 

0.1302 

-0.0779 

0.1311 

0.1352 

-0.0626 

0.1335 

0.1307 

A2104V 

222240 

Bulkhead 

3504 

0 

390 

-0.0623 

0.0965 

0.1018 

-0.0636 

0.0963 

0.1023 

-0.0388 

0.0953 

0.0912 

A2116V 

414367 

Bulkhead 

3504 

0 

702 

-0.0088 

0.0986 

0.0877 

-0.0134 

0.1011 

0.0903 

0.0070 

0.1003 

0.0891 

A2237V 

416419 

Bulkhead 

3504 

135 

848 

-0.0148 

0.0968 

0.0868 

-0.0159 

0.0991 

0.0889 

0.0047 

0.0931 

0.0826 

V2014V 

230461 

Keel 

2952 

0 

85.5 

-0.0696 

0.2213 

0.2056 

-0.0701 

0.2173 

0.2024 

-0.0590 

0.2126 

0.1956 

A3565V 

231696 

Bulkhead 

2952 

-81 

316.7 

0.0184 

0.1749 

0.1558 

0.0165 

0.1699 

0.1513 

0.0290 

0.1675 

0.1506 

V2016V 

242399 

Keel 

2544 

0 

116 

-0.0534 

0.2102 

0.1923 

-0.0428 

0.2082 

0.1884 

-0.0169 

0.2038 

0.1812 

V2026V 

312302 

Keel 

1992 

0 

55 

-0.0931 

0.2254 

0.2162 

-0.0922 

0.2224 

0.2134 

-0.0751 

0.2185 

0.2047 

A2310V 

320746 

Bulkhead 

1536 

0 

177 

-0.0155 

0.2194 

0.1949 

-0.0121 

0.2102 

0.1866 

0.0103 

0.2047 

0.1816 

V2034V 

330759 

Keel 

1152 

-135 

193.3 

-0.1679 

0.1899 

0.2247 

-0.1601 

0.1898 

0.2200 

-0.1442 

0.1874 

0.2095 

V2032V 

330764 

Keel 

1152 

0 

196 

-0.1232 

0.1697 

0.1859 

-0.1183 

0.1669 

0.1813 

-0.1067 

0.1593 

0.1699 

V2035V 

330769 

Keel 

1152 

135 

193.3 

-0.0106 

0.1894 

0.1681 

-0.0090 

0.1823 

0.1618 

0.0009 

0.1692 

0.1500 

V2019V 

340167 

Keel 

672 

0 

196.5 

-0.1747 

0.2229 

0.2510 

-0.1727 

0.2222 

0.2494 

-0.1327 

0.2391 

0.2423 

V2018VI 

340992 

Bulkhead 

672 

0 

364.8 

-0.1509 

0.2152 

0.2329 

-0.1492 

0.2176 

0.2338 

-0.1076 

0.2369 

0.2306 

V2020V 

350052 

Keel 

288 

0 

211 

-0.2916 

0.2676 

0.3507 

-0.2760 

0.2592 

0.3355 

-0.2477 

0.2657 

0.3219 

A2413V 

350220 

Bulkhead 

288 

0 

273.4 

-0.1822 

0.2169 

0.2510 

-0.1625 

0.2104 

0.2356 

-0.1245 

0.2219 

0.2255 

Sum(E(X)) 

-2.89030 

5.77170 

6.03520 

-2.75070 

5.72490 

5.93370 

-1.93630 

5.75280 

5.73580 

Sum(E(XA2)) 

0.55106 

1.17304 

1.35421 

0.50179 

1.15183 

1.29863 

0.37537 

1.16764 

1.21171 

Mean 

-0.09634 

0.19239 

0.20117 

-0.09169 

0.19083 

0.19779 

-0.06454 

0.19176 

0.19119 

Variation 

0.00940 

0.00216 

0.00483 

0.00861 

0.00205 

0.00431 

0.00863 

0.00222 

0.00397 
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Shot  2  Russell's  Error  vs.  Longitudinal  Position 


Figure  55.  Russell’s  Comprehensive  Error  Trend  as  a  Function  of  Longitudinal 

Position 


2.  Velocity  Plots 

It  is  useful  to  correlate  the  errors  summarized  above  using  actual  velocity 
comparison  plots.  The  same  sensor  velocity  plots  shown  for  Shotl  and  Shot  3  are  shown 
below  in  Figures  56  through  60.  Several  more  velocity  plots  from  Shot  2  are  included  in 
Appendix  A. 

Figure  56  shows  keel  sensor  V2000V  located  above  the  sonar  dome.  This  sensor 
has  shown  little  correlation  in  the  Shot  1  and  Shot  3  comparisons,  the  Shot  2  comparison 
is  no  exception.  Though  the  best  overall  correlation,  the  Shot  2  prediction  at  this  sensor 
showed  the  worst  correlation  of  all  the  sensors  compared  for  all  three  shots. 
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SHOT  2 


Grid  120217-vz  (V2000V) 

ICeel  (x=5328  y=0  z=82) 


Time  (msec) 


Ship  Shock:  Trial 
Simulation  ( 1VI EDI l_J  IVl) 


Simulation  (COARSE) 
Simulation  (DEN SE) 


Figure  56. 


Keel  Sensor  V2000V 


Figure  57  shows  bulkhead  sensor  V2124V.  This  sensor  is  located  higher  in  the 
structure,  making  the  prediction  more  dependent  on  the  structural  model  and  structural 
damping.  Though  showing  signs  of  being  slightly  over  damped,  this  comparison  shows 
fairly  good  phase  correlation.  There  is  drift  apparent  in  the  late  time  measure  response 
data  which  explains  one  reason  why  this  sensor  comparison  had  one  of  the  highest 
comprehensive  error  factors  of  the  Shot  2  sensors  compared. 

Further  aft  and  lower  on  the  ship  is  keel  sensor  V2016V  shown  in  Figure  58.  The 
velocity  comparison  at  this  sensor  shows  exceptional  peak  response  agreement,  followed 
by  a  relatively  good  correlating  late  time  response  magnitude,  but  with  a  small  phase  shift 
in  all  three  cases.  Minor  measured  data  drift  is  apparent  in  the  late  time  response, 
contributing  some  to  the  error  factor  calculated  for  this  sensor. 
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Vertical  Velocity  (ft/s)  Vertical  Velocity  (ft/s) 


SHOT  2 

Grid  222060-vz  (V2124V) 


Bulkhead  (x=3504  y=-375  z=39Q) 


Time  (msec) 

—  Ship  Shock:  Trial  - Simulation  (COARSE) 

—  Simulation  (MEDIUM)  -  Simulation  (DENSE) 

Figure  57.  Bulkhead  Sensor  V2124V 

SHOT  2 

Grid  242399-vz  (V2Q16V) 


Keel  (x=2544  y=0  z=l  16) 


Time  (msec) 

Ship  Shock  Trial  -  Simulation  (COARSE) 

Simulation  (MEDIUM)  -  Simulation  (DENSE) 

Figure  58.  Keel  Sensor  V2016V 


86 


The  furthest  aft  sensor  compared  is  keel  sensor  V2020V  shown  in  Figure  59. 
This  comparison  showed  the  second  worst  correlation  of  all  of  the  Shot  2  sensors 
compared.  In  fact,  for  Shot  2  only  two  sensor  predictions  did  not  meet  the  acceptance 
criteria  of  0.28  Russell’s  comprehensive  error  factor  or  lower.  Sensor  V2000V  shown  in 
Figure  56  and  this  sensor  were  those  sensors.  As  shown,  little  correlation  exists  both  in 
early  and  the  late  time  response.  The  predicted  initial  peak  magnitude  is  approximately 
half  that  which  was  measured  for  all  three  fluid  meshes.  The  late  time  predicted  response 
continues  to  underestimate  the  magnitude  and  there  is  significant  phase  error. 
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SHOT  2 

Grid  350052-vz  (V2020V) 
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Simulation  ( IVI EDI IJ  M) 
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Figure  59.  Keel  Sensor  V2020V 


The  last  sensor  comparison  shown,  mast  sensor  A2237V  shown  in  Figure  60,  had 
the  best  correlation  of  all  of  the  comparisons  for  all  three  shots.  The  early  time 
correlation  between  predicted  and  measured  response  was  nearly  perfect.  Exceptional 
correlation  was  shown  in  the  late  time  response  correlation.  Only  a  minor  phase  error 
exists  and  the  model  appears  to  be  slightly  over-damped.  The  comprehensive  error  factor 
for  this  sensor  comparison  was  below  0.9  in  all  three  cases,  with  the  majority  of  the  error 
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due  to  the  phase  error.  As  in  the  previous  two  shots  compared,  significantly  better 
correlation  was  observed  when  using  measured  data  from  accelerometers  than  with 
velocity  meter  data. 


SHOT  2 

Grid  416419-vz  (A2237V) 


Mast  (x=3504  y=  135  z=848) 


O  lOO  200  300  400  500 

Time  (msec) 


- Ship  Shock  Trial  - Simulation  (COARSE) 

- Simulation  (MEDIUM)  -  Simulation  (DENSE) 

Figure  60.  Mast  Sensor  A2237V 

E.  QUALITATIVE  COMPARISON 

In  conducting  the  fluid  mesh  density  comparisons,  several  trends  emerged  which 
were  prevalent  in  all  three  shot  comparisons.  In  all  cases,  the  predicted  response 
correlation  was  better  when  comparing  sensors  mounted  higher  up  in  the  structure. 
Figure  61  shows  the  error  trend  in  relation  to  the  vertical  position  on  the  ship.  This  plot 
was  developed  by  averaging  all  nine  shot  comparisons  at  each  vertical  region  of  the  ship. 
This  plot  shows  the  significant  improvement  in  correlation  in  the  higher  regions  of  the 
structure.  This  improvement  can  be  contributed  to  a  good  fluid  and  structural  model, 
accurate  proportional  structural  damping  coefficient,  and  less  measured  data  drift.  The 
decrease  in  correlation  in  the  lower  regions  may  be  due  to  the  fact  that  the  response  is 
generally  more  dramatic  in  this  region  amplifying  any  error  between  the  predicted  and 

measured  response.  An  additional  source  of  error  may  be  due  to  the  fact  that  the  model, 
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though  accounting  for  the  liquid  loading  mass,  does  not  model  the  fluid-structure 
interface  in  the  water  compensated  fuel  oil  tanks  which  line  the  majority  of  the  lower 
hull.  To  do  this  would  require  an  extraordinary  amount  of  effort  due  to  the  complex 
geometries  inherent  to  the  inner  regions  of  these  tanks.  The  fluid  in  these  tanks  most 
certainly  has  an  effect  on  the  response  of  the  actual  ship. 


Vertical  Comprehensive  Error  Trend 


Russell's  Comprehensive  Error 


Figure  61.  Average  Russell’s  Comprehensive  Error  Factor  Vertical  Trend 

One  significant  reason  for  this  improved  correlation  in  the  higher  regions  of  the 
ship  could  be  due  to  the  fact  that  many  of  the  sensors  used  in  the  comparisons  in  these 
regions  were  accelerometers.  Comparisons  made  using  integrated  accelerometer  data 
showed  significantly  better  correlation  between  the  predicted  and  measured  response  than 
those  comparisons  made  using  velocity  meter  data.  As  shown  in  several  of  the  plots, 
much  of  this  error  can  be  contributed  to  the  drift  in  the  velocity  meter  data.  While  some 
integrated  accelerometer  plots  required  minor  drift  correction,  the  majority  of  the  plots 
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showed  no  apparent  drift.  In  sharp  contrast  to  this,  all  of  the  velocity  meter  plots  showed 
substantial  late  time  drift  and  required  considerable  manipulation  of  the  data  to  make  a 
valid  comparison  with  the  predicted  velocity  response.  Table  23  shows  the  average 
comprehensive  errors  for  each  shot  based  on  whether  the  measured  data  was  from  an 
accelerometer  or  a  velocity  meter.  Each  value  is  an  average  of  all  three  meshes  compared 
for  each  shot. 


Table  23.  Accelerometer  vs.  Velocity  Meter  Correlation  Comparison 


SENSOR 

AVERAGE  COMPREHENSIVE  RUSSELL’S  ERROR  FACTOR 

Shot  1 

Shot  2 

Shot  3 

Overall 

Accelerometer 

0.1845 

0.1434 

0.1910 

0.1730 

Velocity  Meter 

0.2269 

0.2196 

0.2315 

0.2260 

As  shown  in  Table  23,  the  difference  in  correlation  between  accelerometer  data 
and  velocity  meter  data  is  significant  with  an  overall  average  comprehensive  error  factor 
difference  of  over  0.05.  It  can  be  argued  that  at  least  some  of  this  improvement  can  be 
contributed  to  many  of  these  accelerometers  being  mounted  higher.  While  this  is  true, 
Tables  17,  20,  and  22  show  that  in  other  regions  of  the  ship  where  accelerometers  were 
used  there  is  significant  improvement  over  the  surrounding  velocity  meter  correlations. 

Overall,  using  the  comparison  of  the  fluid  mesh  displaying  the  best  correlation  for 
each  shot,  the  average  Russell’s  magnitude,  phase,  and  comprehensive  error  factors  were 
-0.0401,  0.2051,  and  0.1995  respectively.  Out  of  a  total  of  90  comparisons  made  for  all 
three  shots,  84  (93.3%)  predictions  were  considered  acceptable,  11  (12.2%)  were 
considered  excellent,  and  6  (6.7%)  were  considered  unacceptable.  Shot  2  showed  the 
best  correlation  (RCmean  =  0.1912)  overall  while  Shot  3  showed  the  least  (RCmean  = 
0.2114). 

F.  EXTENDED  DEPTH  FLUID  MODEL 

In  Part  B  of  this  chapter,  it  was  shown  that,  while  a  cavitation  depth  model 
yielded  excellent  predictions,  a  model  having  a  depth  twice  the  cavitation  depth  yielded 
poor  predictions.  The  comparison  error  factors  doing  a  120  msec  comparison  using  this 
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2X-cavitation  depth  model  are  repeated  and  averaged  in  Table  24.  The  results  of  the 
cavitation  depth  model  comparison  were  shown  previously  in  Table  10. 


Table  24.  2X  Cavitation  Fluid  Mesh  Russell’s  Error  Comparison  (120  msec) 


Sensor 

RM 

RP 

RC 

V2000V 

-0.6022 

0.3148 

0.6022 

V2002V 

-0.3407 

0.1664 

0.3360 

V2007V 

-0.4032 

0.2202 

0.4071 

V2008VI 

-0.3683 

0.2736 

0.4066 

V2009VI 

-0.3466 

0.1847 

0.3481 

V2010V 

-0.3472 

0.1959 

0.3533 

V2011VI 

-0.3287 

0.1901 

0.3365 

V2012VI 

-0.3097 

0.1485 

0.3044 

V2013V 

-0.3048 

0.1543 

0.3027 

V2014V 

-0.2080 

0.1621 

0.2337 

V2016V 

-0.2350 

0.1687 

0.2564 

V2018VI 

-0.3806 

0.2692 

0.4132 

V2019V 

-0.3924 

0.2887 

0.4317 

V2020V 

-0.4265 

0.2606 

0.4430 

V2026V 

-0.2738 

0.2147 

0.3084 

V2032V 

-0.4255 

0.1587 

0.4024 

V2034V 

-0.3917 

0.1392 

0.3684 

V2035V 

0.4147 

0.1442 

0.3891 

V2108V 

-0.3856 

0.1992 

0.3846 

V2124V 

-0.3149 

0.1734 

0.3186 

V2125V 

-0.3355 

0.2027 

0.3474 

Average 

-0.3193 

0.2014 

0.3664 

These  two  models  compared  previously  had  the  same  number  of  hull  normal 
elements.  This  resulted  in  the  2X-cavitation  depth  model  having  a  hull  normal  spacing 
between  elements  in  the  outer  mesh  twice  that  of  the  cavitation  depth  model.  As  has 
been  shown  in  the  previous  comparisons,  this  large  of  variation  in  the  mesh  density 
would  certainly  have  an  effect  on  the  resulting  predictions. 

To  test  whether  or  not  increasing  the  mesh  density  of  this  2X-cavitation  depth 
model  would  improve  the  prediction  correlation,  a  similar  test  was  conducted  using  the 
3ld  generation  model.  This  model  was  constructed  to  twice  the  actual  cavitation  depth  but 
having  a  normal  fluid  element  spacing  equal  to  the  dense  cavitation  depth  model 
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compared  earlier.  The  resulting  120  msec  error  comparison  for  this  model  is  shown 
below  in  Table  21.  This  comparison  shows  only  the  sensors  compared  in  Reference  [21]. 


Table  25.  2X  Cavitation  Depth  Fluid  Mesh  Russell’s  Error  Comparison  (120  msec) 


Sensor 

RM 

RP 

RC 

V2000V 

-0.4389 

0.2537 

0.4492 

V2002V 

-0.1658 

0.2096 

0.2368 

V2007V 

-0.3123 

0.2393 

0.3487 

V2008VI 

-0.2991 

0.2304 

0.3346 

V2009VI 

-0.2467 

0.1875 

0.2746 

V2010V 

-0.3325 

0.1950 

0.3416 

V2011VI 

-0.3105 

0.1649 

0.3116 

V2012VI 

-0.2939 

0.1500 

0.2925 

V2013V 

-0.2929 

0.1417 

0.2884 

V2014V 

-0.2436 

0.1790 

0.2679 

V2016V 

-0.2643 

0.1649 

0.2760 

V2018VI 

-0.3270 

0.1535 

0.3201 

V2019V 

-0.3503 

0.1795 

0.3489 

V2020V 

-0.3600 

0.1945 

0.3626 

V2026V 

-0.2825 

0.2074 

0.3106 

V2032V 

-0.3472 

0.1742 

0.3443 

V2034V 

-0.2974 

0.1637 

0.3008 

V2035V 

-0.3268 

0.1958 

0.3376 

V2108V 

-0.2996 

0.1749 

0.3074 

V2124V 

-0.2938 

0.1529 

0.2935 

V2125V 

-0.3055 

0.1674 

0.3088 

Average 

-0.30431 

0.18475 

0.31698 

It  should  be  noted  that  because  this  model  had  a  depth  equal  to  twice  the 
minimum  cavitation  depth,  its  depth  was  17  feet  less  than  the  model  used  in  the  previous 
study.  Additionally,  as  shown  previously,  the  3rd  generation  model  displayed  measurable 
improvement  over  the  1st  generation  model.  However,  the  previous  study  used  fdtered 
data  which  was  shown  to  reduce  the  error  factors  significantly.  Assuming  that  these 
variations  essentially  canceled  each  other  out,  the  resulting  error  was  significantly 
reduced  using  the  denser  fluid  mesh. 

The  resulting  error  factors  while  significantly  reduced,  were  still  significantly 
higher  than  the  comparable  cavitation  depth  model.  A  side  by  side  500  msec  error 

comparison  of  these  two  models  is  shown  in  Table  26. 
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Table  26.  Russell’s  Error  Comparison:  Cavitation  Depth  Model  vs.  2X  Cavitation 

Depth  Model  (500  msec) 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Mesh:  Cavitation  Depth  i 

1  Dense  Mesh:  2x  Cavitation 

X 

y 

z 

RM 

RP 

RC 

RM 

RP 

RC 

V2000V 

120217 

Keel 

5328 

m 

0.2979 

-0.4112 

0.3215 

V2002V 

142489 

Keel 

4656 

m 

0.2150 

-0.1692 

0.2354 

V2007V 

210430 

Keel 

4080 

m 

0.2766 

-0.3631 

0.3075 

V2009VI 

210808 

Bulkhead 

4080 

0.1996 

-0.3088 

0.2332 

V2008VI 

210894 

Bulkhead 

4080 

0.2270 

-0.3268 

0.2764 

A8516V 

210993 

Bulkhead 

4080 

0.1217 

0.2502 

-0.2536 

0.3552 

A4701V 

211371 

Bulkhead 

4080 

0.1921 

-0.3064 

0.2283 

V2107V 

211904 

Bulkhead 

4080 

-0.0295 

0.2161 

-0.2907 

0.3048 

V2108V 

212196 

Bulkhead 

4080 

0.2348 

0.2089 

-0.3463 

0.2387 

V2010V 

220589 

Keel 

3504 

0.2549 

0.2730 

-0.3880 

0.2768 

V2012VI 

221102 

Bulkhead 

3504 

-0.0777 

0.1935 

0.1848 

-0.3195 

0.2137 

V2011VI 

221188 

Bulkhead 

3504 

0.2398 

0.2468 

-0.3784 

0.2799 

0.4171 

V2013V 

221601 

Bulkhead 

3504 

-0.0342 

0.1970 

0.1772 

-0.3202 

0.2338 

0.3514 

V2124V 

222060 

Bulkhead 

3504 

-0.0872 

0.1882 

0.1838 

-0.3196 

0.1945 

0.3296 

A2104V 

222240 

Bulkhead 

3504 

0.1760 

0.1571 

-0.2873 

0.1712 

0.2964 

V2125V 

222436 

Bulkhead 

3504 

0.2264 

0.2298 

-0.3494 

0.2673 

0.3898 

V2014V 

230461 

Keel 

2952 

0.2501 

0.2344 

-0.2798 

0.2355 

0.3194 

A3565V 

231696 

Bulkhead 

2952 

0.0646 

0.2003 

0.1865 

-0.1765 

0.1842 

0.2261 

V2016V 

242399 

Keel 

2544 

0.2085 

0.2005 

-0.2804 

0.1971 

0.3037 

V2026V 

312302 

Keel 

1992 

0.2312 

-0.2819 

0.2225 

0.3183 

A2310V 

320746 

Bulkhead 

1536 

0.1635 

0.1716 

-0.3346 

0.1806 

0.3370 

V2034V 

330759 

Keel 

1152 

0.1537 

0.1608 

-0.3223 

0.1813 

0.3277 

V2032V 

330764 

Keel 

1152 

0.2103 

0.2307 

-0.3700 

0.3740 

V2035V 

330769 

Keel 

1152 

0.2233 

0.2434 

-0.3755 

0.3947 

V2019V 

340167 

Keel 

672 

0.2335 

-0.3444 

0.3703 

V2018VI 

340992 

Bulkhead 

672 

-0.1212 

0.1945 

-0.3204 

0.3319 

V2020V* 

350052 

Keel 

288 

0.2343 

-0.3804 

0.3877 

A2413V 

350220 

Bulkhead 

288 

-0.1315 

0.2233 

-0.3701 

0.3795 

A2116V 

414367 

Bulkhead 

3504 

0.1616 

-0.2645 

0.2753 

A2109V 

414953 

Bulkhead 

4059 

0.2028 

-0.1999 

0.2466 

A2240V 

416269 

Bulkhead 

3504 

-135 

848 

0.1552 

-0.2610 

0.2679 

A2237V 

416419 

Bulkhead 

3504 

E1EI 

0.1750 

-0.2626 

Sum(E(X)) 

-1.82320 

6.90470 

6.76420 

-9.96280 

7.33730 

11.02690 

Sum(E(XA2)) 

0.34947 

1.52398 

1.47020 

3.21065 

1.75589 

3.89648 

Mean 

-0.05698 

0.21577 

0.21138 

-0.31134 

0.22929 

0.34459 

Variation 

0.00792 

0.00110 

0.00130 

0.00351 

0.00237 

0.00312 

*350msec  comparison  due  to  long  time  sensor  data  inaccuracy. 


The  cavitation  depth  fluid  model  produced  significantly  better  predictions  than  the 
2X-cavitation  depth  model.  In  fact,  few  of  the  2X-cavitation  depth  model  predictions 
even  met  the  acceptable  criteria.  Though  displaying  significant  amounts  of  error  both  in 
phase  and  in  magnitude,  the  majority  of  the  comprehensive  error  for  the  2X-cavitation 
depth  model  is  due  to  the  magnitude  error.  This  is  in  sharp  contrast  to  the  cavitation 
depth  model  where  the  majority  of  the  comprehensive  error  is  due  to  the  phase  error.  The 
average  phase  error  of  the  2X-cavitation  fluid  model  prediction  is  only  slightly  greater 


93 


than  the  cavitation  depth  fluid  model  prediction.  By  comparison,  the  magnitude  error  of 
the  2X-cavitation  model  is  nearly  six  times  that  of  the  cavitation  depth  model. 

Of  additional  interest  is  the  relatively  small  increase  in  the  overall  error  of  the  500 
msec  comparison  compared  with  the  120  msec  comparison.  Here  the  majority  of  the 
increase  in  the  comprehensive  error  is  due  to  the  increase  in  the  phase  error.  This 
indicates  that  the  majority  of  the  magnitude  error  occurs  in  the  early  time  response,  where 
the  peak  response  occurs.  To  help  visualize  this,  Figures  62  and  63  show  the  velocity 
plot  comparisons  for  sensors  A2104V  and  V2034V  respectively. 

SHOT  3 

Grid  222240-vz  (A2104V) 


Bulkhead  (x=3504  y=0  z=39()) 


Time  (msec) 

Ship  Shock:  Trial  - Simulation  (lx  Cavitation) 

Simulation  (2x  Cavitation) 


Figure  62.  Bulkhead  Sensor  A2104V 
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SHOT  3 

Grid  330759-vz  (V2034V) 


KLeel  (x=l  152  y=-  1  3 5  z=193) 


Ship  Shock:  Trial  -  Simulation  (lx  Cavitation) 

Simulation  (2x  Cavitation) 


Figure  63.  Keel  Sensor  V2034V 

As  is  shown  in  these  plots,  the  larger  fluid  mesh  prediction  significantly 
underestimates  the  initial  kick  and  recovery  of  the  structure.  This  portion  of  the  response 
is  due  primarily  to  the  incident  pressure  wave.  This  gives  evidence  of  losses  occurring  in 
the  incident  pressure  wave  as  the  DAA  boundary  is  moved  further  away  from  the  hull. 
Dr.  John  DeRuntz,  author  of  the  USA  code,  suggests  that  these  losses  could  be  due  to  the 
extremely  large  fluid  mesh  causing  dispersion  of  the  shock  wave  front  [Ref.  28]. 

G.  LATERAL  AND  LONGITUDINAL  VELOCITY  COMPARISON 

While  the  majority  of  the  response  to  underwater  shock  occurs  in  the  vertical 

direction,  a  significant  portion  of  the  structural  response  can  occur  in  the  lateral 

(athwartship)  and  longitudinal  (fore  and  aft)  directions  depending  on  the  shot  geometry. 

Previous  simulations  using  LS-DYNA/USA  have  shown  poor  correlation  when 

comparing  lateral,  or  athwartship,  and  longitudinal,  or  fore  and  aft,  velocity  response 

[Ref.  9].  To  determine  how  well  simulations  using  the  DDG  81  model  coupled  with  the 

fluid  models  predicted  these  responses  1 1  comparisons  in  the  lateral  direction  and  6 

95 


comparisons  in  the  longitudinal  direction  were  perfonned  using  measured  data  from  Shot 
3.  All  sensors  used  in  these  comparisons  were  accelerometers.  The  locations  of  these 
accelerometers  are  shown  in  Tables  27  and  28. 


Table  27.  Lateral  (Athwartship)  Accelerometer  Nodal  Locations 


Sensor 

NODE 

X(Node) 

Y(Node) 

Z(Node) 

Compartment/Location 

Compt#/Area 

A2001A 

120217 

5328.0 

0.0 

82.0 

PASSAGE 

4-022-0-L 

A2102A 

212058 

4080.0 

0.0 

390.0 

CIC  PROJECTION  ROOM 

1-126-0-C 

A2110A 

414953 

4059.1 

0.0 

722.8 

RR  #1 

03-128-0-C 

A2117AI 

414367 

3504.0 

0.0 

702.0 

RR  #2 

03-142-0-C 

A2238AI 

416419 

3504.0 

135.0 

848.0 

PORT  MAST  LEG 

MAST 

A2241A 

416269 

3504.0 

-135.0 

848.0 

STBD  MAST  LEG 

MAST 

A2105AI 

222240 

3504.0 

0.0 

390.0 

CIC  ANNEX 

1-126-0-C 

A2123A 

232685 

2952.0 

0.0 

502.0 

ACCESS  TRUNK 

1-220-0-T 

A2015AI 

230461 

2952.0 

0.0 

85.5 

AMR  #2 

4-220-0-E 

A2311A 

320746 

1536.0 

0.0 

177.0 

A/C  MCHY  &  PUMP  RM 

5-300-01-E 

A2033A 

330764 

1152.0 

0.0 

196.0 

FAN  ROOM 

3-362-0-E 

A2021A 

350052 

288.0 

0.0 

211.0 

STEERING  GEAR  RM 

4-442-0-E 

Table  28.  Longitudinal  (Fore-Aft)  Accelerometer  Nodal  Locations 


Sensor 

NODE 

X(Node) 

Y(Node) 

Z(Node) 

Compartment/Location 

Compt#/Area 

A2103F 

212058 

4080.0 

0.0 

390.0 

CIC  PROJECTION  ROOM 

1-126-0-C 

A2111F 

414953 

4059.1 

0.0 

722.8 

RR  #1 

03-128-0-C 

A2118F 

414367 

3504.0 

0.0 

702.0 

RR  #2 

03-142-0-C 

A2239F 

416419 

3504.0 

135.0 

848.0 

PORT  MAST  LEG 

MAST 

A2242F 

416269 

3504.0 

-135.0 

848.0 

STBD  MAST  LEG 

MAST 

A2106F 

222240 

3504.0 

0.0 

390.0 

CIC  ANNEX 

1-126-0-C 

1.  Lateral  Velocity  Response 

Table  29  lists  the  error  factors  for  each  of  the  1 1  lateral  velocity  comparisons.  As 
is  shown,  none  of  the  comparisons  showed  good  correlation.  None  of  the  comparisons 
yielded  a  comprehensive  error  factor  below  0.30.  Errors  in  both  phase  and  magnitude 
were  significant.  The  simulation  predictions  underestimated  the  response  magnitude  in 
all  cases. 

Similar  to  the  vertical  velocity  comparisons,  both  the  furthest  forward  and  furthest 
aft  sensors  displayed  the  least  correlation.  Though  none  of  the  fluid  models  yielded 
acceptable  correlation,  the  dense  fluid  mesh  model  did  show  a  significant  improvement 
over  the  medium  fluid  mesh  model.  The  coarse  mesh  showed  measurably  less  correlation 
than  the  medium  mesh. 
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Table  29.  Russell’s  Error  Factor  for  Lateral  Velocity  Response  Comparisons 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Mesh  | 

|  Medium  Mesh  I 

]  Coarse  Mesh 

X 

z 

RM 

RP 

RC 

RM 

RP 

RC 

RM 

RP 

RC 

A2001A 

120217 

Keel 

El 

-0.3568 

0.4254 

0.4921 

-0.4251 

0.4039 

0.5197 

-0.4587 

0.4203 

0.5513 

A2102A 

212058 

Bulkhead 

O 

-0.1838 

0.2899 

0.3042 

-0.2130 

0.3189 

0.3399 

-0.2265 

0.3223 

0.3491 

A2110A 

414953 

Bulkhead 

o 

-0.1760 

0.3575 

0.3532 

-0.2097 

0.3798 

0.3845 

-0.2156 

0.3889 

0.3941 

A2105AI 

222240 

Bulkhead 

o 

-0.2601 

0.3485 

0.3854 

-0.2663 

0.3746 

0.4073 

-0.2691 

0.3725 

0.4072 

A2117AI 

414367 

Bulkhead 

o 

-0.3132 

0.3175 

0.3952 

-0.3142 

0.3485 

0.4158 

-0.3084 

0.3551 

0.4168 

A2238AI 

416419 

Bulkhead 

-0.3312 

0.3548 

0.4302 

-0.3316 

0.3815 

0.4480 

-0.3248 

0.3883 

0.4487 

A2241A 

416269 

Bulkhead 

-0.3113 

0.3366 

0.4063 

-0.3082 

0.3658 

0.4239 

-0.2954 

0.3741 

0.4224 

A2015AI 

230461 

Keel 

mm 

-0.2241 

0.3200 

0.3462 

-0.2403 

0.3245 

0.3578 

-0.2428 

0.3212 

0.3568 

A2311A 

320746 

Bulkhead 

o 

-0.1668 

0.3080 

0.3104 

-0.1794 

0.3038 

0.3127 

-0.1937 

0.3040 

0.3194 

A2033A 

330764 

Keel 

mm 

-0.2036 

0.3764 

0.3793 

-0.2063 

0.3780 

0.3816 

-0.2102 

0.3771 

0.3826 

A2021A 

350052 

Keel 

o 

-0.3904 

0.4259 

0.5120 

-0.4054 

0.4261 

0.5212 

-0.4226 

0.4295 

0.5340 

Sum(E(X)) 

-2.91730 

3.86050 

4.31450 

-3.09950 

4.00540 

4.51240 

-3.16780 

4.05330 

4.58240 

Sum(E(XA2)) 

0.83632 

1.37459 

1.73653 

0.94147 

1.47225 

1.89571 

0.98772 

1.50961 

1.96127 

Mean 

-0.26521 

0.35095 

0.39223 

-0.28177 

0.36413 

0.41022 

-0.28798 

0.36848 

0.41658 

Variation 

0.04840 

0.07657 

0.09743 

0.05440 

0.08175 

0.10623 

0.05721 

0.08390 

0.11012 

Figures  64  through  69  show  6  of  the  1 1  lateral  velocity  500  msec  comparison 
plots.  The  remaining  five  lateral  velocity  comparison  plots  are  included  in  Appendix  A. 
Figure  64  shows  the  comparison  for  sensor  A2001A  located  above  the  sonar  dome. 
Previous  comparisons  in  this  location  have  shown  poor  correlation.  This  trend  continues 
with  the  lateral  velocity  comparison.  The  simulation  prediction  in  this  case  drastically 
underestimates  the  magnitude  of  the  response.  Interestingly,  while  the  peak  lateral 
velocity  measured  was  approximately  6  ft/sec,  the  corresponding  vertical  velocity 
measured  was  only  4.8  ft/sec  indicating  that  the  primary  response  is  in  the  lateral 
direction  at  this  location. 

Figure  65  shows  the  comparison  of  sensor  A2015AI  located  at  midships.  The 
predicted  response  accurately  predicts  the  initial  peak  magnitude,  but  the  error  grows 
considerably  during  the  initial  late  time  response  despite  showing  similar  damping 
characteristics  after  200  msec. 

Moving  further  aft  is  the  velocity  plot  comparison  for  bulkhead  sensor  A2311A 
shown  in  Figure  66.  As  in  the  sensor  A2015AI  comparison,  the  initial  peak  velocity  is 
captured  quite  well  by  the  simulations.  The  majority  of  the  poor  correlation  is  due  to  the 
response  after  this  peak  where  the  simulation  predicted  response  is  out  of  phase  and  the 
magnitude  of  the  measured  late  time  response  actually  increases  while  the  predicted 
response  continues  to  decay. 
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Figure  64.  Keel  Sensor  A2001A 

SHOT  3 

Grid  23046  1  -vy  (A2015AI) 
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Time  (msec) 

Ship  Shock  Trial  -  Simulation  (COARSE) 

Simulation  (MEDITJIVl)  -  Simulation  (DENSE) 

Figure  65.  Keel  Sensor  A2015AI 
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SHOT  3 

Grid  320746-vy  (A231  1A) 


Bulkhead  (x=1536  y=0  z=177) 


Ship  Shock:  Trial  - Simulation  (COARSE) 

Simulation  (MEDIUM)  - Simulation  (DENSE) 


Figure  66.  Bulkhead  Sensor  A2311A 

The  furthest  aft  sensor  compared  was  keel  sensor  A2021A  shown  in  Figure  67. 
As  with  the  furthest  forward  sensor,  little  correlation  between  the  predicted  and  measured 
response  is  shown  throughout  the  response  time  compared.  The  simulation  prediction 
underestimates  the  early  time  peak  velocity  and  poor  phase  correlation  exists  during  the 
remainder  of  the  response  time.  The  measured  response  shows  significant  drift 
contributing  to  the  error  during  the  late  time  response. 

Figure  68  shows  the  comparison  for  bulkhead  sensor  A2105AI  located  forward  of 
midships  and  higher  in  the  structure  above  the  waterline.  The  magnitude  of  the  predicted 
response  is  comparable  to  the  measured  response,  although  the  existence  of  high 
frequencies  in  the  predicted  response  which  are  not  present  in  the  measured  response  is 
evident.  Though  out  of  phase,  the  predicted  response  shows  very  good  late  time 
correlation  to  the  measured  response  in  terms  of  damping. 

Sensor  A2238AI  located  on  the  port  mast  leg  is  shown  in  Figure  69.  As  with 
sensor  A2105AI,  the  predicted  response  contains  higher  frequencies  than  the  measured 
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Athwart  Velocity  (ft/s) 


response.  The  simulation  prediction  accurately  predicts  the  early  time  response,  but  is 
largely  out  of  phase  for  the  remainder  of  the  response  compared. 
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Figure  67.  Keel  Sensor  A2021A 
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Figure  68.  Bulkhead  Sensor  A2105AI 
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SHOT  3 

Grid  416419-vy  (A2238AI) 


Mast  (x=3504  y=135  z=848) 


Time  (msec) 

- Ship  Shock:  Trial  - Simulation  (COARSE) 

- Simulation  (TvlEDIUTvI)  -  Simulation  (DENSE) 

Figure  69.  Mast  Sensor  A2238AI 
2.  Longitudinal  Velocity  Response 

.Unlike  the  vertical  and  lateral  comparisons,  only  a  small  number  of  sensors  in  the 
same  general  proximity  were  provided  for  the  longitudinal  velocity  comparisons.  Table 
30  summarizes  the  longitudinal  velocity  comparisons  performed  using  the  six  sensors  for 
which  data  was  provided.  As  with  the  lateral  velocity  comparison,  the  longitudinal 
velocity  correlation  is  generally  poor  overall.  The  dominant  source  of  error  in  all  the 
longitudinal  comparisons  is  in  the  phase  although  there  is  significant  magnitude  error  in 
most  of  the  cases. 


Table  30.  Russell’s  Error  Factor  for  Longitudinal  Velocity  Response  Comparisons 


Sensor 

NODE 

Bulkhead 

Keel 

Location  (in) 

Dense  Mesh 

Medium  Mesh 

Coarse  Mesh 

X 

V 

z 

RM 

RP 

RC 

RM 

RP 

RC 

RM 

RP 

RC 

A2103F 

212058 

Bulkhead 

4080 

0 

390 

0.2939 

0.3728 

0.4207 

0.2889 

0.3756 

0.4199 

0.3049 

0.3831 

0.4339 

A2111F 

414953 

Bulkhead 

4059 

0 

723 

-0.0820 

0.4068 

0.3678 

-0.1638 

0.3999 

0.3830 

-0.1472 

0.4004 

0.3781 

A2106F 

222240 

Bulkhead 

3504 

0 

390 

0.3068 

0.3791 

0.4322 

0.3135 

0.3762 

0.4340 

0.3366 

0.3700 

0.4433 

A2118F 

414367 

Bulkhead 

3504 

0 

702 

0.3102 

0.3776 

0.4331 

0.2484 

0.3730 

0.3972 

0.2623 

0.3740 

0.4048 

A2239F 

416419 

Bulkhead 

3504 

135 

848 

0.1800 

0.3719 

0.3662 

0.1230 

0.3744 

0.3493 

0.1413 

0.3721 

0.3527 

A2242F 

416269 

Bulkhead 

3504 

-135 

848 

0.2771 

0.3745 

0.4129 

0.2073 

0.3522 

0.3621 

0.2257 

0.3502 

0.3692 

Sum(E(X)) 

1.28600 

2.28270 

2.43290 

1.01730 

2.25130 

2.34550 

1.12360 

2.24980 

2.38200 

Sum(E(XA2)) 

0.39264 

0.86932 

0.99123 

0.32838 

0.84587 

0.92226 

0.36764 

0.84496 

0.95231 

Mean 

0.21433 

0.38045 

0.40548 

0.16955 

0.37522 

0.39092 

0.18727 

0.37497 

0.39700 

Variation 

0.05560 

0.11387 

0.12997 

0.04857 

0.11081 

0.12096 

0.05376 

0.11070 

0.12495 

01 


Though  there  are  few  comparisons  and  the  correlations  of  all  of  the  fluid  mesh 
simulations  were  relatively  poor,  the  medium  density  fluid  mesh  showed  the  best 
correlation,  followed  by  the  coarse  mesh,  and  then  the  dense  mesh.  Figures  70  through 
75  show  the  longitudinal  velocity  comparisons  for  all  six  sensors  compared.  It  is 
important  to  note  the  magnitude  of  the  longitudinal  response  when  analyzing  these  plots. 
In  comparison  to  the  vertical  response,  the  longitudinal  response  is  an  order  of  magnitude 
smaller.  For  this  reason,  very  minor  errors  in  the  magnitude  of  the  simulation  predictions 
are  magnified. 

Figure  70  shows  the  longitudinal  velocity  comparison  plot  for  sensor  A2103F 
located  forward  of  midships,  1 1  feet  above  the  waterline.  This  plot  shows  how  the 
simulation  predictions  over  estimate  the  magnitude  of  the  response  for  the  majority  of  the 
response  time  compared.  Despite  being  dramatically  dissimilar  initially,  the  predicted 
and  measured  responses  appear  to  die  out  in  a  similar  fashion  indicative  of  accurate 
damping  in  the  structural  model. 


SHOT  3 

Grid  21205  8-vx  (A2103F) 


Bulkhead  (x=4080  y=0  z=390) 


Time  (msec) 

Ship  Shock:  Trial  - Simulation  (COARSE) 

Simulation  (MEDIUM)  - Simulation  (DENSE) 

Figure  70.  Bulkhead  Sensor  A2103F 
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The  comparison  of  bulkhead  sensor  A21 1  IF  is  shown  in  Figure  71,  located  higher 
up  in  the  structure  and  slightly  aft  of  sensor  A2103F.  Unlike  for  sensor  A2103F,  the 
simulation  predictions  for  this  sensor  underestimated  the  magnitude  of  the  longitudinal 
velocity  during  the  early  time  response.  This  is  the  only  longitudinal  velocity  comparison 
in  which  the  simulation  predictions  underestimated  the  velocity  magnitude.  Similar  to 
the  previous  comparison,  the  predicted  and  measure  response  decay  identically. 


SHOT  3 

Grid  4  1 4953-vx  (A21  1  IF) 


j  Time  (msec) 

Shock:  Trial  Data  Simulation  (COARSE)  I 

-  Simulation  (MEDIUM)  -  Simulation  (DENSE) 

Figure  71.  Bulkhead  Sensor  A2111F 

Figures  72  through  75  show  the  comparisons  of  sensors  located  just  forward  of 
midships.  In  all  of  these  cases,  the  simulation  predictions  overestimate  the  velocity 
magnitude  during  the  initial  response  time.  As  in  the  previous  two  comparisons,  the 
measured  and  predicted  responses  die  out  identically.  Although  the  correlation  is  poor, 
the  damping  coefficients  used  for  the  model  appear  to  be  accurate. 
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Figure  72.  Bulkhead  Sensor  A2106F 
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Figure  73.  Bulkhead  Sensor  A2118F 
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Longitudinal  Velocity  (ft/s)  Longitudinal  Velocity  (ft/s) 
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Figure  74.  Bulkhead  Sensor  A2242F 
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Figure  75.  Bulkhead  Sensor  A2239F 
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VII.  Conclusions  and  Recommendations 


Creating  a  virtual  underwater  shock  environment  has  proven  to  be  a  valid  method 
of  predicting  the  response  of  a  surface  ship  to  a  variety  of  underwater  explosions.  Three 
separate  underwater  explosions  were  investigated  in  this  report,  all  of  which  produced 
acceptable  vertical  velocity  response  predictions  overall.  These  underwater  explosions 
originated  from  abeam  the  ship  to  off  the  bow,  slightly  forward  of  the  ship.  While  the 
accuracy  of  these  predictions  falls  off  in  regions  of  the  ship  furthest  from  the  shot,  the 
predictions  in  these  regions  still  remain  acceptable  with  only  a  few  exceptions.  The 
accuracy  of  the  simulation  predictions  was  shown  to  have  little  sensitivity  to  stand-off 
distance. 

The  depth  of  the  fluid  model  had  the  largest  impact  on  the  accuracy  of  the 
predictions.  A  fluid  model  having  a  depth  equal  to  the  actual  cavitation  depth  at  the 
region  closest  to  the  charge  resulted  in  the  most  accurate  predictions.  Small  deviations  of 
only  a  few  feet  significantly  impacted  the  accuracy  of  the  predictions. 

Similarly,  a  fluid  model  having  a  depth  equal  to  twice  the  cavitation  depth  yielded 
unacceptable  predictions,  significantly  worse  than  the  cavitation  depth  fluid  model.  This 
disparity  appears  to  be  primarily  due  to  losses  occurring  in  the  incident  pressure  wave 
front.  Fluid  models  having  depths  significantly  less  than  the  cavitation  depth  produced 
unacceptable  results  as  well. 

The  fluid  model  mesh  density  was  shown  to  have  a  significant  effect  on  the 
simulation  predictions.  For  shocks  originating  from  abeam  the  ship,  a  fluid  mesh  having 
a  higher  mesh  density  was  shown  to  provide  substantially  better  predictions.  For  the  shot 
off  the  bow,  closer  to  being  an  end-on  shot  than  a  beam  shot,  the  coarser  mesh  was 
shown  to  produce  the  most  accurate  predictions.  It  is  theorized  that  this  disparity  is  due 
to  the  tangential  propagation  of  the  shock  wave  and  the  resultant  fluid  particle  velocities. 
In  this  direction,  the  coarser  mesh  has  more  preferable  fluid  element  aspect  ratios  than  the 
higher  density  meshes. 

Small  variations  in  the  quality  of  the  fluid  mesh  were  shown  to  have  measurable 
but  minor  effects  on  the  accuracy  of  the  overall  simulation  predictions.  Although  this 
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study  did  not  investigate  the  effects  of  gross  variations  in  the  quality  of  the  mesh,  minor 
improvements  in  the  quality  of  the  mesh  were  shown  to  result  in  overall  improvement, 
but  mixed  results  locally.  Variables  such  as  the  size  of  the  fluid  mesh,  the  quality  of  the 
measured  data,  and  environmental  factors  appear  to  outweigh  any  benefit  gained  by 
making  small  improvements  in  the  quality  of  the  mesh. 

The  correlation  of  the  simulation  predictions  was  found  to  be  substantially  better 
when  being  compared  with  accelerometer  measured  data  than  with  velocity  meter 
measured  data.  Substantial  drift  is  apparent  in  nearly  all  of  the  velocity  meter  data.  Data 
collected  from  accelerometer  displayed  significantly  less  drift  even  after  being  integrated. 
This  suggests  that  a  substantial  amount  of  the  calculated  error  may  be  due  not  to  error  in 
the  simulation  predictions,  but  in  the  measured  data  itself. 

The  study  conducted  in  this  paper  has  merely  helped  to  prop  open  the  door  for 
further  research  using  virtual  underwater  shock  environments.  By  validation  of 
underwater  explosion  simulations,  the  possibilities  for  the  future  use  of  virtual 
underwater  shock  environments  are  endless.  Some  recommendations  for  future  study 
include: 

1.  In  addition  to  producing  significantly  better  predictions  for  beam  shots, 
the  dense  fluid  mesh  simulations  required  considerably  more  time  to  process,  especially 
as  the  size  of  the  mesh  was  increased.  This  is  counterproductive  to  simulating  more  of 
the  overall  structural  response  to  underwater  explosions.  The  model  created  in  the  course 
of  this  study  is  easily  modified  to  model  a  wide  range  of  underwater  environments.  It  is 
recommended  that  this  model  be  used  to  study  the  effects  of  reducing  the  size  of  the 
charge  while  maintaining  the  same  shock  factor.  This  would  reduce  the  size  of  the  fluid 
model  required,  allowing  for  a  denser  mesh  to  be  used  without  having  a  large  time 
requirement  to  process  the  simulation. 

2.  This  study  used  models  having  a  depth  equal  to  the  cavitation  depth  at  the 
location  on  the  keel  closest  to  the  charge.  As  the  wave  propagates  outward  the  actual 
cavitation  depth  increases.  This  fluid  model  may  be  utilized  to  study  the  effects  of 
contouring  the  depth  of  the  fluid  mesh  to  correspond  with  the  actual  cavitation  depth 
throughout  the  length  of  the  ship. 


107 


3.  This  study  investigated  blue  water  underwater  explosions  where  bottom 
reflection  is  not  a  consideration.  The  evolution  of  the  existing  threats  in  recent  years  has 
placed  an  increased  emphasis  on  littoral  operations.  In  shallow  water  reflection  of  the 
shock  wave  off  the  ocean  bottom  and  the  bubble  pulse  cannot  be  ignored.  Recent 
enhancements  to  the  codes  utilized  in  this  study  allow  the  capability  to  simulate  time 
dependent  bottom  reflection  and  bubble  pulse.  It  is  recommended  that  the  simulations  be 
extended  out  in  time  to  predict  the  structural  response  to  these  phenomena. 

4.  Recent  studies,  including  the  study  conducted  in  this  paper,  have  focused 
primarily  on  the  velocity  and  acceleration  response  of  the  ship  in  order  to  validate  the  use 
of  underwater  shock  simulation.  The  purpose  for  subjecting  ship  systems  to  an 
underwater  shock  is  to  ensure  that  the  systems  can  survive  the  combat  shock 
environment.  It  is  suggested  that  underwater  shock  simulation  be  carried  a  step  further 
by  using  simulation  to  conduct  damage  surveys  of  ship  systems  when  tested  at  design 
level  using  more  severe  shock  simulations. 
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APPENDIX  A:  DDG  81  VELOCITY  RESPONSE  PLOTS 


A.  SHOT  1 

The  following  velocity  plots  are  from  the  Shot  1  fluid  mesh  density  comparison. 
The  Russell’s  error  factors  following  each  figure  caption  are  from  the  coarse  fluid  model 
comparison.  This  fluid  model  mesh  showed  the  best  correlation  overall  for  this  shot. 
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Figure  76.  Keel  Sensor  V2002V:  (RM=0.0772,  RP=0.1755,  RC=0.1699) 
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Figure  77. 
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Keel  Sensor  V2007V:  (RM=-0.1178,  RP=0.1760,  RC=0.1877) 
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Figure  78.  Bulkhead  Sensor  A2109V:  (RM=-0.0679,  RP=0.0848,  RC=0.0962) 
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Figure  79.  Bulkhead  Sensor  V2013V:  (RM=-0.1366,  RP=0.1819,  RC=0.2016) 
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Figure  80.  Bulkhead  Sensor  A2104V:  (RM=-0.2235,  RP=0.1184,  RC=0.2241) 
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Figure  81.  Bulkhead  Sensor  A2116V:  (RM=-0.0679,  RP=0.0848,  RC=0.0962) 
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Figure  82.  Bulkhead  Sensor  A2237V  (RM=-0.0906,  RP=0.1203,  RC=0.1334) 
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Figure  83.  Keel  Sensor  V2034V:  (RM=0.0296,  RP=0.2233,  RC=0.1996) 
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Figure  84. 


Keel  Sensor  V2019V:  (RM=-0.1463,  RP=0.2447,  RC=0.2526) 
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Figure  85.  Keel  Sensor  A2413V:  (RM=0.0171,  RP=0.3649,  RC=0.3238) 

B.  SHOT  2 

The  following  velocity  plots  are  from  the  Shot  2  fluid  mesh  density  comparison. 
The  Russell’s  error  factors  following  each  figure  caption  are  from  the  dense  fluid  model 
comparison.  This  fluid  model  mesh  showed  the  best  correlation  overall  for  this  shot. 
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Figure  87. 


Keel  Sensor  V2007V:  (RM=-0.0245,  RP=0.2092,  RC=0.1866) 
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Figure  88.  Bulkhead  Sensor  A2109V:  (RM=0.1170,  RP=0.1647,  RC=0.1790) 
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Figure  89.  Bulkhead  Sensor  V2013V:  (RM=-0.1049,  RP=0.1420,  ROO.1565) 
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Bulkhead  Sensor  A2104V:  (RM=-0.0388,  RP=0.0953,  RC=0.0912) 
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Figure  91.  Bulkhead  Sensor  A2116V:  (RM=0.0070,  RP=0.1003,  R00.0891) 
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Figure  92.  Mast  Sensor  A2240V:  (RM=-0.0626,  RP=0.1335,  RC=0.1307) 
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Figure  93.  Keel  Sensor  V2034V:  (RM=-0.1442,  RP=0.1874,  RC=0.2095) 
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Figure  94.  Keel  Sensor  V2019V:  (RM=-0.1327,  RP=0.2391,  ROO.2423) 
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Figure  95.  Keel  Sensor  A2413:  (RM=-0.1245,  RP=0.2219,  RC=0.2255) 
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C.  SHOT  3 

1.  Vertical  Velocity  Plots 

The  following  vertical  velocity  plots  are  from  the  Shot  3  fluid  mesh  density 
comparison.  The  Russell’s  error  factors  following  each  figure  caption  are  from  the  dense 
fluid  model  comparison.  This  fluid  model  mesh  showed  the  best  correlation  overall  for 
this  shot. 
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Figure  96.  Keel  Sensor  V2002V:  (RM=0.0762,  RP=0.2303,  RC=0.2150) 
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Figure  97.  Keel  Sensor  V2007V:  (RM=-0.1198,  RP=0.2883,  RC=0.2766) 
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Figure  98.  Bulkhead  Sensor  A2109V:  (RM=0.1353,  RP=0.1846,  RC=0.2028) 
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Figure  99.  Bulkhead  Sensor  V2013V:  (RM=-0.0342,  RP=0.1970,  RC=0.1772) 
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Figure  100.  Bulkhead  Sensor  A2104V:  (RM=0.0207,  RP=0.1760,  RC=0.1571) 
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101.  Bulkhead  Sensor:  (RM=0.0353,  RP=0.1789,  RC=0.1616) 
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Figure  102.  Mast  Sensor  A2237:  (RM=0.0613,  RP=0.1877,  R00.1750) 
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Figure  103.  Keel  Sensor  V2034V:  (RM=-0.0966,  RP=0.1537,  RC=0.1608) 
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Figure  104.  Keel  Sensor  V2019V:  (RM=-0.1212,  RP=0.1830,  RC=0.1945) 
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Figure  105.  Bulkhead  Sensor  A2413V:  (RM=-0.1315,  RP=0.2150,  ROO.2233) 
2.  Lateral  (Athwartship)  Velocity  Plots 
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Figure  106.  Bulkhead  Sensor  A2102A 
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Figure  107.  Bulkhead  Sensor  A2110A 
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Figure  108.  Bulkhead  Sensor  A2117AI 
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Figure  109.  Mast  Sensor  A2241A 
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Figure  110.  Keel  Sensor  A2033A 
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APPENDIX  B:  FLUID  MODELING  IN  TRUEGRID 


This  appendix  describes  the  complete  fluid  modeling  process  from  start  to  finish 
using  TrueGrid.  Most  of  the  commands  used  in  the  modeling  process,  many  of  which  are 
highlighted  in  bold  in  this  appendix,  have  been  explained  in  the  following  appendix, 
Appendix  C.  Further  guidance  is  given  in  Reference  [5].  The  intent  of  these  appendices 
is  not  to  serve  as  a  tutorial,  but  rather  as  a  guide  for  the  more  experienced  TrueGrid  user 
to  develop  a  similar  high  quality  fluid  model.  It  is  suggested  that  first  time  users  of 
TrueGrid  begin  by  using  one  of  several  tutorials  available  for  download  online  or  attend 
the  three  day  training  course  offered  by  XYZ  Scientific  Applications,  Inc.  More 
information  on  this  and  other  TrueGrid  services  provided  can  be  found  at 
www.truegrid.com. 


A.  STRUCTURAL  MODEL 

To  begin  the  fluid  modeling  process,  the  structural  model  must  first  be  imported 
into  TrueGrid.  TrueGrid  has  the  capability  to  import  NASTRAN,  NEUTRAL,  LS- 
DYNA,  and  DYNA3D  format  models  directly  [Ref.  5].  This  is  accomplished  using  the 
readmesh  command  followed  by  the  type  of  file,  the  filename,  the  exclude  command, 
and  then  endpart.  The  exclude  command  is  necessary  to  prevent  the  fluid  nodes  merged 
with  the  structural  nodes  to  be  redefined  when  the  fluid  mesh  is  output.  Since  the 
structural  model  is  being  defined  separately  in  LS-DYNA,  it  is  not  desirable  to  have 
duplicate  nodes.  Similarity,  since  the  structural  properties  are  being  defined  separately, 
all  materials  and  springs  should  be  deleted  from  the  model  to  avoid  being  output.  This  is 
done  using  the  delmats  and  delspds  commands  respectively. 

readmesh  nastran  ddg81shot3.nas  exclude  endpart 
merge 

delmats  1 :346; 
delspds  1:416; 
delspds  1 : 1 00; 
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Rectangular  voids  were  left  where  the  two  shafts  exit  the  hull  (stem  tubes)  on  the 
DDG  8 1  model.  Leaving  these  as  voids  would  create  an  extremely  challenging  surface  to 
merge  the  fluid  model  to.  To  fix  this  problem,  two  plates  are  constructed  to  cover  the 
voids,  leaving  a  smooth  surface  flush  with  the  remainder  of  the  hull.  Two  phases  are 
involved  in  the  TrueGrid  modeling  process,  the  part  phase  and  the  merge  phase.  The  part 
phase  is  where  each  part  is  constructed  separately  while  the  merge  phase  is  where  these 
parts  are  merged  together.  While  no  parts  are  constructed  in  the  merge  phase,  many  pre¬ 
cursors  to  building  new  parts,  such  as  defining  a  face  set,  are  done  in  this  phase. 
Construction  of  the  stern  tube  covers  is  quite  simple  and  is  done  in  the  block  part  phase. 
The  starboard  stern  tube  cover  is  shown  in  the  example  below. 

c  PART  2:  STARBOARD  STERN  TUBE  COVER 
block  1  8;-1;1  4;1 ,1520000e+03  1 ,4880000e+03  ; 

-1 35;9.9440002e+01  1 ,2893800e+02  ; 
pb  2  1  2  2  1  2  xyz  1488.0  -202.125  95.4375 

pb  2  1  1  2  1  1  xyz  1488.0  -135.0  60.875 

pb  1  1  2  1  1  2  xyz  1 152.0  -196.375  128.938 

pb  1  1  1  1  1  1  xyz  1 152.0  -135.0  99.440002 

endpart 

Following  completion  of  these  covers,  the  boundary  nodes  of  the  covers  should  be 
checked  to  ensure  that  they  have  merged  properly  with  the  original  hull  nodes. 
Improperly  merged  nodes  are  called  “cracks”.  It  is  advised  that  a  small  merging 
tolerance  such  as  0.25  inches  or  less  be  set  using  the  stp  command  during  the  early  stages 
of  modeling.  This  merging  tolerance  is  only  applicable  in  the  merge  phase.  No  tolerance 
is  recognized  in  the  part  phase.  This  results  in  cracks  removed  by  an  excessive  merge 
tolerance  reappearing  in  the  part  phase.  As  the  fluid  mesh  is  extended  outward,  cracks  in 
the  inner  layers  that  were  only  0.5  inch  or  less  may  grow  in  excess  of  several  inches.  If 
cracks  are  present,  the  part  can  be  reconstructed,  adding  more  partitions  to  attach  each 
individual  node  to  the  hull.  An  alternative  is  to  use  the  pn  command  and  attach  the 
unmerged  nodes  manually  in  the  merge  phase.  Although  this  command  is  not  to  be 
abused,  in  this  case  it  is  perfectly  acceptable. 
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B. 


INNER  LINER 


Construction  of  the  fluid  model  begins  first  by  wrapping  the  wetted  hull  structure 
with  two  individual  fluid  mesh  layers.  This  allows  the  proper  inner  liner  thickness  to  be 
ensured,  creates  transition  regions  at  the  corners  of  the  hull,  and  transitions  the  structural 
mesh  lines  to  the  test  waterline.  This  is  the  most  complex  portion  of  the  fluid  model, 
consisting  of  37  individual  parts  merged  together  into  the  final  two-layer  inner  liner. 

The  inner  layer  of  the  inner  liner  consists  of  32  parts;  the  port  and  starboard  lower 
hull  parts,  the  sonar  dome  part,  the  lower  and  upper  stern  parts,  the  port  and  starboard 
upper  hull  parts,  and  several  “seams”  which  connect  these  parts  into  one  continuous  fluid 
mesh  layer.  To  begin,  curves  must  be  created  using  the  lp3  command  along  the  keel 
seam,  the  sonar  dome  seams,  the  bow  seam,  and  the  stern  seams.  These  curves  must  then 
be  manually  modified  so  as  to  produce  curves  approximately  seven  inches  from  the  hull 
and  45°  off  from  the  original  curves  created.  These  curves  will  be  utilized  to  attach  the 
outer  edges  of  the  individual  parts  that  will  eventually  form  the  inner  layer  of  fluid  mesh 
on  the  hull.  Figure  111  shows  an  example  of  a  portion  of  these  curves. 
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Prior  to  construction  of  the  inner  mesh  layer,  the  face  elements  of  the  hull  must 
first  be  selected  using  the  manual  fset  method.  This  involves  using  the 
“pick;>sets;>faces”  option  in  the  environmental  window  of  TrueGrid.  The  port  and 
starboard  face  sets  are  separated  at  the  keel  and  bow  seam.  Only  the  face  elements  up  to 
three  mesh  lines  from  the  mesh  line  closest  to  the  test  waterline  on  the  ship  are  selected. 
The  sonar  dome  structural  face  elements  on  the  bottom  of  the  structure  and  the  lower 
stern  face  elements,  up  to  one  structural  mesh  line  from  the  test  waterline,  make  up 
separate  face  sets.  Three  face  sets  are  defined  to  create  each  upper  port  and  starboard 
fluid  mesh  layer.  These  face  sets  consist  of  the  structural  face  elements  not  picked  in  the 
lower  hull  face  sets  up  to  the  structural  mesh  line  closest  to  the  test  waterline.  This  is  to 
facilitate  transitioning  the  structural  mesh  lines  to  the  waterline  at  the  top  of  the  fluid 
mesh.  Figure  1 12  provides  some  clarification  as  to  exactly  how  these  face  sets  should  be 
defined.  Each  face  set  is  differentiated  by  color. 


Figure  112.  Face  Sets  Used  to  Generate  Inner  Fluid  Mesh  Layer  [from  Ref.  21], 
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The  port  and  starboard  lower  hull  mesh  layers  are  constructed  first  using  the 
corresponding  face  sets  defined.  This  is  accomplished  by  using  the  blude  command  to 
extrude  the  fluid  mesh  through  the  structural  mesh  and  the  bb  and  normal  commands  to 
ensure  the  mesh  is  orthogonal  to  the  structure  and  no  more  than  seven  inches  in  thickness. 
The  lower,  bow,  and  stem  outer  edges  of  the  mesh  are  attached  to  the  corresponding 
curves  generated  earlier.  This  process  is  repeated  for  the  sonar  dome  and  the  lower  stern 
region. 

To  maintain  orthogonality  and  fill  in  the  gaps  left  at  the  seams  between  these 
parts,  seam  fluid  mesh  parts  are  constructed  in  the  block  part  phase  to  fonn  a  wedge 
shaped  fluid  mesh  in  these  transition  regions.  For  the  keel  and  sonar  dome  seams  the 
wedge  shaped  parts  are  two  elements  thick  along  their  length  while  the  bow  seam  part  is 
four  elements  thick.  Pyramid  shaped  parts  are  necessary  to  properly  transition  the 
junction  of  the  sonar  dome  and  bow  seams  and  terminate  the  keel  seam  at  the  stern. 

These  are  constructed  in  the  block  part  phase  as  well.  Figure  1 13  shows  an  exploded 
view  of  the  extruded  parts  constructed  so  far  as  well  as  the  seam  parts. 


Figure  113. 


Inner  Mesh:  Lower  Port  and  Starboard  Hull,  Sonar  Dome,  and  Seams. 
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The  next  step  is  to  extend  the  inner  fluid  mesh  layer  to  the  waterline.  The  lower 
port  and  starboard  liner  was  only  constructed  up  to  three  elements  below  the  structural 
mesh  line  closest  to  the  waterline  because  this  structural  mesh  line  does  not  remain  at  a 
constant  height.  It  veers  upward  in  the  aft  section  of  the  ship,  extending  much  higher 
than  the  test  waterline.  This  requires  that  two  “step-downs”  must  be  taken  along  the 
structural  mesh.  To  aid  in  this  transition,  triangular  shell  elements  replaced  the 
quadrilateral  shell  elements  on  the  structural  model  on  each  side  at  two  locations  where 
they  could  transition  from  above  the  test  waterline  to  below  the  test  waterline.  This  only 
serves  as  an  aid  in  making  the  transition  to  the  test  waterline.  The  transition  is 
accomplished  using  eight  separate  parts  for  each  side  of  the  hull  and  one  part  for  the 
stern.  While  this  could  be  accomplished  using  fewer  parts,  using  separate  parts  allows 
for  more  control  over  the  final  shape  and  orthogonality  of  the  mesh. 

These  parts  use  a  curve  developed  by  first  extruding  a  single  seven  inch  thick 
nonnal  mesh  layer  from  the  hull  up  to  and  beyond  the  test  waterline.  A  surface  is  fonned 
on  the  outer  boundary  of  this  part  by  creating  a  face  set  out  of  the  boundary.  A  curve  can 
then  be  created  to  extend  the  inner  liner  to  the  desired  height  at  the  outer  boundary  by 
using  the  lp3  feature  to  pick  points  along  this  surface.  It  is  then  necessary  to  manually 
change  the  coordinates  of  the  curve  extending  it  to  the  desired  height,  while  at  the  same 
time  keeping  it  on  the  surface.  At  this  point  it  is  necessary  to  exit  out  of  TrueGrid  and 
save  the  file.  The  file  can  then  be  opened  with  any  text  editor,  the  part  and  surface  used 
to  generate  this  curve  removed,  keeping  the  curve  generated  with  them.  The  file  is  then 
saved  and  reloaded.  Since  the  wetted  hull  is  symmetric,  this  curve  can  be  transformed 
over  to  the  opposite  side  by  using  the  rzx  transformation  command  which  reflects  the 
curve  about  the  xz-axis. 

The  upper  portion  of  the  inner  liner  is  then  constructed  by  attaching  the  outer  top 

edge  of  the  individual  parts  to  this  curve.  The  structural  mesh  line  transitions  from  above 

to  below  the  test  waterline  at  three  locations  along  each  side  of  the  hull.  This  requires 

constructing  three  tetrahedron  elements  on  each  side  of  the  hull.  While  in  most  instances 

these  would  not  be  desirable,  they  are  required  to  achieve  a  proper  waterline  and  to  allow 

for  construction  of  the  outer  fluid  mesh  at  the  waterline.  These  elements  can  be  deleted 

later  after  the  file  is  output  into  LS-DYNA.  Their  location  near  the  waterline  make  their 
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deletion  transparent  to  the  simulation  results.  The  parts  used  to  make  the  aft  section  of 
the  inner  mesh  were  made  primarily  in  the  block  part  phase,  while  the  parts  having  more 
than  two  elements  were  constructed  using  the  blude  part  phase.  Figure  114  shows  the  aft 
section  of  the  upper  inner  liner. 


Figure  114.  Transition  to  Test  Waterline:  Aft  Section  of  Upper  Inner  Liner 

The  forward  three  quarters  of  the  upper  mesh  consists  of  two  parts.  An  initial  part 
is  first  extruded  from  hull  to  fonn  the  majority  of  the  fluid  mesh  in  this  region.  A  single 
element  thick  wedge  shaped  part  is  then  extruded  from  the  top  of  this  fluid  mesh.  This 
transitions  the  structural  mesh  line,  which  is  at  a  constant  height  of  254  inches  for  the 
majority  of  this  region,  to  the  test  waterline  of  258  inches.  These  parts  are  shown  in 
Figure  115. 

The  inner  mesh  layer  is  completed  by  extruding  a  fluid  mesh  layer  from  the  upper 
stern  region,  transitioning  it  down  towards  the  test  waterline.  The  stem  seam  is  then 
transitioned  using  a  wedge  shaped  stem  part  in  the  same  manner  in  which  the  keel  and 
bow  seams  were  transitioned.  Figure  116  shows  the  completed  inner  layer  of  the  mesh. 
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Figure  115.  Transition  to  Test  Waterline:  Forward  Section  of  Upper  Inner  Fluid 

Mesh 


Figure  116. 


Completed  Inner  Mesh  Layer  of  Inner  Liner 
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Upon  completion  of  the  inner  mesh  layer,  the  fluid  mesh  has  been  transitioned  to 
the  test  waterline  everywhere  except  in  the  stern  region.  This  is  due  to  the  large 
transition  distance  from  the  structural  mesh  line  to  the  test  waterline  here.  To  avoid 
extremely  non-orthogonal  inner  mesh  elements,  this  transition  is  accomplished  using  two 
layers  of  fluid  mesh. 

The  second  layer  of  inner  liner  is  eight  inches  and  thickness  and  consists  of  five 
parts;  a  port  and  starboard  mesh,  a  sonar  dome  mesh,  a  stern  mesh,  and  a  bow  mesh.  The 
sonar  dome  is  first  extruded  from  the  inner  mesh  sonar  dome  constructed  previously. 
The  port  and  starboard  meshes  are  then  extruded  from  their  respective  inner  meshes  and 
their  respective  half  of  the  inner  mesh  keel  seam  and  the  entire  sonar  dome  seam  on  their 
respective  side.  In  addition,  the  side  of  the  inner  mesh  stem  seam  is  encompassed  in  this 
extrusion.  Curves  must  be  generated  at  the  desired  waterline  to  which  the  upper  outer 
edge  of  the  port  and  starboard  meshes  can  be  attached.  These  parts  merge  together  with 
each  other  and  the  sonar  dome  outer  layer  to  form  a  continuous  outer  layer  of  the  inner 
mesh.  A  single  mesh  layer  is  extruded  from  the  stem,  transitioning  the  upper  boundary  to 
the  waterline.  The  final  part  is  the  bow  constmcted  in  the  block  part  phase.  This  part 
extends  the  bow  out  eight  inches  and  connects  the  port  and  starboard  mesh  layers  at  the 
bow.  The  completed  outer  layer  of  the  inner  liner  is  shown  in  Figure  117. 
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Figure  117.  Completed  Outer  Mesh  Layer  of  Inner  Liner. 

C.  INNER  MESH 

After  completion  of  the  inner  liner,  the  fluid  mesh  construction  process  moves  to 
extending  the  fluid  out  to  the  desired  depth.  This  process  is  accomplished  using  multiple 
fluid  mesh  layers.  The  first  layer  has  been  termed  the  “Inner  Mesh”.  The  purpose  of  this 
mesh  is  to  gradually  “smooth”  the  contours  of  the  ship.  This  greatly  simplifies 
construction  of  the  outer  fluid  mesh  layers.  This  inner  mesh  extends  out  a  nominal  five 
feet  from  the  outer  boundary  of  the  inner  liner. 

Before  beginning,  it  is  highly  advised  to  use  the  labels  cracks  feature  to  identify 
any  improperly  merged  nodes.  If  these  cracks  are  allowed  to  remain,  they  will  only 
multiply  as  the  fluid  model  is  expanded  outward.  These  cracks  should  be  removed  either 
by  reconstruction  of  the  parts  in  which  they  appear  or  by  using  the  pn  command. 
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The  mesh  extrusion  process  used  in  the  inner  liner,  though  challenging  for  the 
beginner,  was  relatively  simple.  Due  to  the  contours  of  the  hull,  constructing  the  inner 
mesh  requires  a  much  more  complicated  extrusion  process.  To  accomplish  this,  a  good 
conceptual  understanding  of  how  the  extrusion  process  works  is  necessary.  The  mesh 
extrusion  process  involves  first  projecting  the  mesh  being  constructed  to  the  predefined 
outer  surface  contour.  The  mesh  is  then  projected  to  the  inner  surface  through  which  the 
mesh  is  extruded.  This  method  presents  a  problem  in  these  areas  as  “shadow  regions” 
exist.  These  shadow  regions  can  be  illustrated  with  a  flashlight  shining  onto  the  model  as 
shown  in  Figure  118. 


Figure  118.  Fluid  Model  Projection  Illustration  Showing  “Shadow  Regions”  Created 

by  the  Contours  of  the  Ship  Model 


The  flashlight  can  be  thought  of  as  a  desired  final  outer  boundary  position  of  a 
fluid  mesh  line  which  begins  on  the  previously  defined  surface.  From  this  illustration, 
one  can  see  that  by  simply  moving  the  flashlight  to  a  different  location  these  shadow 
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regions  would  be  eliminated.  While  this  is  true,  the  resulting  mesh  would  be  highly  non- 
uniform  with  poor  aspect  ratios  due  to  the  large  separations  created  in  these  transition 
areas.  In  these  areas  it  is  desirable  to  transition  a  mesh  line  which  begins  in  the  shadow 
region  to  the  location  of  the  flashlight  as  shown  in  Figure  118. 

These  shadow  regions  can  be  reduced  or  avoided  entirely  by  projecting  from  all 
three  sides  of  the  hull.  The  best  way  to  accomplish  this  is  by  using  separate  parts  for 
each  side  in  conjunction  with  curves  and  internal  projections.  These  curves  and  internal 
projections  can  be  thought  of  as  essentially  using  fiber-optic  cables  to  guide  the  projected 
mesh  lines  from  the  inner  mesh  to  the  desired  location  on  the  outer  boundary.  The  end 
result  is  a  high  quality  mesh  with  a  minimal  amount  of  non-uniformity. 

The  inner  mesh  was  constructed  using  five  separate  parts;  a  port  side  mesh,  a 
starboard  side  mesh,  a  sonar  dome  mesh,  a  keel  mesh,  and  a  small  lower  bow  plug  mesh. 
These  parts  are  all  created  in  the  blude  part  phase.  Before  beginning  this  extrusion 
process,  it  is  first  necessary  to  define  what  the  outer  boundary  of  this  mesh  will  be  when 
completed.  It  is  desirable  to  have  a  convex  surface,  but,  at  the  same  time,  it  is  also 
desirable  to  not  have  any  large  variations  in  the  inner  mesh  spacing  nonnal  to  the  hull. 
This  requires  an  outer  boundary  which  is  not  perfectly  cylindrical.  To  form  this  surface  a 
number  of  cubic  spline  curves  are  formed  using  the  csp3  command.  These  curves  are 
constructed  so  that  they  follow  the  contours  of  the  hull  at  a  nominal  distance  of  five  feet 
from  the  inner  liner.  These  curves  essentially  trace  the  outline  of  the  desired  outer 
boundary  of  this  inner  mesh.  These  curves  are  then  combined  to  form  surfaces  using  the 
blend4  function  in  TrueGrid  and  the  surfaces  combined  to  form  a  single  surface  using  the 
sds  command.  Figure  119  shows  the  curves  and  surfaces  used  to  form  the  outer 
boundary  of  the  inner  mesh. 
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Figure  119.  Curves  and  Surfaces  Used  to  Form  Outer  Boundary  of  Inner  Mesh 


In  addition  to  these  curves  and  surfaces,  several  additional  curves  and  surfaces  are 
created  to  serve  in  order  to  control  the  mesh.  Controlling  the  mesh  in  this  manner  is 
necessary  to  create  a  uniform  mesh  at  the  outer  boundary. 

The  construction  of  the  port  and  starboard  inner  mesh  is  very  complex.  The 
commands  used  to  generate  the  starboard  inner  mesh  are  shown  below: 


blude  4  stbdfc2  1  2  3;1  %inner;1  2  3;-80  4897  5600;-450  -390;-75  176 
258; 

pb  1  2  1  1  2  1  xyz  -69.286003  0.0  210.75 
pb  1  2  2  1  2  2  xyz  -77.0  -212.2  258.0 
pb  1  1  1  1  1  1  xyz -80.0  -1.1 37e-07  130.0 
pb  3  2  3  3  2  3  xyz  5595.0703  -15.323369  258.0 
pb  3  2  1  3  2  1  xyz  5589.0005  -7.280057  -3.0631487 
pbs  3  113  11  xyz  806.1 

pb  2  2  1  2  2  1  xyz  4897.7222  -3.4121273  -13.485395 

pb  2  2  3  2  2  3  xyz  4894.4614  -173.06313  257.99969 

pbs  2  112  11  xyz  812.21 

pbs  3  13  3  13  xyz  89.1 

pbs  112  112  xyz  12.1 

mseq  i  200  50 

mseq  k  40  0 

insprt  16  12 

pb  3  2  2  3  2  2  xyz  5594.8252  -10.205372  6.3662529 
pbs  3  1  2  3  1  2  xyz  9.1 
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insprt  12  12 

pbs  2  1  2  2  1  2  xyz  8.2 

pbs  2  112  11  xyz  812.1 

pb  2  2  1  2  2  1  xyz  -56.667458  -6.0664253  210.01097 

pb  2  2  2  2  2  2  xyz  -57.072742  -12.096737  210.1 1931 

pb  2  2  3  2  2  3  xyz  -61 .072933  -21 3.76886  257.99969 

pb  1  2  2  1  2  2  xyz  -69.286003  0.0  210.75 

pbs  2  1  3  2  1  3  xyz  808.394.91 

pb  1  1  2  1  1  2  xyz  -80.0  -1.1 37E-07  1 30.0 

pbs  4  2  34  2  3  xyz  1.1319.2 

pbs  3  2  3  3  2  3  xyz  1 .9.1 

insprt  1  2  2  35 

insprt  12  3  2 

dei  1  3;  1  2;  3  4; 

pbs  4  13  4  13  xyz  1 15.1 

pbs  3  2  3  3  2  4  xyz  1 .23.3 

pbs  4  2  4  4  2  4  xyz  1.78.2 

pbs  4  2  34  2  3  xyz  1.78.1 

pb  4  1  4  4  1  4  xyz  863.36  -382.3  [%dftaft+940*%trim/5672] 

pb  3  1  3  3  1  4  xyz  8.16E+02  -3.79381 07E+02  [%dftaft+893*%trim/5672] 

pbs  3  2  1  3  2  1  xyz  1.120.2 

pbs  4  2  2  4  2  2  xyz  1 .642.3 

pbs  4  2  1  4  2  1  xyz  1.1300.1 

pbs  3  2  2  3  2  2  xyz  1.50.2 

pbs  3  113  11  xyz  802.212.18 

pbs  4  114  11  xyz  802.212.13 

pbs  5  1  4  5  1  4  xyz  800.473.106 

pb  5  1  3  5  1  3  xyz  4.8934E+03  -2.4791 3E+02  [%dftaft+%trim*4970/5672- 
16] 

pbs  6  13  6  13  xyz  1 15.99 

pb  4  1  2  4  1  2  xyz  8.5765637E+02  -30  -3.3454285E+01 

pb3  1  23  1  2  xyz  815 -30  -21.867 

mseq  i  0  5  1  5-10 

mseq  k  1  0  0 

edge  3  2  4  4  2  4  1.2 

edge  4  2  4  5  2  4  1 .2 

edge  52  4  6  24  1.2 

edge  1  2  3  2  23  1.2 

edge  2  2  3  3  2  3  1.2 

curs  3  1  4  6  1  4  801 

curs  1  1  3  3  1  3  801 

curs  4  13  6  13  115 

curs  4  2  3  6  2  3  112 

curs  111112  12 

curs  112  113  12 

curs  2  12  5  12  8 
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curs  2  1  1  5  1  1  812 
edge  12  12  2  1  1.3 
edge  2  2  1  3  2  1  1.3 
edge  3  2  1  4  2  1  1.3 
edge  4  2  1  52  1  1.3 
edge  5  2  1  6  2  1  1.3 
curs  5  1  1  6  1  1  806 
curs  2  2  2  6  2  2  5 
curs  5  12  6  12  9 
curs  6  1  3  6  1  4  89 
curs  6  1  2  6  1  3  89 
edge  6  2  1  6  22  1.1 
edge  6  2  2  6  2  3  1.1 
edge  6  2  3  6  2  4  1.1 
edge  1  2  2  1  23  1.3 
curs  6  1  4  6  2  4  40 
curs  6  1  1  6  2  1  94 
curs  6  1  2  6  2  2  43 
curs  5  1  1  5  2  1  98 
curs  5  1  2  5  2  2  45 
curs  2  1  1  2  2  1  500 
insprt  12  2  9 

pbs  3  2  3  3  2  3  xyz  1 .491 .3 

pbs  3  2  2  3  2  2  xyz  1 . 1 96.2 

pbs  3  2  1  3  2  1  xyz  1.1989.1 

pbs  3  13  3  13  xyz  14.6 

pbs  3  1  1  3  1  1  xyz  14.127 

pbs  3  1  2  3  1  2  xyz  14.124 

curs  3  1  1  3  1  3  14 

insprt  12  3  5 

pbs  4  2  3  4  2  3  xyz  1.64.4 

pbs  4  2  2  4  2  2  xyz  1 .734.3 

pbs  4  2  1  4  2  1  xyz  1.734.2 

pbs  4  114  11  xyz  17.95 

pbs  4  13  4  13  xyz  17.1 

pbs  4  1  2  4  1  2  xyz  17.98 

curs  4  1  14  13  17 

sfi  7  8;1  2;-1;sd  19 

sfi  7  8;1  2;-2;sd  15 

sfi  2  7;1  2;-1;sd  31 

sfi  2  7;1  2;-2;sd  35 

sfi  5  8;  1  2;-4;sd  999 

sfi  1  5;1  2;-3;sd  999 

sfi  1  8;-1;1  4;sd  803 

sfi  1  8;  -2;  1  4;sd  1 

res  1  1  1  5  2  4  j  1 
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as  5  1  1  8  2  2  j  1  9 

das  7  1  2  8  2  4  j  [150/(%inner+2)]  10 

tmei  1  5;1  2;2  3;50  0  1 

tmei  5  8;  1  2;1  4;50  0  1 

relaxi  7  8;-2;1  4;50  0  1 

relaxi  2  7;-2;2  3;50  0  1 

relaxi  7  8;-1  ;1  4;50  0  1 

relaxi  2  7;-1;2  3;50  0  1 

bb  7  1  1  8  2  1  17; 

bb  2  1  1  7  2  1  21; 

endpart 

The  port  mesh  is  constructed  in  a  similar  fashion.  Key  points  to  highlight  in  this 
example  are  the  use  of  the  tme  and  relax  interpolation  methods  and  the  use  of  parameter 
functions.  The  mesh  must  be  interpolated  properly  to  generate  a  high  quality  mesh.  This 
requires  a  fair  amount  of  practice  as  some  interpolation  methods  work  better  than  others 
depending  on  the  topology  of  the  mesh  and  the  way  in  which  the  boundaries  are 
controlled.  The  parameter  functions  used  are  prefaced  by  the  percent  sign.  These 
represent  previously  defined  parameters  created  using  the  para  command.  Using  these 
parameters  and  some  simple  FORTRAN  expressions  in  the  model  is  key  to  allowing 
automation  of  the  mesh  generation  process.  Using  some  creative  modeling  techniques, 
major  modifications  can  be  made  to  the  fluid  model  by  changing  just  a  few  of  these 
parameter  values. 

The  sonar  dome  and  keel  were  constructed  separately  to  allow  for  greater  control 
of  the  mesh  in  these  challenging  areas.  To  avoid  having  several  cracks  between  these 
individual  parts,  the  bb  command  was  used  extensively.  This  command  allows  a  face  of 
a  part  to  be  defined  as  a  master  side  in  the  part  phase.  Subsequent  parts  can  define  their 
coinciding  face  using  the  same  boundary  number,  creating  a  slave  side.  When  the  parts 
are  merged  together,  slave  and  master  sides  are  glued  together  with  their  nodes  properly 
merged  throughout.  Without  using  this  feature,  construction  of  the  inner  mesh  in  this 
manner  would  not  be  feasible.  Figure  120  shows  the  parts  used  to  generate  this  inner 
mesh  after  being  merged  together. 
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Figure  120.  Completed  Inner  Mesh 
D.  TRANSITIONAL  MESH 

The  inner  mesh  is  encapsulated  by  another  fluid  mesh  which  extends  the  fluid 
mesh  out  another  five  feet,  continuing  the  transition  to  a  cylindrical  outer  boundary.  This 
mesh  consists  of  three  parts;  the  mesh  encapsulating  the  inner  mesh  and  a  bow  and  stern 
cone. 

To  extrude  this  mesh  from  the  inner  mesh,  the  face  elements  of  the  outer 

boundary  of  the  inner  mesh  must  be  selected  and  saved  as  a  surface.  In  the  previous  parts 

of  the  mesh  it  was  necessary  to  perform  this  process  manually  using  the  lasso.  A  key 
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aspect  in  automating  the  mesh  generation  process  is  to  have  this  process  performed 
automatically.  This  is  done  by  using  the  fset  surface  command  in  combination  with  the 
surface  used  to  initially  form  the  outer  surface  of  the  inner  mesh.  This  command  selects 
all  of  the  faces  in  the  model  which  have  a  user  specified  number  of  nodes  within  a  user 
defined  distance  from  the  surface. 

fset  innrmesh  =  surface  803  2  3 

By  creating  surfaces  between  curves  on  each  end  of  the  inner  mesh,  bow  and  stern 
face  sets  can  be  created  in  a  similar  fashion.  It  should  be  noted  that  the  parts  from  which 
the  face  sets  are  to  be  defined  must  be  displayed  for  this  option  to  work. 

The  transitional  mesh  is  much  simpler  to  extrude  then  the  inner  mesh  since  the 
contours  of  the  inner  mesh  are  convex  in  nature.  To  form  the  outer  surface  of  this  mesh, 
curves  are  used  to  form  a  surface  very  similar  to  that  which  was  done  for  the  inner  mesh. 
The  difference  in  constructing  this  mesh  is  that  the  surface  must  have  a  bow  and  stern 
portion.  To  simplify  construction  of  the  outer  mesh  which  will  eventually  encapsulate 
this  transitional  mesh,  this  surface  is  made  to  be  spherical  at  each  end.  The  curves  and 
surfaces  used  to  form  this  outer  boundary  are  shown  in  Figure  121. 


Figure  121. 


Curves  and  Surfaces  Used  to  Form  Outer  Boundary  of  Transitional  Mesh 
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One  notable  aspect  of  the  transitional  mesh  is  its  ability  to  transition  the  inner 
mesh,  which  has  a  set  waterline  of  258  inches,  to  a  different  waterline  which  would  result 
from  different  loading  conditions  of  the  ship.  Though  not  used  in  the  simulations  in  this 
paper,  the  ability  to  model  different  loading  conditions  is  designed  for  future  study.  This 
is  accomplished  by  first  defining  a  large  planar  surface  between  two  linear  curves  defined 
using  parameters  inside  of  functions  that  extrapolate  the  defined  waterline  out  several 
thousand  inches  from  the  structure.  A  series  of  linear  curves  are  generated  at  a  height  of 
258  inches  and  equal  to  the  transitional  mesh  in  width  at  each  location  where  an  outer 
waterline  node  of  the  transitional  mesh  is  attached.  These  curves  are  projected  to  the 
waterline  surface  using  the  projcur  command.  The  outer  waterline  nodes  of  the 
transitional  mesh  are  then  connected  to  these  floating  curve  points  using  the  pbs 
command.  The  completed  transitional  mesh  showing  the  three  individual  parts  merged 
together  is  shown  in  Figure  122. 


E.  OUTER  MESH 

The  outer  mesh  consists  of  a  single  mesh  having  a  prolate  spheroid  shaped  outer 
boundary  contour.  This  portion  of  the  fluid  mesh  is  one  of  the  easiest  portions  to  extrude 
due  to  the  convex  contours  of  the  transitional  mesh.  The  most  complex  aspect  of  this 
mesh  is  in  defining  the  surfaces  and  curves  which  form  the  outer  boundary. 

This  outer  mesh  is  designed  to  be  easily  expanded  or  reduced  in  size  by  using 
parameters.  This  requires  that  curves  and  surfaces  be  defined  using  these  parameters. 
The  outer  boundary  of  the  mesh  is  then  attached  to  these  floating  surface  and  curve  points 
rather  than  fixed  coordinates.  This  allows  the  outer  mesh  to  move  with  these  curves  and 
surfaces  which  are  modified  by  only  a  few  parameter  inputs.  Some  care  should  be  used 
when  connecting  to  surfaces  and  curves  as  the  numbering  scheme  of  interpolated  curves 
and  surfaces  changes  as  the  size  changes.  Spheres  and  circular  curves  maintain  their 
numbering  scheme  regardless  of  size  making  there  use  preferable  if  possible.  The 
completed  outermesh  is  shown  in  Figure  123. 
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pbs  2  2  1  2  2  1  xyz  25.1.15 
pbs  2  3  1  2  3  1  xyz  25.21.55 
pbs  2  1  1  2  1  1  xyz  25.89.55 
pbs  2  1  2  2  1  2  xyz  44.1 
pbs  2  3  2  2  3  2  xyz  44.353 
pbs  2  2  2  2  2  2  xyz  44.177 
curs  2  1  2  2  3  2  44 
curs  2  1  1  23  1  127 
insprt  12  11 
pbs  2  1  2  2  1  2  xyz  43.1 
pbs  2  3  2  2  3  2  xyz  43.275 
pbs  2  2  2  2  2  2  xyz  43.138 
pbs  2  1  1  2  1  1  xyz  141.240 
pbs  2  3  1  2  3  1  xyz  141.1 
pbs  2  2  1  2  2  1  xyz  129.55.13 
curs  2  1  1  23  1  141 
curs  2  1  2  2  3  2  43 
insprt  12  2  9 

pb  3  3  2  3  3  2  xyz  143.01845  417.35803  258.0 

pb  3  1  2  3  1  2  xyz  143.01845  -417.35803  258.0 

pb  3  2  2  3  2  2  xyz  1 .341 5643E+02  -5.380081 5E-01  -1 .5499544E+01 

curs  3  1  1  3  3  1  241 

insprt  114  5 

pb  4  2  2  4  2  2  xyz  5461 .3589  -0.7900944  -1 1 9.76952 

pb  4  3  2  4  3  2  xyz  5453.7617  281 .53503  257.97253 

pb  4  1  2  4  1  2  xyz  5453.7622  -281 .53503  257.97256 

curs  4  1  1  4  3  1  49 

curs  6  2  2  6  3  2  41 

curs  6  1  2  6  2  2  42 

mseq  i  50  2  20  10  0 

curs  1  2  2  6  2  2  40 

curs  1  2  1  6  2  1  301 

curs  2  3  1  5  3  1  401 

curs  2  1  1  5  1  1  402 

res  1  1  1  6  3  2  k  %ores 

sfi  1  6;  1  3;  -2;sd  30 

sfi  -1  0  -6;  1  3;  1  2;sd  999 

sfi  1  6;  -1  0  -3;  1  2;sd  999 

Sfi  1  6;1  3;-1  ;SD  1 

tmei  1  6;  1  3;  1  2;  50  0  1 

relaxi  1  6;1  3;-2;20  .5  1 

relaxi  1  6;1  3;-1;20  1  1 

endpart 
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Figure  123.  Completed  Outer  Mesh 
F.  DAA  BOUNDARY 

The  orpt  command  should  be  used  prior  to  defining  the  DAA  boundary.  This 
command  ensures  that  the  face  sets  are  defined  with  their  normal  vectors  pointing  away 
from  the  point  specified  by  the  orpt  command.  Using  a  point  located  in  the  center  of  the 
fluid  mesh  ensures  that  this  occurs  at  all  locations  on  the  boundary. 

The  DAA  boundary  is  defined  by  defining  a  face  set  on  the  submerged  outer 
boundary  of  the  outer  mesh  using  the  surface  generated  to  form  the  outer  boundary  of  the 
mesh  and  the  fset  command.  Once  defined,  all  other  face  sets  should  be  deleted  using  the 
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delset  command.  This  avoids  having  these  sets  output  into  LS-DYNA.  The  fluid  model 
is  output  into  LS-DYNA  by  using  the  following  command: 

Isdyna  keyword  write 

The  completed  coupled  fluid  structure  model  is  shown  in  Figure  124. 


Figure  124.  Coupled  Fluid-Structural  Model 
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APPENDIX  C :  TRUEGRID  MODELING  COMMANDS 


This  appendix  contains  a  brief  description  of  many  of  the  commands  used  in  the 
fluid  mesh  modeling  process  in  TrueGrid  outlined  in  Appendix  B.  A  complete 
description  of  all  of  the  commands  used  in  the  fluid  mesh  modeling  process  plus  many 
more  that  were  not  used  can  be  found  in  Reference  [5],  Commands  have  been  grouped 
according  to  what  phase  of  the  modeling  process  they  are  used  in. 


IMPORT/EXPORT  COMMANDS 

1 .  READMESH:  This  command  is  utilized  to  import  a  NASTRAN,  LS-DYNA, 
DYNA3D,  or  NEUTRAL  mesh  into  TrueGrid.  When  used  in  conjunction 
with  the  EXCLUDE  command  the  nodes  of  parts  generated  in  TrueGrid 
which  have  merged  with  nodes  on  the  imported  mesh  will  not  be  redefined. 
Material  properties  of  imported  NASTRAN  files  must  be  redefined  if  being 
exported  in  another  format. 

2.  DELMATS/DELSPDS:  These  commands  delete  the  material  and  spring 
properties  of  the  imported  model.  TrueGrid  does  not  support  conversion  of 
material  and  spring  properties  from  NASTRAN  to  LS-DYNA.  These 
properties  must  be  deleted  to  allow  the  model  to  be  exported  in  LS-DYNA 
format. 

3.  LSDYNA  KEYWORD  WRITE:  This  exports  the  model  in  LS-DYNA 
keyword  fonnat.  The  exported  file  is  saved  in  the  working  directory  as  a 
trugrdo  file. 

PARAMETERS  &  PROGRAMMING 

4.  PARA:  Parameters  can  be  given  a  numerical  value  using  this  command,.  To 
use  this  command  a  parameter  is  named  and  then  followed  by  its  numerical 
value.  The  parameters  are  used  in  the  file  by  preceding  the  name  of  the 
parameter  with  a  percent  sign. 
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PARA  daafwd  900  daaaft  954  daachg  900; 

curd  37  arc3  seqnc  rt  5449  %daafwd  258  rt  5449  0  [258-%daafwd] 
rt  5449  [0-%daafwd]  258; 

Parameters  can  allow  automation  of  the  mesh  generation  process  when  used  in 
combination  with  fonnulas  and  basic  FORTRAN-like  code. 

5.  IF/ELSEIF/ELSE/ENDIF:  TrueGrid  allows  basic  control  statements  much 
like  that  of  FORTRAN.  This  allows  automation  of  the  mesh  generation 
process  when  used  in  conjunction  with  parameters. 

if  (%outerk>20)  then 
para  ores  1 ; 
else 

para  ores  0.96 
endif 

CURVE  COMMANDS 

Note:  All  curve  commands  begin  with  the  CURD  command. 

6.  CPCD:  This  command  allows  the  creation  of  a  curve  by  perfonning  a 
transformation  to  a  previously  defined  curve.  Since  the  wetted  surface 
structure  is  symmetric  about  the  xz-plane,  many  curves  were  reflected  about 
this  plane  using  the  RZX  command  in  conjunction  with  the  CPCD  command. 
Translation,  reflection,  and  rotation  transformations  in  all  three  dimensions 
are  possible  using  this  command. 

7.  PROJCUR:  This  command  projects  a  curve  onto  a  previously  defined 
surface.  This  is  a  very  useful  command  used  alone  or  in  conjunction  with  the 
other  curve  commands,. 

8.  CSP3:  This  command  allows  the  creation  of  cubic  spline  curves.  Extremely 
useful  for  defining  the  outer  boundary  curves  as  well  as  the  mesh  control 
curves  in  the  sonar  dome  and  skeg  regions.  This  command  allows  the  user  to 
choose  a  natural  end  derivative,  specify  an  end  derivative,  or  match  an 
existing  curve’s  end  derivative. 
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9.  ARC3:  This  command  allows  the  creation  of  a  curve  through  a  circular  arc 
defined  by  three  points.  This  is  a  useful  command  for  defining  spherical 
contours  of  the  outer  mesh. 

10.  LP3:  The  most  basic  curve  command.  This  command  defines  a  curve  by  a 
series  of  coordinates  usually  selected  by  picking  a  node  or  surface  point  with 
the  computer  mouse.  This  command  allows  modification  of  the  coordinates 
which  is  particularly  useful  when  manually  creating  the  initial  curves  for  the 
inner  liner. 

11.  *PT:  Following  the  above  commands  by  PT  allows  the  coordinates  to  be 
replaced  with  surface  or  curve  points.  This  is  particularly  useful  when 
automating  the  mesh  generation  process. 


SURFACE  COMMANDS 

Note:  All  surface  commands  are  preceded  by  the  SD  command. 

12.  PLAN/SP/CY:  These  commands  form  an  infinite  planar,  spherical,  and 
infinite  cylindrical  surface  respectively. 

13.  BLEND3/4:  This  command  is  used  to  blend  three  or  four  curves  together  to 
form  a  surface.  Very  useful  to  form  the  contoured  surfaces  of  the  inner  and 
transitional  meshes  as  well  as  the  surfaces  required  to  form  many  of  the  face 
sets.  Care  should  be  used  as  the  curves  are  combined  in  the  order  in  which 
their  nodes  are  defined.  If  a  curve  defined  left  to  right  is  combined  with  a 
curve  defined  right  to  left  the  resulting  surface  will  cross  over  itself. 

14.  RULE3D:  This  command  combines  a  ruled  surface  between  two  three 
dimensional  curves.  The  same  rule  as  in  the  BLEND  case  applies  concerning 
the  order  of  the  curve  nodes. 

15.  SDS:  Combines  multiple  surfaces  into  one  surface. 

16.  FACESET:  Creates  a  surface  out  of  a  previously  defined  face  set. 

17.  R3DC:  Creates  a  surface  by  rotating  a  3D  curve  a  specified  angle  about  an 
arbitrary  axis. 
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PART  COMMANDS 


18.  BLOCK:  This  command  is  the  most  basic  part  generating  process  in 
TrueGrid.  This  command  allows  the  user  to  generate  parts  in  the  part  phase 
by  defining  nodal  coordinates  (cartesian)  and  through  projections  and 
curve/edge  attachments. 

19.  BLUDE:  This  command  extrudes  a  set  of  polygons  through  a  defined  face  set 
by  following  the  mesh  lines  formed  by  the  block  part.  This  command  requires 
a  thorough  understanding  of  the  block  command  from  which  it  is  based.  The 
block  part  must  completely  cover  the  face  from  which  it  is  being  extruded. 
This  command  is  essential  in  developing  a  properly  coupled  fluid-structure 
model. 

20.  PB/PBS:  The  PB  command  attaches  the  vertices  of  a  block  part  to  a  set  of 
coordinates.  The  PBS  command  attaches  the  vertices  of  a  block  part  to 
surface  or  curve  points.  Careful  use  of  the  PBS  command  is  necessary  to 
automate  the  mesh  generation  process. 

21.  CURS:  This  command  attaches  an  edge  of  the  block  part  to  a  previously 
defined  curve.  The  end  nodes  of  the  block  part  edge  are  attached  to  the  curve 
first  with  the  remaining  nodes  evenly  distributed  along  the  curve. 

22.  EDGE:  Similar  to  the  CURS  command  only  this  command  attaches  an  edge 
of  the  block  part  to  an  edge  of  a  surface. 

23.  SFI:  This  command  projects  the  face  of  a  block  part  to  a  previously  defined 
surface.  All  surrounding  vertex  and  edge  attachments  should  be  made  prior  to 
using  this  command. 

24.  DEI:  This  command  deletes  a  region  of  the  block  part. 

25.  RES/AS/DAS:  These  commands  are  used  to  distribute  the  nodes  along  the 
boundary  of  the  block  part.  The  first  or  last  spacing,  both,  or  the  relative 
spacing  in  a  given  direction  can  be  specified. 

26.  BB:  This  command  defines  a  boundary  interface  which  can  be  used  to  attach 
two  separate  parts  or  to  apply  a  transformation  to  a  face  set  of  the  same  block 
part. 
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27.  NORMAL:  Used  in  conjunction  with  the  BB  command,  this  command 
translates  the  slave  boundary  face  of  a  block  part  out  a  specified  distance 
nonnal  to  the  master  boundary  face. 

28.  INSPRT:  This  command  is  used  to  insert  a  partition  in  the  block  part  being 
created.  This  is  done  primarily  in  order  to  add  additional  boundary  constraints 
to  the  block  part 

29.  MSEQ:  This  command  is  used  to  add  elements  to  or  remove  elements  from 
the  block  part  being  created.  Often  it  is  easier  to  begin  the  construction 
process  having  only  a  few  elements  as  a  large  amount  of  elements  “clutters” 
the  picture  making  attachments  more  difficult  to  see.  This  command  can  be 
used  to  add  the  desired  number  of  elements  to  the  block  part  just  prior  to 
attaching  an  edge  to  a  curve  or  surface  edge,  adding  a  partition,  and/or 
projecting  to  a  surface. 

30.  TMEI/RELAXI:  These  are  commands  used  to  create  a  uniform  mesh  through 
interpolation  methods.  These  commands  should  be  issued  as  the  last 
commands  in  the  block  part  construction.  Failure  to  use  these  or  other 
interpolation  methods  will  result  in  an  improper  extrusion  and/or  a  low  quality 
mesh. 

31.  ENDPART:  This  commands  ends  the  part  being  worked  on.  No  further  part 
construction  can  occur  after  this  command  has  been  issued. 

32.  MERGE:  This  command  ends  the  part  phase  and  merges  the  individual  parts 
together. 


FACE  SET  COMMANDS 

33.  FSET  “...”  =  LB6:  This  command  allows  a  face  set  to  be  defined  by 
manually  selecting  a  group  of  face  elements  on  a  part  while  in  the  merge 
phase  using  the  “pick”  -  “sets”  -  “faces”  option.  The  faces  are  then  selected 
using  a  lasso  to  circle  the  outer  nodes  of  the  desired  faces. 

34.  FSET  “...”  =  SURFACE:  This  command  defines  a  face  set  by  selecting  all  of 
the  faces  with  a  user  defined  number  of  nodes  within  a  defined  distance  from 
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a  surface.  This  command  is  necessary  to  automate  mesh  generation  processes 
involving  mesh  extrusion. 

35.  FSET  “...”  =  SET  “. . This  command  defines  a  face  set  as  another  face  set. 

36.  FSET  OR  SET  This  command  adds  the  face  set  elements  from 

one  face  set  to  another  face  set. 

37.  ORPT:  This  command  allows  the  vector  normals  of  the  face  set  to  be 
controlled.  For  USA,  the  vector  nonnals  should  be  pointing  outward. 

38.  DELSET  This  command  deletes  a  face  set.  Prior  to  exporting  the 

model  into  LS-DYNA,  all  of  the  face  sets  defined  should  be  deleted  with  the 
exception  of  the  face  set  defining  the  DAA  boundary.  This  eliminates 
unnecessary  data  being  exported  into  the  output  file. 


MERGING  COMMANDS 

39.  STP:  This  command  is  used  to  set  the  merging  tolerance  between  nodes.  Care 
should  be  used  in  using  this  command,  particularly  early  in  building  the  fluid 
model,  as  it  is  only  active  in  the  merge  phase.  Using  a  smaller  tolerance  in  the 
early  stages  is  advised  to  avoid  large  cracks  appearing  later  on  in  the  building 
process.  Too  large  of  a  tolerance  at  any  stage  may  result  in  nodes  merging 
together  undesirably. 

40.  PTOL:  This  command  sets  the  merging  tolerance  for  nodes  within  a  specified 
part.  Must  be  followed  by  the  STP  command  to  become  active. 

41. BPTOL:  This  command  sets  the  merging  tolerance  for  nodes  between  two 
parts.  Must  be  followed  by  the  STP  command  to  become  active. 

42.  PN:  This  command  allows  nodes  to  be  repositioned  in  the  merge  phase. 
Useful  to  eliminate  cracks  in  the  mesh.  Having  to  use  this  command 
excessively  is  often  indicative  of  problems  in  the  construction  of  the  mesh. 

43.  LABELS  CRACKS:  This  command  is  actually  one  of  the  diagnostics  tools 
that  TrueGrid  has.  It  is  used  to  show  improperly  merged  nodes,  cracks,  in  the 
mesh  based  on  a  user  specified  maximum  allowable  angle  between  nodes. 
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44.  LABELS  L0C3D:  This  command  locates  a  solid  element  in  the  mesh.  This  is 
useful  to  find  elements  which  have  been  identified  by  either  TrueGrid  or  LS- 
DYNA  as  being  illegal  elements. 

45.  DELEM  LB/LS:  This  command  allows  blocks  or  shells  to  be  deleted  in  the 
mesh.  The  port  and  starboard  stern  tube  covers  were  deleted  using  this 
command. 

GRAPHICS  COMMANDS 

46.  POSTSCRIPT:  With  this  set  to  “on”,  the  current  frame  is  output  to  the 
working  directory  as  a  postscript  file  each  time  the  DRAW  command  is  issued 
in  the  FILL  mode.  Using  the  right  mouse  button  allows  the  user  to  specify 
resolution,  enable/disable  lighting,  control  color,  etc. 

47.  SET  PCOLOR:  This  command  defines  the  color  of  a  part  by  setting  the  levels 
of  red,  blue,  and  green. 

48.  PEXP:  This  command  allows  parts  to  be  “exploded”  a  set  distance  from  its 
actual  position.  This  is  useful  in  showing  how  individual  parts  are  constructed 
in  a  model. 
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APPENDIX  D:  USA/LS-DYNA  INPUT  DECKS 


Examples  of  the  USA  and  LS-DYNA  input  decks  are  shown  in  this  appendix. 
These  input  decks  consist  of  four  files:  the  LS-DYNA  “Keyword”  file,  and  the  USA 
FLUMAS,  AUGMAT,  and  TIMINT  files.  References  [6]  and  [7]  should  be  consulted  for 
explanation  of  each  input  shown.  Confidential  parameters  pertaining  to  the  actual  shot 
geometry  have  been  replaced  by  “XXXX”  in  each  file. 


A.  LS-DYNA  KEYWORD  FILE 


*KEYWORD 

*TITLE 

SHOT3:  COUPLED  SHIP  AND  FLUID  MODE,  DAA1  BOUNDARY 
*CONTROL_TERMINATION 
0.51,0,0,0,0 

*CONTROL_TIMESTEP 
4.0E-6, 0.9, 0,0.0, 0.0, 1,0 
*CONTROL_PARALLEL 
1,0,1 

*DEFINE_CURVE 

1 

0..4.0E-6 
0.51 ,4.0E-6 

*DATABASE_HISTORY_NODE 
*DATABASE_NODOUT 
*DATABASE_BINARY_D3PLOT 
*DATABASE_BINARY_D3THDT 
*DATABASE_EXTENT_B  I  NARY 
*BOUNDARY_USA_SURFACE 
1,1,0 

*INITIAL_DETONATION 

-1, XXXX, XXXX, XXXX, 0.0 

XXXX, XXXX, XXXX, XXXX, XXXX, 2947 

$  FLUID  NODES 

*NODE 

"NOTE:  FLUID  MODEL  NODE  COORDINATE  DEFINITIONS 
OUTPUT  FROM  TRUEGRID  ARE  ENTERED  HERE" 


$SHIP  NODES 
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"NOTE:  STRUCTURAL  MODEL  NODE  COORDINATE 
DEFINITIONS  ARE  ENTERED  HERE" 

*MAT_ACOUSTIC 

*PART 

1 

1,1,90 

*SECTION_SOLID 

1,8 

*ELEMENT_SOLID 

"NOTE:  FLUID  ELEMENTS  ARE  DEFINED  HERE  USING 
TRUEGRID  OUTPUT" 

*SET_SEGMENT 

"NOTE:  DAA  BOUNDARY  FACESET  DEFINED  IN  TRUEGRID 
IS  DEFINED  HERE" 

$  SHIP 
*PART 

*SECTION_BEAM 

*ELEMENT_BEAM 

*PART 

SPRING-ELEMENT 

*SECTION_DISCRETE 

*MAT_SPRING_ELASTIC 

*ELEMENT_DISCRETE 

*DEFINE_SD_ORIENTATION 

*ELEMENT_MASS 

*MAT_ELASTIC 

*CONSTRAINED_NODAL_RIGID_BODY 

*SET_NODE_LIST 

*PART 

*SECTION_SHELL 

*ELEMENT_SHELL 

*LOAD_SHELL_ELEMENT 

*DEFINE_COORDINATE_SYSTEM 

$  RAYLEIGH  DAMPING 

*DAMPING_GLOBAL 

*DAMPING_PART_STIFFNESS 

$  BEAM 

$SHELL 

*END 
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B.  USA  INPUT  FILES 
1.  FLUMAS 

Note:  The  number  of  nodes  (NSTRF)  and  number  of  DAA  elements  (NGENF) 
are  fluid  model  specific. 


DDG81  SHOT  3  COUPLED  ANALYSIS  -  FLUMAS 


flunam  geonam  strnam  daanam  $FLUNAMGEONAM  GRDNAM  DAANAM 


F  F  FT 
T  F  F  F 
F  FT  F 
F  F  FT 
FTF  F 
F  F  F  F 
F  F  F  F 
F  F 

0  xxxxxx  0  xxxx 
0  0  0 

0.9345E-4  60687.6 
2 

258.  0.  0.1. 

14.7  386.4 
0 
0 


$  PRTGMT  PRTTRN  PRTAMF  CALCAM 
$  EIGMAF  TWODIM  HAFMOD  QUAMOD 
$  PCHCDS  NASTAM  STOMAS  STOINV 
$  FRWTFL  FRWTGE  FRWTGR  FRESUR 
$  RENUMB  STOGMT  ROTGEO  ROTQUA 
$  PRTCOE  STRMAS  SPHERE  ROTSYM 
$  OCTMOD  CAVFLU  FRWTFV  INTCAV 
$  BOTREF  MASREF 
$  NSTRC  NSTRF  NGEN  NGENF 
$  NBRA  NCYL  NCAV 
$  RHO  CEE 
$  NVEC 

$  DEPTH  CXFS  CYFS  CZFS 
$  PATM  GRAVAC 
$  NSRADI 
$  NSORDR 


2.  AUGMAT 

Note:  The  number  of  nodes  (NSTR),  number  of  degrees  of  freedom  (NSFR),  and 
last  DAA  element  number  (JSTOP)  are  fluid  model  specific. 


DDG81  SHOT3-COUPLED  ANALYSIS-AUGMAT 


strnam  flunam  geonam  prenam  $  STRNAM  FLUNAM  GEONAM  PRENAM 


F  F  F  F 
F  F  FT 
F  F  F  F 
F  F  F  F 
11 

xxxxxx  xxxxxx  3  3 
1 

0  1  xxxx  1 


$  FRWTGE  FRWTST  FRWTFL  LUMPFM 
$  FLUSKY  DAAFRM  SYMCON  DOFTAB 
$  PRTGMT  PRTTRN  PRTSTF  PRTAUG 
$  MODTRN  STRLCL  INTWAT  CFADYN 
$  NTYPDA 

$  NSTR  NSFR  NFRE  NFTR1 
$  NSETLC 

$  NDICOS  JSTART  JSTOP  JINC 
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3. 


TIMINT 


Note:  XC,  YC,  and  ZC  define  the  charge  location.  XS,YS,  and  ZS  define  the 
location  of  the  closest  point  on  the  DAA  boundary  to  the  charge.  SLANT  is  the  distance 
from  the  charge  to  the  closest  point  on  the  DAA  boundary.  XV,  YV,  and  ZV  is  the 
coordinates  of  the  fluid  node  one  fluid  element  away  from  the  closest  structural  node  to 
the  charge.  All  of  these  values  have  been  omitted  for  confidentiality. 


DDG81  SHOT3  COUPLED  ANALYSIS-TIMINT 


prenam  posnam 
resnam 

FTF  F 
F  F  F 
1 

0.0  4.0E-6 
T  F  F  F 
FTF  F 
F  F  F  F 
1 

0. 

xxxx  xxxx  xxxx 
xxxx  xxxx  xxxx 
201 
1.0. 

5.  IE-5 
1 

XXXX  xxxx  xxxx 
99999  99999 
0  0  0  0 

F  F  F  F 

XXXX  XXXX  XXXX 

F 


$  PRENAM  POSNAM 
$  RESNAM  WRTNAM 
$  REFSEC  FLUMEM  XXXXXX  XXXXXX 
$  INCTSR  CENINT  BUOYAN 
$  NTINT 

$  STRTIM  DELTIM 

$  EXP  WAV  SPLINE  VARLIN  PACKET 
$  HYPERB  EXPLOS  DOUBDC  VELINP 
$  BUBPUL  SHKBUB  XXXXXX  XXXXXX 
$  NCHARG 
$  HYDPRE 
$  XC  YC  ZC 
$  XS  YS  ZS 
$  JPHIST 
$  PNORM  DETIM 
$  DTHIST 
$  CHGTYP 

$  WEIGHT  SLANT  CHGDEP 
$  NSAVER  NRESET 
$  LOCBEG  LOCRES  LOCWRT  NSTART 
$  FORWRT  STBDA2  ASCWRT  XXXXXX 
$  XV  YV  ZV 
$  DISPLA 
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